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Large wind tunnel facility GROWIAN at Kaiser Wilbelm-Koog, Elbe 
Peninsula. Power: 3000 kW above 12 in/sec Wind speed. 

Rotor diameter: 100.4 m 

Tower height: 100 m 

First operation: Pall 1982, see report on page 



PREFACE 


VI 


At the present time within the energy research program of the 
Federal Minister for Research and Technology, about 50 projects for 
obtaining energy are being supported. Most funds are being expended 
on the development, building and testing of large facilities in the 
MW range. Long term projects Include the introduction of wind energy 
to the power network. 

The seminar was carried out on March 23-24, I 98 I. The purpose of 
this seminar was to exchange information among project leaders in the 
area of wind exploitation. In addition, based on this experience, 
future possible projects and directives were discussed within the frame' 
work of the energy research and technology project. 

During the seminar , first of all , reports were presented about the 
status of large wind facility projects. After this, two special prob- 
lems were addressed: 

- building of large rotor blades 

- control of oscillation in large facilities 

Work on the development, building and testing of large wind energy 
facilities are at present stressing these two tests. According to our 
present knowledge, which was again confirmed in the seminar, the tech- 
nical problems in the design of large rotor blades and the control of 
oscillations in large facilities are certainly solvable. 

Testing of these large facilities is being performed in foreign 
countries in addition to West Germany. The purpose of this is to pro- 
vide operational safety, especially with regard to optimum rotor blades 
and benign oscillation behavior . In addition, initial data is being 
selected about the economics of such facilltiesv The seminar of March 
23-24, 1981 showed clearly that the program on energy research and 
technology of the BMFT for mastering wind energy by large facilities 
will result in the addition of a small percentage- to our energy produc- 
tion. 










The present reports are derived from ongoing research and devel- VII 
opment work and some of the reports will be of a preliminary nature. . 

The BMPT has agreed that this information should be published as a 
special Jul-Spez report of the nuclear research facility Jueilch, 
accessible to the general publie. At the same time, we are making a 
small contribution to the present energy controversy. We should like 
to mention that in the meantime an additional seminar on thfe exploita- 
tion of wind energy was carried out by the energy research' BLE project 
office at KPA. The topic of the seminar which was carried out on May 
18-19, 1981 was ’'meteorological assumptions for the use of wind energy" 
(will also be published as Jul-Spez report no. 132). In addition, the 
results of the projects of the BMPT is being presented during regular 
seminar and status report wind seminars which will take place every 
two years. 

We would like to thank all participants of the seminar for the 
interesting contributions and discussions about tasks during the exploi- 
tation of wind energy. We would like to thank the firm MAN-Neue Techno- 
logie, Muenchen, Which made a substantial contribution to the first 
large wind facility in West Germany. We thank them for their hospital- 
ity. We would also like to thank Herr Dr. R. Cuntze, MAN-Neue Techno- 
logie. We thank Mrs. M.Abt, energy research project manager for 
excellent work in the processing of the manuscripts for the Juel-Spez 
report . 
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ROWIAN ROTOR BLADE 


MANUFACTURING, DEVELOPMENT, BUILDING AND TESTING 





In January 1981, the Growian test rotor blade was finished after 
about 1-1/2 years development and construction time. This report dis- 
cusses the manufacturing technologies and gives the experience made 
during the building and testing of this component. 

BOUNDARY CONDITIONS 

The following boundary conditions were adapted for the Growian 
two-blade rotor, a tower with simple support and the use of laminar 
profiles, assuming a design lifetime of 20 years. 

specific facilities requirements - load assumptions 

- blade geometry 

- lifetime 

—permissible blade tip deformations 
-^permissible eigen frequencies 

.®pseS-flc profile requirements - good shape maintenance 

- low surface roughness 

- close tolerance surface wavlness 

use of tested and cheap manufactur- 
ing technologies 

extensive quality control 

maintenance of close manufacturing 
tolerances 

high requirements for welding seams 

maintenance of the same mass 
distribution for the two rotor blades 

resistant corrosion protection 

DESIGN CONCEPT 

A mixed construction method was selected as the design concept, /3 
which has an untwisted steel box spar as the load supporting element . 

The external blade geometry is fixed by a GPK (glass fiber reinforced 


specific manufacturing 
requirements 
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carbon) sandwich design, which is connected to the spar with fittings. 
This shell design method comes from glider construction and makes it 
possible to satisfy the requirements for the small surface roughness 
and waviness of the laminar profiles. 



blade cross-section 


MANUFACTURING CONCEPT ' 

Within this project, it is intended to build three rotor blades. 
Therefore, manufacturing technologies were selected which would require 
the minimum amount of special equipment such as molds, pressing machines 
etc . This measure, however, requires special skills by the workers in 
order to maintain the close construction tolerances. 

With these assumptions, the following manufacturing concept was 
developed: 

spar - welded construction 

^ ^ complexity 

- reduction of the welding seam 
to a minimum 
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blade skin - GPK (classified with reinforced 

carbon) sandwich design 

- manual lamination method for the 
mold cdvities 

- production of the largest possible 
shell section 

- simple installation 


fittings - simple welding designs 

- variable design in order to equalize 
tolerances between the spar and 
the skin 


By dividing this into the three component groups, it was possible 
to have parallel manufacturing of individual components according to 
the following manufacturing plan. 




MANUFACTURIN3 


P! Ari< 


^i^iClNAL PAG" 

-'NO '.VMiTr pwoVographt 


The spar box was a welded design which represents a reliable 
and tested manufacturing technique with appropriate quality control. 
The welds have to be done very carefully because the continuous 
strength of the basic material is greatly Influenced by them. The 
welding instruction of the spar was adjusted to the manufacturing 
requirements and was divided into several spar segments. By making 
the steel sheets in the form of half shells with a length of up to 
7.5 ni, it is possible to greatly reduce the number of welded seams. 
In this way, it is possible to control the distortions of the compo- 
nents which occurred during welding. 



spar section 

A 40 m long berth was used for assembly where the individual spar 
sections were attached and measured. The spar was welded outside of 
the berth in order to make accessible the connection points from all 
three sides. The construction tolerances and the welding qualities 
were, therefore, maintained using this production method. 

It was difficult to bend the steel sheets into half shells without 
special pressing tools. Deviations can easily occur which are made 
more difficult by the difficult handling of the sheets. The fittings 
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Assembly jig 


for accepting the external skin on the ribs in the end box were also 
made as welding construction and then bolted to the spar. The connec- 
tion surfaces to the external skin are adjustable so that it is possible 
to eq.:alize tolerances between the spar and the skin. The fittings are /J 
assembled in the berth. The external skin is measured later on at six 
measurement points. 

The shells of the blade external skin are made according to the 
proven plastic manufacturing technique. All of the shell components 
are laminated in negative molds. The separation plane of the shells are 
the profile chords and in the internal region after R 27. 5 j these are 
also at the spar trailing edge. In this way one obtains approximately 
equal shell sections which make sure that the hardening time of the 
resin system is not exceeded during lamination with the appropriate use 
of personnel. In addition, the center of the mold has to be accessible 
in the depth direction from the edge of the mold which is another bound- 
ary condition for the dimensions of the mold. 

After an appropriate ti’aining time, it was possible to manufacturer 
the individual shells without many problems. About 10 workers had to 

GhIo.’NAL PAGE 15 

BLACK AND WHiTC PHOrOGRAPH 






El 

IHi 


!»•' Tyro 

1 / 




iW 







ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


lamination of the shells 


simultaneously handle the matrix of each half shell in orde"* to not 
exceed the hardening time. After taking the shell parts out of the 
molds, examinations showed that the permissible surface roughness and 
waviness had to be maintained and was even less than required in the 
external blade region. The required accuracy of the shell weights of 
2 % could not be achieved, there were excesses of up to 


shell set at the blade station R15.7 




no 


rotor blade R50, RIO . 85 

The assembly procedure of the external skin to the spar has 
proven Itself. Matching difficulties between the skin and the fittings 
did not occur. The examined profile sections and blade twist were 
within the permissible tolerances. 

After manufacturing, the rotor blade was transported from Mainz 
to Munich on trucks. During the transport, the blade was supported 
at stations R15 and on self-controlling trailer vehicles and was 

transported on the Autobahn at a speed of about 40 km/hr. The acceler- 
ations during transport were below Ig for the blade. 

RESULTS 


Comparison between the calculated mass distribution and the mea- 
sured ones showed that there was good agreement in the region between 
RIO . 85 and R40. At the blade tip the measured mass distribution deviates 
from the nomln 3 values. This results in excessive shell weight and 
wei; 7 hed weig^its of the spar section between R45-R50. This is because it 
It was not possible to finally divide the spar section after each wall 
thickness Jump because of the manufacturing process. After the install- 
ation of the blade on the test stand, we measured the bending under 
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eigen weight in the flapping direction at the blade tip. We also 
measured the eigen frequencies and eigen modes during test oscilla- 
tion ineasureraents. The deformation at the blade tip is larger and 
the first flapping eigen frequency is 3 * 25 /? lower than calculated. 

The first deflection eigen frequency agrees well with the calculated 
one. 

STRUCTURAL TEST GROWIAN ROTOR BLADE 


Pour static load tests, one test band oscillation test and a dyna- 
mic load test at 50 LW were performed with the Growian rotor blade 
prototype. 

A preliminary evaluation of the measured results shows good agree- 
ment between calculated and measured values.; Deviations from the 
theoretical data occurred at blade station R 27.5 which can be attributed 
to a deviation already found during the manufacturing. At the lower 
shell K27.5-RI5.7 of the profile end box, delamlnatlons. were found dur- 
ing the tests which can be attributed to manufacturing errors. 

The tests were concluded with a fracture test which was performed 
with up to 17055 of the safe load. We did not increase the load further 
because local failure of the structure was established at several points 
of the test specimen. 

TEST RESULTS 

Test stand oscillation test 

The test stand oscillation test was used for determining the eigen 
frequencies and eigen modes of the rotor blade. I6 eigen modes were 
determined using the phase resonance method. The rotor blade was 
excited to perform oscillations using 2 electrical-dynamic exciters 
and the oscillations were measured using 80 'accelerometers distributed 
over the blade. Figure 1 shows the test configuration with the measure- 
ment container of the DPVLR-Goettlngen. 


20 



ORIGIN'Ai: 

BLACK AND WHITE PHOTOGRAPH 



test setup for test stand oscillation test 
Measurement results 


A3 


The measured mass distribution of the blade including the force 
Introduction fittings for load test as well as the DMS (strain gauge) 
cabling were considered over the spar in the determination of the 
theoretical values. 


calculated 

Mode measured value value 

1st flapping eigen frequency 0.923 0.956 

1st deflection eigen frequency 1.27 1.27 

2nd flapping eigen frequency 2.39 2.6^ 

2nd deflection eigen frequency 3-55 3-64 

1st torsion eigen frequency 21.67 26.19 


Load tests 


The load of the test specimen was done using load equipment with 
hydraulic cylinders. The configuration of the cylinders was designed 
so that at blade stations R35 and R27.5 the bending moments were 





test configuration load case 1.3 


simulated as closely as possible. The intersection forces deviated 
by as much as 10?^ from the theoretical values at the other blade sta- 
tion. During tests, the additional strains due to the loading equip- 
ment and the cylinder loads were measured at 500 measurement points. 
Strains due to the eigen weight of the blade were measured only in load 
case 2.5 by rotating the test specimen by l80°. 


The loading of the test specimen is done using two cylinders at 
R^l (80 kN) and R13-6 (^1 kN). 

Measurement results: /15 


The measured strains for blade cross section R 21 .k (undisturbed 
i cross section) have somewhat srr.ailer values than the theoretical ones, 

t For cross section R 27.5, the measured values are higher compared with 

‘ theoretical ones by a factor of 1.56. This can be attributed to a 

1 22 ORIGINAL PAQf 

black and white photograph 




turning experiment 


collapse of the welding seam and excessive edge displacement, which 
is not permissible for the two other rotor blades. The strains mea- 
sured at tbo outer skin are also below the theoretically determined 
ones . 


The blade deformation .t the tip was 
in the flapping direction 1.17^ m 

in the deflection direction O.O 63 m 


For determining the bearing friction moments, in the load case /I 6 
1.3 of the test configuration, at lOOif nominal load the rotor blade 
was rotated with a hydraulic cylinder by ±2.5° at the blade control 
connection, around the blade longitudinal axis. 100 alternating loads 
with a frequency of 0.15 Hz were applied. The measurements at the 


ORIGINAL PAGP 

BLACK AND WHITE PHOTOGRAPH 


23 



beginning and end of the tests showed a slight reduction in the fric^ 
tion moment* The blade displacement moments were 
Mjj = 26.9 kNm in the initial position and 
Mjj = 44.92 kNm in the final position 

This data is within the theoretically predicted range. 

Load case 2.5 

Figure 4 shows the test configuration for load cast 2.5. The 
load is applied to the test specimen by three hydraulic cylinders at 
stations R39 (133 kN) R24.6 (107 kN) and RI 5.6 (126 kN). 

Measurement results : 

The statements for load case 1.3 apply for this case. Here again 
the stresses were increased due to structural deviations at R27.5 by 
a factor of 1,56. 

The blade deformation at the tip was: 
in the flapping direction 0.792 m 
in the deflection direction 1.929 m 

Load case 2.5 

After load case 2 . 5 , a dynamic loading test with 50 alternating 
loads was performed. As a lower load, we selected 20^ and as an upper 
load we selected 100^ of the nominal load of load case 2.5. The mea- 
surement at the beginning and at the end of the test did not show any 
differences in the strains and deformations. The delamination points 
around the end box along the under side of the outer skin between R27.5 
RI 5 . 7» however, increased j otherwise, no damage found . 

Fracture test 

The fracture test was carried out using the test configuration of 
load case 2.5. When 170^ of the case load was reached, the test was 
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terminated. At this time, at station R27.5, we found local instabili- 
ties in the structure, that is, the strains measured in the steel spar 
were no longer linear. In addition, there was a local buckling of the 
blade skin at R27.5 along the blade underside. A 150 mm hole in the 
skin was the reason for this , which had to be di’illed in order to exam- 
ine measurement points over the spar at the beginning of the test. 
Additional damage in the skin was not found (photograph page l8). 

Measurement results: 


The stresses in the spar and in the outer skin achieved for 170 % 
of the safe load are shown in Figures 10-12. 

The following was measured at the blade tip 
flapping direction 3.235 m 

deflection direction 0.888 m measured 

After unloading of the specimen, a residual deformation of about 
0.400 ra was found at the blade tip. 

CONCLUSIONS 

A rough evaluation Of test results shows good agreement with calcu- 
lated results. Since the rotor blade will be subjected up to 108 load 
alterations during operatioh, it was dimensioned for continuous strength . 
However, we can thereby explain the large separation between the safe 
load and the static fracture load (J'v2). In addition, tests have 
Shown that the blade skin made of glass fiber plastic will follow the 
spar deformations without tearing and the strain of the skin are by 
a factor of 10 smaller than the measured fracture strain of the laminate. 
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General remarks ; 

The Voith-Qetrlehe KQ Heldenhelm, West Germany, Division for 
Aviation Technology already in 197^ planned the development, manufac- 
turing and testing of a wind energy aonverter of medium power and slse 
for use as a fuel saver for Island operation and network operation. 

It was primarily planned for rosions with low infrastructure in the 
third world. 

We consider an axial flow equal pressure turb.lne with a gondola 
rotor in the leeward side. It had 2 PVW I'otor blades with the shell 
design which are positioned over a wind flag (no side wheel gondola 
deflection device), Also, we considered a steel tube mast supported 
with guy wires . 

The star-shaped reinforced Goncrete foundation was not anchored 
to the ground . 

As a location for the test, we used the terrain at (Huetter) in 
Stoetten in the Swabian Alps, in the Ooeppingen district . 

Professor Hlrich Huetter Stuttgart/Kirchhelm (PWE) advised us 
about the aerodynamic design and the basic concept. 

The firm MBS Munich was considered for the development and supply- 
ing of the high performance rotor blades after a long search. 

is as follows and was maintained over the 

rotor rpm 37 l/nvin 
rotor performance 31^ 
synchroneous generator 350 KVA 


The main data of the WEC 
years ;; 

rot or blade diameter 52 
nominal speed 8, 5 
mast height 3Q m 


2S 



1. Goal of ET 4104 A (facility) ; 
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According to the 52 m diameter size and the assumed DFVLR load 
assumption (see Figure 1) and a century gust of 60 m/sec, we consider 
the following; 

- completely automatic facility, including startup and shutdown 

- self-saving (fail-safe) machine Installation 

- a synchronous generator which switches in and out of a three- 
phase current network automatically (20 KV, 50 Hz) 


The following was expected: 

- long lifetime 

- high availability 

- operational safety 

- easy asembly (if possible without excessively large cranes) 

- easy maintenance with long intervals 

- small weight with high strength 

- competitive price 

- national manufacturing components for the third world 

- environmentally sound with regard to appearance and noise 

- movability of the entire isntallation (other test field) 

- safe against unauthorized removal 

/2l 

2 . Goal of ET 43^3A (rotor blade development) ;' 

- the recommendations of Prof. Huetter included the following 
for the Voith facility with a fast running coefficient of 


^ ^ = 11-85 (-e Figure la) 

or for c -optimum =0.45 according to an extremely high fast running 
number of - 16 

(see Figure 2) and very slender rotor blades resulted with an area 
loading of 

316000/52^ XT - 149 W/m^ 


- we consider the fact that the masses of the completed gondola 
should be made as small as possible in order to not exceed the 
total eigen frequency of the WEC with an operating frequency 
of 37 X 2 = 1.211 Hz 

- the rotor blade mass was limited to about 1000 kg per blade, 
and the slender profile family of the blade cross sections were 
specified by Prof. Huetter. 

These parameters led to an extremely light hollow blade design 
with support foam inserts and transverse bolt attachment. 

The initial plan to divide the blade into two lengths and to 
flange them together as required had to be given up because of oscill- 
ation problems. 

/22 

3 . Goal of the ET 4344 A (telemetry) ; 

The behavior of the two slender rotor blades has to be monitored 
especially when penetrating the wake of the mast and excessively high 
blade deflections have to be controlled by turning off the facility. 

The transfer of the measured current of strain gauge strips and 
accelerometers cannot be done using slip rings. 

. Solution of the problems of ET ^[10^ A (facility) ; 

Reference Hoffman 

k .1 Structure 

Prom a number of possibilities (see Figure 3) we finally decided 
for the prototype shown in Figure 4: 

- pendulum hub according to Prof .Huetter (semi-cardan suspension), 
steel sheet welded construction with the possibility for self- 
controlled blade angle adjustment (in order to equalize shear 

28 - . ^ ' ' 




flow) and to 1) look it (see Figure 5). 

- pendulums bearing pair from prestressed double conical roller 
bearings (over dimensioned because of the fact that the period- 
ic rotational motions were only small) using continuous grease 
lubrications 

- rotor blade bearings, one prestressed and coiranerclal crossed 
roller bearing each with a measured friction value under full 
load of about 0. 003 with continuous grease lubrication and 
double seals against water entry 

- adjustable control rpds with commercial spherical joint -head and 

continuous grease lubrication I 

- hub covering with trim disks as well as solar shells for mea- 
surement transducers mounted on the outer surface (ET 

telemetryl ^^^ ^ ^ ^ ^ ^ ^ ^ ^ , 

■ ' ■ ■ r 

- the completed pendulum hub is connected to the rotor hollow shaft 1 

using a Hirth gear (see Figure 6)t The shaft is suppox’*ted at ' 

two places with a throat bearing which is developed as a cylin- 
der full roller bearing which can be displaced in the axial 
direction and a prestressed conical roller bearing for the oper- 
ational axis movements which occur during operation. Both bear- 
ings have oiled lubrication 

-the conical wheel set 70/17 teeth or i - ^.12 has a reducing effect 
on the aerodynamic thrust because of its axial tbrust, that is, 
the conical roller bearings are unloaded in this way. 

'■ I 

.... - : • ' ■ : ' ' ^ ■ ■ ■ . ■ , . ' 

- accepts the rotor blade adjustment inside /23 | 

(see Figure 7) • ^ I 

, ■ ■ ■ . . ' ■ ■ I 

- a differential system fol^ rotor blade adjustment is also provided I 

by a hydrostatic locking or damping of the piston rod. This oould | 

be important if blade angle deflections occur caused by pulsating | 

■ ■■ V 

^ 


aerodynamic forces. 


a screw spring also provides blade adjustment in the "flag” 
direction which is integrated with the adjustment mechanics, 
that is, the blades then have no lift, the WEG facility will 
turn off itself in emergencies . ’ 

the position information of the control rod or rotor blades 
which is important for the controller is done with a potentio- 
meter (see Figure 8). 

the gondola hydraulics (see Figure 9) is supplied by a wing 
cell pump with 0 stroke control . The electrical motor obtains 
it current with slip rings (see Figure 9)» 

the oil supply and removal for the differential system Inside 
the rotor shaft is done with an electrical hydraulic control 
slide of our own design (see Figure 10). 

the heatable oil container is installed in the gondola housing. 
All of the hydraulic lines are, therefore, short . The pump 
supply is done with a natural slope. 


the hydraulic unit is attached outside on the housing and is 
easily accessible (see Figure 11 ) . 

the ’gondola housing is flanged to a rotatable bearing (see Figur 
12) . The worm gear wheel is connected with a fixed mast, the 
gondola housing is centered using a conventional crossed roller 
bearing using the mast. 

A three stage unlimited worm gear supported in the gondola hous- 
ing gives a gear reduction ratio of i - ^ - 16.67 

3 

After the worm gear shaft there is a combined commercial spur 
wheel worm gear which has 1 - 21^.8. The drive motor has 1^40 
rpm . For a total i„^„ = 21^.8 x 16.67 = 3*581, in other words, 

■■ ■ ■ , ■ ■ ■ ' " .;:g,©b .. ■ . . , - . 



the gondola would have a theoretical rpm pf 0#4 rpm. Between 
the worm gear shaft i and the spur gear worm gear there is an 
electrically controlled lamella coupling which is closed when 
power is available. This provides absolute self-control for the 
entire gondola operating range. 

When there is power failure, erroneous gondola positions, etc., 
or oblique loads for the unfavorable case of rotor blades in the 
flag position, a gondola torque could be created for very high wind 
speeds which would rotate the gondola out of the oblique incident 
wind position. This is possible because the electromagnetic 
coupling would not be switched in in this case. 

three electrically controlled hydraulic pistons are distributed 
along the circumference of the gondola and mast flange for oper- 
ating the plier breaks . The brakes are turned on when the 
facility is operated and are only opened during a positioning 
c ommand or in the event of power failure. 

in order to Introduce current for the two gondola motors for the 
control voltages , receptacles, etc., a slip ring transmission 
is required (see Figure 13) which is below the gondola rotating 
gear arid housing, arid which belongs to the fixed mast . 

The measured currents are not run through slip rings but are 
directed directly to the controller, recording devices, etc., 
with a 20 conductor flat cable. This winding sleeve makes it 
possible to have ±5 gondola rotations before the connector con- 
nection would be opened. This means that already after about 
three gondola rotations a restoration command will be given and 
this would result in a turning back under favorable conditions . 

two pedestal planes are provided for assembly and maintenance 
work in gondola 2 (see Figure l4). A safety conductor goes from 
level 5m to the gondola plane. 


a v:lnjg wheel anemometer and a combined wind flag are located 
above the gondola bousing as well as the required lightning 
protection. 

additional lightning slip rings are provided for protecting the 
main bearings between the rotor blades and the rotor shaft, 
between the rotor shaft and the gondola housing and between 
the gondola housing and the mast . 


The mast is 'connected several times with iron reinforcement of 
the concrete foundation and the conventional grounding wires 
are used. 

The long tube shaft which goes downward (see Figure 15) between 
the conical shaft and the spur wheel shaft is supported with 
two V-belt connections and can be displaced in tbfe ’longitudinal 
direction. In the central region (see Figure 15) it is sup- 
ported with a pendulum ball bearing with respect to the mast in 
order to hold the critical bending rotation rate above the 
operating rate of 152.7 rpm. 

The spur wheel gear, the normal Volth design 1 = 9 •82, is coup- 
led through a jointed shaft (see Figure l6) with the 350 KVA 
synchronous generator. 

2 ring! tensioning rings between the generator shaft end and the 
coupling flange transmit the maximum torque = 2.5 x operating 
torque at about 300 kW and n = 1500 rpm. In this way, possible 
short Impacts from the WEC transmission are blocked. 

A fail-safe plier brake below the spare wheel gear output shaft 
(see Figure l6) provides for fixing of the entire shaft assembly 
Including the rotor and the generator in the event of danger or 
during repairs . 

the generator has a tachomacbine for monitoring rpm ^ two addi- 
tional toothed disks with electrical transducers measure the 


there is a large assembly door at the base of the mast in the 
cylindrical cover which allows one to remove the generator or 
fc he spur wheel gear or the brake. Three small inspection holes 
are also provided. 

in order not to have to separate the electrical force and con- 
trol lines inside the mast base during tipping of the mast, a 
special clamping box is attached at the foundation flange. 

there is a ladder with protection provided for climbing the mast 
up to the upper gondola assembly platform. 

the total mast including the two rotor blades is 38,800 kg. The 
overall eg of the facility is at the height of about -20 m. 



load assumptions 
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FIGURE 2 (CONTINUED) 

KEY: 1 — closed after (Illegible); 2 — tungsten granulate; 3 — section; 

^ — detail; 5 — cross-section of blade; 6 — Voith axes; 7 — profile chord; 
8 — sandwich region; 9 — mode separation plane; 10— Voith axes: axis 
system for installation of blade to rotor bead; 11 — manufactured 
according to FV; 12 — (illegible); 13 — all laminates laminated with 
mixture of (illegible); 1^< — glued with mixtures from (illegible); 

15 — left (Illegible) ... (others illegible)... 
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FIGURE 3 (continued) 

ABBREVIATIONS: ( 6/1 ) crossed rope operation 

( 6/2 ) Jointed shaft 
( 8/2 ) two stage 
(12/2 ) telescopic tube 
^2/^1 ) reinforced 

(12/5 ) elllptlcal/cyllndrlcal tube 

(13/12) guy ropes/rods 

(15/10)+(15/12 ) hydraulic cylinder 

(18/12) stable mast/subcrltlcal opertlon 

( 18 / 13 ) mast frequency below excitation frequency 

PATENT APPLICATIONS: ( 1/1 )-( 1/14) blade flange 

( H/6 )-( 4/1 i|) complete card, suspension 
( 6/1 ) rope operation 

(14/3 )-(14/12) buckled mast 
(15/3 )-(15/10) pulling device 
(15/12)-(15/13) pulling device 
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Section C-D 
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FIGURE (continued) 

Key: 1 — device for attaching transducers (for wind direction and 

wind speed); 2 — platform not shown on this drawing; 2 — grid support 
only at the Indicated points; after this configuration of grid 
support rods; 3 — view Z; 4 — section; 5 — view; 6 — rail complete; 

7 — grid complete; 8 — grid configuration; 9 — design of upper pedestal 
with rail ... (others Illegible),.. 
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FTilllRE 15 (continued) 
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Key: 1 — detail; 2 — detail; 3”-detall; — detail; 5 — upper shaft 

connection; 6 — middle shaft connection; 7 — lower shaft connection; 
8— markings In the plane; 9— after installation (illegible) 
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FIGURE 19 (continued) 
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Key: 1 — alarm indication: 

power out 

phase out 

frequency 

control voltage 

control voltage T5 

excess current 

excess temperature 

excess RPM 

reverse power 

emergency shut down 

excess current 

fracture warning 

fracture triggering 

L measurement absent 

safety monitoring 

erroneous current FI 

excess current 

computer excess current 

DC tower head position motor 

excess current control oil pump 

excess current heating oil 

excess current heating anemometer 

excess current power supply 

UTW 

excess current coupling head position motor 
UTW magnetic valve 

excess current magnetic ventilation main bearing 

excess current brake 220V fuse 

excess current brake battery supply 

brake block 

mechanical disturbance 

2 — main distribution; 3 — input distribution; 4 — later; 5 — computer; 

6 — fixing brake; 7 — coupling head position motor; 8 — immersion body 
with thermostat for oil; 9 — wind measurement device with heater 
power supply; 10 — magnetic valve UTW; 11 — power supply gondola; 

12 — control pump for oil; 13 — head positioning motor; IH — magnetic 
valve main bearing brake; 15 — no. 24 measurement; l6 — from main distri- 
bution; 17 — power operation 
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THE CONCEPTUAL DESIGN OF THE 
HIGH SPEED RATIO WIND ROTOR FOR THE 
62 M DLAMETER HORIZONTAL AXIS VOTTH WIND ENERGY CONVERTER 


W. Weber 

Volth Gctriebe KG, Federal Republic of Germany 


Summary 


The Volth wind turbine la a development of a borleontal axis turbine with two blades 
fully constructed In composite mateiialsL It Is designed either to supply electric power with 
a 265 kW generator to an electric grid or to be an Independent electricity source for specific 
tasks In highly developed countries and In remote areas. For this reason the conceptual 
philosophy of this wind turbine goes In a slightly different direction to other developments of 
wind energy converters actually designed and In operation, a situation, which may be explained 
In the first chapter. 


60 


ORIGINAL PAGE 13 
Of POOR QUALITY 


1. Specific power instaildtlon 

According to the measuretwiUs ui medo wind velocity distribution in Germany in 
/!/ as sho*<n in Fig. 1, there is an evident difference between the coastal 
region with hi^ mean wind velocity and Inland areas with lower wind areas. The 
Voith wirtd turbine is designed specially for inland areas with a mean wind velo- 
city of 4-5 m/sec. This wind distribution profile is not restricted to Germany. 
It is essentially also found in other countries. Therefore the Volth wind tur- 
bine is usable worldwide mainly in inland areas of most continental couitries. 
For this region the right specific power installation should be analyzed. 

The nujor distribution to the solution of this problem and the resulting concep- 
tion of the wind rotor shape we obtained from U. Hotter. 

According to his studies and the calculation based on different wind power 
measurements on several weather stations in Germany it is possible to describe 
the variation of the yearly wind power output from a specific power installation 
factor 

Pr 

IT = 

A 

with Pp _ f 3 ted power 

and A = rotor disk area 

The results of calculations in /2/ demonstrate that a clear dependence exists of 
the zero operation time/year ori this specific power installation factor T . 

Fig. 2 gives a relative feature of this dependence on the mean wind velocity 

^10 measured at a heigth of 10 m. ^his figure shows that there is a need of 
low power installation factor "S , if the wind turbine Installatia*) point lies in 
low wind regions of s coentry in order to get a fairly continuous delivery of 
wind energy. The Voith WEC 52/265.8, therefore with a factor of » « 125 •/«>2 
has about the same operation time/year as the German Project GROKIAN with r » 

300 w/m2 on an installation point in the costal area of the North Sea. The go- 


I 

neral trend In many wind power deve'opments is towards high factors and In- 
stallation points with high nean wind velocities. 
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2. The choice of the speed ratio 

The idea of the designed concept is to test a high-speed wind turbine with an 
efficient blade performance. According to the extensive studies of U. Mutter for 
the Volth wind turbine at the testfleld In Stotten - the same place where the 
former • 34 were in operation from 1956 till 1964 - the maximum extractable wind 
power In Stdtten and on comparable sites In central Germany is possible at a re- 
presentative mean wind velocity of V^q - 5,3 b,/s at 10 m height as seen In 
Fig. 3. 

For this wind speed the portion of the total energy output reaches a maximum, if 
the frequency distribution of the wind velocity after /!/ represented by diffe- 
rent mean wind velocites is taken into account. Even thou^ the low wind velo- 
cities reach the hipest frequency distribution for the two stations selected as 
representative for the general - in Fig. 5, they do not get the most position of 
the possible energy output/year, the same situation as with the low frequency 
distribution of the high wind velocity of = 9 m/s and of V^q = 11 m/s. 

Based on the wind velocity of =6,3 m/s and a constant angular velocity 
of o = 3,9 1/sec (n = 37 rev/mm) the Voith wind energy converter has a 
designed tip speed ratio of 

u> R 

X j = = 16 

''10 

with the rotor tip radius R s 26 m 
and reaches the rated power at 

^TR = 12 (ViQ = 8.7 m/s) 

a ratio, which is an upper limit in respect to other wind turbines as mentioned 

in Fig. 4. 
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The consequence of this layout for the construction of the Khole rotor gives 
some adv^tages and on the other hand sotne problems, «^lch are outlined in the 
following context. 

/50 

3. The effects of the high tip speed ratio 
a) The blade shape 

Following the classical way of determining the optimim performance of wind tur- 
bines as published in /3/ and /A/, the designed speed ratio requires a certain 
contour form of the blade. The result for the Volth rotor is a wing blade which 
is extremly slender compared to other blade configurations as shown in Fig. 5 . 


The main geometric parameters are 


r2 

the aspect ratio of one blade flR = 

^T 

the taper ratio V - — * 0,27 

t 

max 

and the solidity factor of the rotor & = 

A 


with 



blade area 

maximum chord length 

chord length at blade tip 


• * 
I 

i 

I 

= 39 


0,01^ 


The blade twist has a range of 


A Og - 70 


in the relevant blade section from 6 m up to 26 m radius. 

One deslgied goal for this rotor is to operate the wind turbine in a 1 p-upper* 
critical working range. In order to realize the required stiffness and mass 
distribution for these frequency and response calculations, the construction of 
the blades is only possible with the use of advanced composite technology. 
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construction technique was solved in close collaboration with the helicopter 
division of the Gernen aircraft manufacturers HB8, of whom we shall hear more in 
this meeting. 


b) The power coefficient 

Going back to the definition of the power i:oefficient OPIGINAL PAGc 13 

OF POOR QU.V.'T'Y 

P 

Cp * 

?/2Vjo A 

with ^ * air density and the former definition of Vjq and A, the power P can 
have kjider constant conditions of the other terms, different definitions: 

In idealized conditions, which means without drag, the power coeficient is in- 
dependent of the lift/drag ratio (L/D) of the profiles; 

Sid 4 f (L/D) 

In Fig. 6 this idealized power coefficient the Voith rotor reaches 

its best values at a mean wind speed 

5,5 m/s fc 5^5 ^5 

But, in general there is a strong Influence of the profile drag or the power 
coefficient as it is published in /5/ and see^ '’ig. 7. 

with the increasing tip speed ratio a increasing L/D ratio is required to get a 
good performance of the rotor. The influence of L/D ratio is usually expressed 
by an efficiency factor rj 


For structural reasons with Increasing bending and torsional moments and there- 
fore necessary growing inertia moments towards the blade root it is not possible 
to have a constant L/D ratio over the rotorblade at all radii. Thus, the Voith 
rotor blade with the given thickness distribution has a profile variation going 
from the thick profile ODE 642 in the inner parts to the thin profile WH 301 in 
the outer parts of the blades, as is demonstrate<^ in Fig. 8. 

The L/D-ratio distribution is similar to the measured values of the Wortmann 
FX-w profile family going from about L/D = 30 in the 6 m section up to L/D * 150 
at the blade tip. The calculation of L/D efficiency factor is given in Fig. 9. 

For the representative radius of r = 0,75 R the L/D efficiency decrease not very 
much in the case of optiinjm pitch angle with A = 0", between 8 a 18 

whicii indicates a good utilisation of the profiles by thi> desigied tip speed 
ratio of A.J s 16. 

From Fig. 7 it can also be seen that the influence of the nuntoer z of blades to 

the power coefficient Cp of the rotor decreases with increasing tip speed 
ratio and therefore does not tremendously effect the Voith rotor with a deslged 
tip speed ratio of Xj - 16. 

The effect of the nuirber of blades z can be expressed in an analogous matter to 
the L/D influence, which means, to calculate an efficiency factor i| 
factor is consequently a function of the tip speed ratio Ay 

r|, , f ( Ay) 

Including all effects, the integrated mean power coefficient of the whole rotor 
becomes finally 
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'^pg * *1 2 J Cp • n L/D dr 

R 0 
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The results are given in Fig. 10 . The Volth rotor will have a range of operation 
between a cut-in wind velocity of _ 4^5 ,4/5 a cut-out wind velocity of 
V^o . 25 m/s. 

The rated power will be reached at _ 0^7 ,4/5 wind velocity. An im- 
pression of the rotor loading conditions give the nean lift coefficient c. of 
the rotor blades. It can be seen that with decreasing tip speed ratio the lift 

coefficient c^ of the blades increases up to the rated p^-*r At the rated 
power and with further increasing wind velocity the mean lift coefficient goes 
back again, while there is no need to use the whole capable wind energy. The 
blade pitch angle should be turned toward the lower angle of attack. 

4 . The consequences of the high performance blades on the construction and the 
aeroelastic behavior 

To conclude the presentation of the designed Voith rotor concept in respect of 
power optimatlon it is necessary to give a survey of the consequences on 
structural and aeroelastic behavior of the blades. The rotor as presented here, 
has in some aspects considerable advantages. 

a) The slender blades with low blade area give a considerable mass reduction 
compared to bigger blades with low ^deslgied tip speed ratio as shown in Fig. 5 . 
Therefore the alternating in-plane bending moments and forces due to the 
gravitational loads, the gyrocopic moments and other mass related 
characteristics are reduced. 

b) In the loading cases of hlc^ wind speed with no rotating blades the blade 
root moments are in the same order as the moments at rated speed and have r»o 
greater problems then the moments at rated speed. 

c) The hlgb rotation speed of the rotor shaft reduce*, th? needed gear ratio to 
the generator and reduces consequently the loss of efficiency due to this gear 
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d) The torsional natural frequency of these slender blades are high, «<hlch :ias a 
positive effect on the aeroelastic stability behavior. 

On the other hand the chosen rotor concept has also realisation difficulties, 
which are to be solved. 

a) The L/D ratio requires an aerodynanvic performance which results in a 
necessity for composite materials with surface conditions of the blades compar- 
able to the wing configurations of modern plastic competition gliders. 


b) The required stiffness and mass distribution for this slender blade is rnly 
possible to realize through the exter.slve use of carbon fibre composite mate- 
rial. Both the components and the handling process are relatively expensive. But 
on the other hand, a mass production of rotor blades in this configuration will 
bring the cost to a lower level. 


The presented development of the Voith rotor will expect some improvements for 
wind power utilisation. A final evaluation will give the operation of the wind 
turbine and the related study of the cost effect! vness of the concept. 
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Fig. 2: Zero operation time/year for different power calculations 
cn 25 weather stations throughout Germany 



Fig. 3. Freguency distribution of mean wind velocity and annual 
energy output distribution per year of two selected 
weather stations in Germany 
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Fig. 4: Comparison of trends in the design of 
power installation on different wind 
turbines 


Fig. 5: The blade geometry of the Voith 

WFC 52/26*'. R wind turbine in relation 
j ’he <1F t40n-1 rotor blade 
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Fig. 3: The power coefficient of the WEC 52/265.8 rotor concept in idealized 
condition without friction 
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Fig. 7: The influence of the lift/drag ratio L/D, and the 
number of blades z on the power coefficient 
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Fig. 8: The distribution of the thickness d/t over the rotor 
blade and the associated L/D rat’o as a function of 
thickness ratio d/t 
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Fig. 9: The efficiency of the L/D 
ratio by different blade 
pitch angles A Be at the 
• erresentative blade radius 
r = 7B " P. 
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1. Introduction 
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Within our research project ET JllO^ A with the title "Development, 
building and testing of a wind energy conyerter WEC 52/265.8", it 
very soon became clear that the design of the rotor and its blades 
was very important . 

New technologies and design methods for bulldilng the rotor blades. 
Rotor blades with a composite fiber design were developed and CPK 
(carbon fiber reinforced plastic) with T 300 fibers were used for the 
primary structure. 

Extensive calculations on strength, stiffness and structural damp- 
ing gave good results using this method. 

However, not all of the required input data were available for 
the calculations so that estimates and assumptions bad to be made. 
Especially coupling influences of the rotor blades on the mast were 
difficult to establish for dynamic behavior. 

For safety reasons, therefore, experiments were made with singled 
rotor blades in test band experiments and also rotors were tested under 
loads in a complete assembly. In this way, one wishes to verify the 
assumed data or to make corrections for the design. 763 

Measurements of rotating rotor blades of WEC ’s are not very con- 
ventional because measured signals have to be transmitted for data 
recording to a ground station on three rotating planes, that is, the 
rotor blade plane, the rotor shaft plane and the gondola rotating plane. 

Analysis shows that the dynamic behavior of the wind turbine rotor 
could be sufficiently covered using l6 measurement points so that the 
solut ion of the measurement task would best be satisfied by using a 
16 channel telemetry facility. 
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2 . Measurement task 


For the oscillation measurements of the two blade wind turbine 
rotor, l6 transducers have to be installed on each rotor blade. In 
order to record measurement signals for deflection and flapping motions 
at the blade tips, one accelerometer transducer each is required. 

For the other measurement points, strain gauges have to be installed. 
Accordingly, l6 analog measurement points are to be measured at about 
the same time. Remote control for measurement point switching from 
one rotor blade to the other will allow on-off switching of the telemetry 
measurement installation from the power supply. 

The measurement point distribution on both rotor blades was esta- 
blished after the conclusion of the EDP calculations of rotor blade 
optimizations and when plotter results were available for the rotor 
blade eigen oscillation modes. A measurement point distribution accord- 
ing to Figure 1 was used. /64 

Because of the geographically exposed location of the wind energy 
converter, and the fact that the rotor blade was partially made of GFK 
(carbon reinforced plastic) fiber material and, therefore, could conduct 
electricity. The telemetry Installation has to be protected against 
lightning. 

In order to provide a high degree of data transmission safety, the 
measurement signals have to be amplified in the rotor head, they have 
to be digitalized and a pulsed code is made for wireless measurement 
value transmission (PCM code). The measurement data which can be trans- 
mitted are to be transmitted to the ground reception station through a 
PCM link. There the measurement data are decoded again, first of all, 
using a line recorder they are optically recorded. In addition, they 
are stored in parallel using an analog tape machine for later evaluation 
with the EDP installation. 

The measurement task was performed in two tests: 



2 . 1 Test stand experiment 


The calibration of the measurement transducers is an important 
part of the test. For this purpose, both rotor blades are statically 
loaded in a test device in the deflection and flapping direction. The 
load is introduced in steps to the rotor blade through a load device 
and the deflections which occur have to be determined. 


The entire telemetry installation in the original version has to 
be tested during measurement value calibration. 

The research task ET 43^3 A for the rotor blade development which 
was performed as a subcontract by MBB, also contains model analysis of 
the finished rotor blade. 

Using the telemetry installation, data recorded in parallel and 
stored on tape are then evaluated later on as a test example by the 
oscillation analysis devices of the firm MBB. 

This evaluation at the same time is used to test the adaption capa 
city of the various devices in order to later on be able to evaluate 
the measurement data of MBB recorded during free field tests by Volth. 

2 . 2 Fr e e field tests 

The dynamic overall behavior of the two blade wind turbine rotor 
where the Incident flow comes from the leeward direction has to be 
determined over the entire operating range . 

3. Components of the telemetry measurement installation 

The following description will give some information about the 
operation and structure of the components of the telemetry measurement 
installation shown in Figure 2. As already mentioned, the rotor blades 
are equipped with strain gauges and accelerometer transducers for 
recording measurements (see Figure 1). 



Both measurement transducer types are switched in a full bridge /66 
circuit and can be treated electrically in the same way. The measure- 
ment transducers are the main current users for the power supply and, 
therefore, their internal resistance is made small. The measurement 
transducers are connected with lightning protection units through 
simply twisted and shielded cables installed in the rotor head. 

In the case of lightning impact in the rotor blades. Induction 
of voltages in the shielded measurement and supply voltage cables can- 
not be avoided. The lightning protection unit is, therefore, dimensioned 
so that induced voltages of up to 40 kV can be removed through combined 
self-extinguishing spark gap and diode networks against the mast. 

In order to switch all the measurement points including their 
power supplies from rotor blade A to B, a mechanical relay system 
protected against vibration is used. 

Switching through a remote control radio device, which is planned 
for turning the circuits on and off of the telemetry measurement 
installation, is operated from the ground station. 

The power supply of the measurement transducers, preamplifiers 
and telemetry units from the measurement station was selected from 
various possibilities. This is located in the rotor head; However, 
because of possible disturbing Influences, slip rings had to be avoided. 
The installation of radial and axial induetlbly. coupled current trans- 
mitters was investigated. Both would have been possible but the /67 

latter version was not possible due to the design. Therefore, we con- 
sidered an autonomous current production using 5 silicon tsolar generator 
panels which are installed ahead of the rotor. Discharge protection 
diodes prevent Influencing of the solar generators . 

A voltage precontroller provides continuous charging of the niekel- 
cadmium buffer battery. For unfavorable weather conditions and a full 
battery, an operating time of at least 24 hours is achieved. 
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The central part of the telemetry measurement Installation con- 
tains the current -supply unit, conversion of the battery voltage and 
this supplies the measurement transducers, preamplifiers, PCM encoder 
and HP transmitter. 

The individual bridge voltages of the measurement transducers are 
amplified in the preamplifiers and are switched to a single channel 
using the time multiplex method and then direeted to the PCM encoder. 

Here there is analog to digital conversion, the measured values are 
represented as digital words, parallel to serial conversion takes place 
and the code is converted. The produced data stream is then frequency 
modulated and is superimposed in the HP transmitter onto a carrier fre- 
quency of a radial link according to the requirements of the German 
post office. The reverse data conversion occurs in the reproduction 
of electronics within the ground reception station. The measurement /68 
data obtained in this way are either directed to the line recorder for 
display or are stored on magnetic tape in parallel. In the EDP facility 
the PCM decoder also has an Interface for the computer for further pro- 
cessing. 

4 . Partial assembly of the- telemetry components 

The telemetry components can be seen in both order housing parts 
in Pigures 3 and 4. The housing parts are arranged diagonally between 
the rotor blades. Although the connectors were designed according to 
MIL-SPEC because the high requirements against contact certainty and 
water tighteners are required. Plgure 5 in the lower part shows the 
control console of the WEC with the telemetry encoder. 

5. Present state of assembly 


The planning of the rotor blade oscillation measurement was 
designed originally for an early production schedule of the entire 
facility. 

This means that partial assembly of the telemetry components was 
already complete during the first quarter of 1980. 


At the present time the individual parts ’.of the rotor head are 
manufactured and are being assembled. /69 

In this connection, we now are carrying out the calibration of 
the strain gauge measurement points l6 which can be seen in Figure 6 on 
a thrust yoke on the rotor shaft. This figure also shows the wiring 
work. 

6 . Planned schedule 


Test band experiment 2.1 will be carried out immediately after 
completion of the rotor blade. The firm MBB has Indicated that this 
would occur at the beginning of the second quarter of 1981 and MBB 
according'to 2.1 is responsible for the development of the rotor blades. 
The outdoors test will probably be able to start during the third 
quarter of 1981. 
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FIGURE 3 
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PIGUPiE 6 (continued): 

Key: 1 — telemetry configuration; 2-rotor blade A; 3^ — measurement 

point; 4 — measurement point; 5 — rotor shaft; 6»-lightning 
protection; 7 — bousing 2; 8^ — PCM telemetry transmitter; 

9— transmission antenna; 10 — solar cell panel; 11 — solar 
cell panel; 12— controller; 13 — battery; l4 — on/off relay | 
switching system; 15— reception antenna; 16 — remote control; 
17^^houslng 1; 18— lightning protection; 20 — measurement 
point 16B; 21— measurement point 1. . .15B; 22— rotor blade B 
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ROTOR BLADE DEVELOPMENT ET 43^3 A /77 

Relnhard Krautwald 778 

SUMMARY 

This report describes the development and manufacturing activity 
performed during the rotor blade development ET 4343 A for wind energy 
facility Stoetten i40707b task during the period between May 1, 1979 to 
February 28, 1981. 

With consideration of the basic data of the overall design such 
as geometry, mass and stiffness distribution, a rotor blade using the 
fiber compound construction method was constructed. We demonstrated 
that the planned goals were achieved using static and dynamic tests on 
a 10. 2 m prototype of the internal blade region. 

In addition, a new kind of molding and working concept is des--" 
cribed which makes it possible to jnanufacture components of this size 
and this is compared with conventional manufacturing technologies with 
a small manufacturing complexity. I 

The report ends with a discussion about the manufacture of the | 

half wing shells of the first 26 m long blade. I 



Status ; February 28, I98I 
Messer scbmitt-Bolkow-Blohm GMBH 
Helicopters and Traffic, Munich 
Project director: Ralf-Thilo Schultz 
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I, Goals 

With consideration of specified design parameters (essentially 
blade geometry and blade stiffness), two rotor blades for the wind 
energy facility "Stoetten 40707 b”are to be developed using the com- 
posite fiber construction method and these are to be manufactured. 

II. Construction 

11. 1 Blade design 

Since the blade mass, mass distribution and stiffness distribution 
are related directly, it is only possible to achieve an approximation 
to the final goal using blade cross section calculations and bar model 
calculations and several iterations. After achieving the required 
values, the blade data are compared with the requirements for the 
total facility and are then verified. 

With consideration -of the loads which occurred during operation, 
the connection elements (expanding bolts and transverse bolts) were 
dimensioned and the blade connection stiffnesses were tested. Ther- 
mally, the blade connection was designed so that the thermal expansion 
Coefriclent of the connection cross section agress with that of the 
steel blade bearing so that there will be no friction between the com- 
ponents . 

11. 2 Construction design 

Based on the high stiffness requirements with the lowest possible 
blade mass, a ’'embedded" blade (foam core) was considered as the design 
principle because the spar construction method for large wings has a 
less favorable ratio of stiffness and mass. In addition, for a full 
foam blade, the fitting and control problems do not exist which would 
occur for spar fittings of small height in the external blade region. 
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The bending stiffness in the flapping deflection directions is /83 
provided for by unidirectional carbon fiber tissue paths (warp, weft 
ratio 9:1) in the blade longitudinal direction. The torsion stiffness 
is provided by glass fiber tissue paths in the ±45° direction with 
respect to the blade longitudinal axis. 

Because of the fact that the blade mold cannot be closed (III. 1.2) 
in each case a half a wing is built up and then both blade halves are 
glued together. The torsion connection is provided by gluing on tor- 
sion shells to the nose edge and the end edge. The foam core acts as a 
local support -structure (bulges) for the shells as well as for accept- 
ing shear forces. Force introduction into the rotor blade bearing is 
done by means of the blade connection already discussed during the 
DNW program (German Holland Wind Tunnel) using transverse bolts and 
expansion screws. Transverse forces are introduced into the bearing 
through a base plate glued to the blade connection using shear fittings. 

The lightning protection system consists of a profile aluminum 

O 

strand (cross section area 50 mm ) which , is Integrated with the blade 
structure. The trailing edge protection consists of a 1 mm thick and 
50 mm wide aluminum foil glued on after manufacture to the trailing edge. 
The aluminum foil is inside the blade contour. Before the part is 
manufactured, a strip is Inserted into the mold for the corresponding 
throughput. The leading edge and trailing edge conductors are connected 
at the blade tip by means of an aluminum cap. As an area protection , 

‘a thin aluminum tissue (Alu-mesh) is naated on along the outer third 
of the blade which is then electrically connected with-'the main con- 
ductors by means of research screws . 

When the design parameters were achieved, the blade was further 
detailed after freezing the design and a set 'of drawings was made. 

The set of drawings Include lamination plans which specify the 
laminate structure and the t vsiaiie deposition to be used as production 
room drawings. 
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II. 3 


Material 


High requirements are placed on the blade for large wind energy 
facilities. This involves lifetime, lightness and stiffness. Because 
of the high specific stiffness in strength CPI (carbon fiber reinforced 
plastic) was adapted as the primary structure with T 300 fibers in the 
form of tissue-prepegs. 


As a matrix, we selected the M 10-resin system qualified for avia- 
tion in Prance. This resin system has a pot time of 60 days at room 
temperature and requires a minimum hardening pressure of only 0.5 bar. 

This extremely long pot time allows continuous manufacturing of the wing 
without intermediate hardening. The materials required for building 
the two rotor blades were for the most part already delivered by 
suppliers. Residual deliveries are expected up to April 1981. 

/84 

III. Manufacturing means 


III.l Lamination molds 
III. 1.1 Surfa;Ge definition 

By dividing the surface of the wings into areas, by accepting 
small errors the slightly spherically curved surface of a section can 
be transformed, into a two-dimensional surface using the tangent condi- 
tion (equal inclinations over the same profile points) and can, there^- 
fore, be developed. The developed surfaces were provided with the bend- 
ing lines which depend on the profile contour and this represents the 
development drawings with two cross section drawings each per section, 
as the basis for manufacturing the blade mold shells . 

III. 1.2 Molding shells 

The blade mold consists of two half shells which cannot be closed 
with 5 sections (sections). Each section is manufactured from 8 nun 
thick sheet metal which is bent, which is then attached to a frame 
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made of mero building block elements and which is anchored on the 

: . ‘ ‘ I 

ground. Because of the thermal expansion which occurred during heat- I 
ing of the shell, the support points between the shell ar.d the frame 
are designed so that thermal uncoupling of both components is provided 
for. The sheet ’;metal shell has so little torsional strength that the 
blade turning angle (twist) can be considered when setting up. 

According to discussions in III. 1.1 and 111.1.2, 11 shell molds 
made of 8 mm thick steel sheet were bent as well as 11 mold supports 
made of mero components and these were assembled. The lltb shell mold 
was required for geometric reasons as an additional shell in the con- 
nection region (circular connection cross section of the blade). 

111. 1.3 Mold heating 

The heating of the mold is done using a 1.5 mm thick CB’K-GPK (car- 
bon reinforced plastic-glass fiber glass reinforced plastic) mat. 

After final acceptance of the assembly of the lamination mold, 
a thin GPK laminate was laminated down to the mold surface as an install- 
ation layer on to which carbon fiber tissue paths could be glued and 
to which a unit could be connected. 

L-y applying electrical current, the carbon fiber tissue heats up 
due to its own resistance. 

The heating unit is a transformer-like unit with Installed con- 
trol device which makes it possible to have heating at oonst ant temper- 
ature during the various phases of the hardening cycle and to control it. 

/83 

The described heating technology was developed and qualified by 
extensive tests during which the temperature distribution and its con- 
sistency over the test time was measured as a means for manufacturing 
wings. 

111.1.4 Blade mode installation 

^ 93 ■■■■ 



The blade mold as well as the mold Installation and heating tech- 
nology represent new technologies. Therefore, the first two mold seg- 
ments for the internal blade region (radius points R = 0.8 to R = 11 m) 
were assembled for manufacturing a 10.2 m long blade segment. 

Set-up description 

Assembly of the shell mold-frame 

The shell molds were screwed together with the various bearing 
heads which could be adjusted in height. The holes were sealed. 

Shell mold-shell mode connection 

The shell molds were screwed together vrithout any torque using 
a metal strip bent according to the profile contour so that the shock 
between the shell mold could be sealed with a round silicon sealing 
string. 

Manufacturing of contour templates 

Using the profile cross section drawings already made available 
for mold manufacturing, wood templates with the nominal profile contour 
were produced for the beginning and end of each shell mold. The posi- 
tions of the profile threading axis was characterized by a hole in 
each case. 

Alignment of the molds 

First of all, the wood templates at the radius poiri^ 0.8m, 6m and 
11m were positioned for the best fit. Using the bearii.g element of the 
frames Which could be adjusted in height, the shell molds were positioned 
so that 

a) the holes (threading axis. of the templates lay along one 
line (measurement with a theodolite). 



b) reference lines marked on the templates were horizontal 
(measured with a level). 
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The profile threading axes of both shell molds have to run par- 
allel to one another both to the side and in the elevation direction. 

Adjustment of the molds for flatness 

Since the shell molds are developed surfaces, that is, non-spher- 
ical surfaces, the bending lines have to be straight lines. Theix‘ 
straightness is tested using a ruler and may be corrected by adjusting 
the bearing element. 

Pixlng of the configuration 

After again testing the conf Igjitatioh , the blade molds were fixed and 
the displacement elements were sealed. 

Pisplacement of the rails 

Driving rails for the alignment port were measured at a defined 
side and elevation separation and, therefore, reference to the floor 
of the plant . 

III.l. 5 Quality assurance measures for the blade mold installation 

Gontrol measures were used for acceptance of the shell mold con- 
tour and for the blade mold installation. The control requirements 
specified the extent of controls, the time of controls and the type of 
controls used. 

The corrected configuration state for manufacture of the mold is 
tested continuously by continuous presence of a controller during the 
set-up work and during a final acceptance of the molds according to the 
control specifications. 
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III. 2 Surfacing device 


A surfacing device vfas constructed and designed for milling the 
blade separation plane. The surfacing device consists of a port made 
of mero building components which runs along two rails and covers 
both blade molds. The milling device is suspended in the port which 
essentially consists of a support with two drive units for two pot 
disks in the form of surfacing tools which runs in the blade chord 
direction. The twisting of the blade plane requires variable inclina- 
tion of the milling device with respect to the horizontal which is 
provided for by a guy track to the side of the blade mold. 
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III. 3 Cutting and .drilling devices 

The devices for working the blade ■connection described in the fol- 
lowing were designed, constructed and used during manufacturing of the 
full scale part. Drive for the drilling and cutting tools' is done 
through a central drive unit through a flexible shaft ^ ^ ^ 

111. 3.1 Cutting device 

The cutting device essentially consists of a guide ring which goes 
around the blade root on whlcj) a GUttrln^ device is which can 

move in the radial and circumferential direction. This means that 
manually it is possible to cut through the blade connection laminate 
which are about 90 mm thick. A diamond saw blade was 
used as the cutting tool. 

The positioning of the guide ring over the component is provided 
for by marks in the mold and its impressions are again found in the 
unmolded part. 

111. 3. 2 Drilli.ig template 

The drilling template defines the position of the holes in the blade 
longitudinal direction and the position of the passage holes in the 



radial direction. The drilling template is designed so that the drill- 
ing device can be attached to it. 

III. 3 . 3 Drilling device 

The drilling device consists of a frame which; can be flanged to 
the drilling template with a drilling device and an advance device. 

The advance is done thrpugh a tooth wheel and tooth rod using a manual 
wheel. With the drilling device, it is possible to drill both the 
pocket holes in the longitudinal blade direction as well as the passage 
holes in the radial direction by flipping over the advance unit by 90®. 

IV. Full scale part manufacture 

In order to be able to compare the theoretical properties (espe- 
cially mass and thickness distribution) with the actually realized com- 
ponent properties under manufacturing conditions, a 10. -2 full scale 
component of the blade internal region was manufactured with the orig- 
inal laminate structure for testing purposes. Relevant parameters, 
such as workability of the prepeg, lamination time, heating rates arid 
alignment accuracy of the fibers were examined and specified again. 

IV, 1 Manufacturing of the laminate shells 

IV. 1.1 Laying in work 

Before starting the laying in work, a separation agent was intro* 
duced into the lamination molds containing the heating unit . The GPK 
(glass fiber reinforced plastic ) tissue prepegs were inserted with ariv 
orientation angle of ±4 5° into the mold which were cut for the torsion 
skin and this was rolled on with teflon rollers . Without intermediate 
hardening and starting at the radius point R = 11 m, the CPK (carbon 
fiber reinforced plastic) unidirectional layers were laid in in the 
same way but in a zero degree direction but spliced in the direction 
of the blade connection. The CPK layers provide for the flapping and 
deflection stiffness of the blade. 



Additional CPK layers in the zero degree direction were applied 
directly to the blade connection in order to stiffen it and in the 90° 
direction for reasons of thermal expansion. 

In order to thicken the belt laminate in the lateral blade region 
and, therefore, to reduce the hole internal pressure and for better 
stress distribution, spliced GPK layers spliced into the CPK strip are 
used . 

In order to have an even larger hollow area, an internal GPK wedge 
laminate is applied so that the internal and external belt position is 
separated and this makes it possible to connect a laminate about 90 mm 
thick using a transverse bolt. The sequences of the prepeg types to 
be used, the numbers of layers and the fiber orientation angles are 
specified in lamination plans . 

IV.1.2 Vacuum structure 

After ending the laying in work, the vacuum facility was installed 
which was used for large area pressing of the laminate and for sucking 
off part of the lamination resin. 

The unit itself consists of a foil which covers the entire laminate 
and is sealed off everywhere and the air underneath it is evacuated. 

IV . 1 . 3 Hardening of the laminate 

Hardening of the reaction resin is done ;exothermally . The charac- 
teristic of the temperature time variation wcs determined already before 
with sample laminates of the same thickness as the blade connection 
laminate. 

Thermal elements were installed and laminated into the heating 
laminates which are connected to the control devices of the heating 
unit . This measure assures that at the beginning there will be a uni- 
form heating of the laminate and after the exothermal reaction, there 
will be no additional energy supply by the heating system. The entire 


heating cycle followed the determined characteristic of the temperature /89 
time variation and was continuously monitored* 

XV . 2 Manufacturing of the foam core 

IV. 2.1 Molding and gluing in of the foam core 

After removing the vacuum unit, the foam core was fitted into the 
completely hardened laminate shells. This consists of about 80 mm 
thick PVC foam plates ;-7hich are cut according to the internal shell' 
contour. Prom R = 6 m to R * 11 m, these were glued in staggered in 
the blade core direction. Prom R = 0.8 m to R = 6m, they were Inserted 
in the form of a spar. Gluing of the core to the laminate was done 
using an expandable cold foamable epoxy gluing system (XB 305^ A + 
hardener XB 305^ B + propellant XB 2972 C) which fills cavities caused 
by cutting inaccuracies. In addition, 1 GPK (glass fiber reinforced 
plastic) rib was laminated in for the static test. 

IV.2.2 Surfacing of the foam core 

Milling of the blade separation planes was done using the surfacing 
device mentioned under III. 2. Tool advance was done manually. 

IV. 3 Assembly 

The surface wing path shells were placed above one another, posi- 
tioned along the blade leading edge and then glued with the cold foaming 
adhesive system already described under III. 2 . 1 . The third shell in 
the form of a ’’wedge" required for geometric reasons (see III. 1.2) was 
matched by hand and was connected with cold hardening CPK wet laminates 
with the wing half shells. The torsion complex was also closed using 
wet laminates at the wing nose and trailing edge. 

IV. 4 Cutting of the connection cross section 

The cutting device discussed under III. 3.1 was positioned on the 
wing and the 90 mm thick laminates were separated with a slight excess. 
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The separation surface obtained in this way was surfaced using an 
additional device with an accuracy of ±0.2 nun to the final tolerance. 

IV. 5 Drilling 
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The drilling template and drilling device (see III. 3-2 and 111.3.3) 
were assembled and positioned as a unit at the connection cross-section. 
Using a hard meoal tool, there was a stepwise drilling of the hole for 
the force introduction elements up to a drilling diameter of 45 mm, 

V. Test and prototype testing 

V.l Sample testing and component test 

Before the prototypes were used, all of the new kinds of materials 
(for example, foam system XB 3054)were subjected to techhological tests 
and the mechanical properties were determined by tension, shear and 
torsion samples. 

The basic properties of the blade connection method are already 
well known from the DNW program (German-Holland Wind Tunnel) . We did 
not make any component tests. 

Vv2 Prototype testing 

In order to test the extent to which the theoretical component 
properties agree with the ones found after manufacturing, static and 
dynamic tests were performed with the 10. 2 m full size component . For 
this purpose, it was screwed horizontally to a 3*3 ton tube-like bearing 
block designed for this purpose which itself was mounted on two 8 mm 
high and 13 m long double T supports. The support block is designed so 
that the original connection elements could be used and, therefore, 
the original clamping conditions could be simulated. 

V.2.1 Dynamic test 

The measurement of the eigen frequencies of the component allows 
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one to determine the extent to which the mass and stiffness distribu- 
tion (especially their coupling) can be maintained with a specified 
tolerance width. 

The measurements were performed with a test band oscillation 
multiple position measurement installation. 

Using electrodynamic exciters, the full scale component was excited 
over the range 0 - 200 Hz and all resonance frequencies up to 200 Hz 
and oscillation modes up to 100 Hz were identified with accelerometers 
and measured. 

/91 

Since completely rigid clamping of the blade cannot be realized in 
practice, the motions of the clamping points were also measured by their 
own measurement transducers (accelerometers) for all of the eigen modes. 
The measurements themselves were performed for two positions of the 
blade because the clamping for horizontal blade oscillations bad a 
larger yield. ^ ^ 

We had the following results: 

1st flapping bending 9*84 Hz calculated 8.6 Hz 

1st deflection bending 9 •99 Hz calculated 8.9 Hz 

2nd flapping bending 31*8 Hz 
2nd deflection bending 36 •9 Hz 
1st torsidn 45,5 Hz 

Further results are contained in MBB/FE173/Trl2 . 

The calculated values are taken from a frequency calculation of 
the firm Volth Getriebe KG. 

The test results show that the goals of the contracting agency 
were exceeded and the mass and stiffness specifications are within 
tolerances . 

V.2.2 Static tests 
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The theoretically determined bending stiffness safety factors 
are to be demonstrated in thfe ‘Static test with the full scale compo- 
nent. 

The clamp comportent was loaded by a resulting force at the blade 
radius points K =. 6m and R - 10.8 m so that the flapping and deflection 
moments were simulated at the same time. In addition, a pure deflection 
test was carried out. The blade deformations which occurred were mea- 
sured using strain gauges and path recorders. A single force was intro- 
duced by means of a wooden scissors in such a way that Its line of 
action passed through the thrust eenter. The load was Increased in 
steps up to a maximum of 130 kN at the radius point 10.8 m. 

Prom the strains, the stiffnesses in the region 2.5 m < R <6 m 
were determined. In the other regions, they could not be accurately 
determined because there only small strains occurred outside due to the 
small bending moment and in the root region because of the high bending 
stiffnesses. 

The determined stiffnesses are above the nominal values. Bending 
lines recorded in this region also show smaller bending compared with 
calculations and therefore higher stiffness. 

V. 2 . 3 Fracture test /92 

The fracture test is intended to determine values for the most cri- 
tical failure case considered to be crumpling of the belt laminate 
along the wing side subject to compression. In this failure mode, the 
belt separates locally at a critical pressure stress because the sup- 
port by the foam is no longer sufficient . In thb test the component 
failed at a load of 130 kN ahead of the rib installed at 6m in the 
region of the area vmich was the least safe for the component according 
to theory (only is valid from moment loads by individual forces at the 
component tip). 
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Safety Factors achieved 

a) Local pressure stability 

i 

For the straiits which occurred at a failure load of P = 130 kN, 
we find a safety factor against crumpling for the entire blade of 
J = 1,36 for the maximum positive gust. This value, therefore, applies 
for the "century gust" load case. 

b) Transverse force 

The smallest safety factor against transverse force thrust was 
found for the test component at the point R = 6.5 m: j = 4,0 

c) Safety against bending moments for the WEC entire blade 

Prom the point R = 9- 7 m up to the root region, for the test com- 
ponent a higher bending moment was simulated as occurs during operation. 

d) Blade connection 

With the test load of I 30 kN, the 3.2 time root moment was 
demonstrated. 

Further results are contained in TN DE 133-1/81 


V. 3 Concluding remarks 

With the test described in Chapter V.2, the eigen frequencies, the 
static component strength and stiffness and, therefore, the achievement 
of the goals has been demonstrated... 
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VI. Concept development 

VI. 1 Component development and release 

According to the prototype qualifications discussed in Chapter V, 
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the drawing sets were reworked and were released for manufacturing 
the first -25.2 m wings. Based on knowledge from the full scale manu- 
facturing (Chapter IV) in the meantime detailed manufacturing and con- 
trol specifications have been developed and released.; 

VI. 2 Manufacturing concepts development- 

VI. 2.1 Installation of foam core 

In the manufacturing of the prototype (IV. 2.1), it was found that 
i ‘ it takes a lot of time to individually fit in the foam plates. In 

addition, the filling of the foam plates was not always perfect because 
the pressure required for gluing ( ) could only be applied in the 

direction of the laminate shells. Base on these problems, the foam 
core was glued in as follows in the full scale wings. 

Approximately 3 mm thick foam strips were glued in in the blade 
longitudinal direction into the laminate shells. In this way, the 
foam core has a well defined distance from the laminate shell and only 
has to be roughly tailored. Because of the defined gap, one obtains a 
perfect gluing of the entire foam core surface, and the shell with the 
cold expanding foam system XB 305^* In addition,- about 10 roughly 
tailored foam plates were glued together into a block using screws. 
Gluing of the blocks among themselves is done in a different manner as 
during prototype manufacturing, also using XB 305^»’ Based on these 
Improvements, possible erroneous gluing can almost be excluded. 

VI. 2. 2 Surfacing of foam core 

During the surfacing work on the prototype (point IV. 2. 2), the 
surfacing tool (rasp-like pot disk) was found to be insufficient 
because of reduced chip performance due to chip size and the reduced 
possible advanced rates. 

A tool was designed for milling work for the large wing and it 
was built. The tool consists of a flat aluminum plate inhere chisels 







can be Installed in four holes. With this tool, because of its good 
cutting performance, the milling time of the large wing could be re- 
duced to about 1/4 of the milling time of the prototype. The QPK 794 
(glass fiber reinforced plastic) rib at R = 6 m which is a disturbing 
element for the foam milling (tool wear) is only installed after mill- 
ing. 

VI. 2. 3 Cutting concept development 

The separation of the connection cross section discussed under 
IV. 4 led to steps in the connection plane due to guidance inaccuracies 
and discontinuous separations which had to be removed with a large 
amount of effort. A new device was constructed for cutting the large 
wing and it was manufactured which allows continuous separation of the 
approximately 90 mm safe blade connection laminates and uses a diamond 
blade which rotates in the blade circiamferential direction. In this way, 
the required reworking can be reduced to a minimum, 

VI. 2.4 Drilling concept development 

During the drilling of the prototype (IV. 5), it was necessary to 
predrill up to four different steps which required a great deal of time. 
Based on extensive research with special firms, a hard metal drilling 
tool was found which allows one to make the 30 mm and 40 mm diameter 
holes with suction removal of the dust in one step with the prescribed 
fitting quality. The drilling device was modified in order to use these 
tools . :■ 

VII . Large wing manufacture 

VII. 1 Manufacturing tool completion 

Lamination molds in the >same way as discussed in III. 1.4 were 
installed from R = 11 m to R = 26 m for manufacturing the large wing, 
as were the driving tracks for the surfacing port which was extended to 
about 30 meters . After optical measurements and evaluation of the 
individual values, the lamination molds were released. After lamination 
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of the mold heating units, electrical connection and functional con- 
trol by heating, the cooled molds were again optically measured in 
order to detect possible displacements of the mold. Only after it was 
found that the remaining dlspls^fements were within tolerances were the 
molds finally released for manufacture. 


VII. 2 Manufacturing of the laminate shells 
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Manufacture of the 26 m long laminate shells was done according to 
the manufacturing of the prototype laminate shells discussed under 
IV. 1 without any manufacturing changes. 


VII. 3 Manufacturing of the foam core 

Up to the time of this report (February 28, 1981), the gluing in 
work and surfacing of the foam cores discussed under VI.2.1 and 2,2 
has occurred including manufacturing of the trim chambers. 

VIII. Source list 


Gatter, W. Oscillation measurements in the Internal section 

of the WEC blade MBB/FE173/STR12 
of June 26, 198O 


Pfeifer, K. TN-DE 133-1/81 

of January 23, 1981 

Krautwald, R. Status report on rotcr blade development ET 434 3A 
for wind energy facility "Stoetten 40707b” 
status : March 30, 1980 
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FIGURE 3. TOTAL VIEW OP THE WEC PROTOTYPE MANUFACTURER 
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FIGURE 7. SKETCH OF THE L«MINATION MOLD (FIRST PASS) 
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FIGURE 11. TEARING OFF OF THE PROTECTIVE FOIL 


ORIGINA 


BLACK AND WHITE PHOTOGRAPH 




ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAP 


FIGURE 13. V BELT AT THE COMNECTION REGION 
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FIGURE 16. GLUIRG OF THE FOAM CORE 
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FIGURE 17. ALIGNMENT TOOL DURING MILLING OF THE 
FOAM SPAR 






FIGURES 18 and 19. ALIGNMENT OF FOAM CORE 
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FIGURE 23 . ROUGH STATE OF THE BLADE CONNECTION 
SURFACE FOR RADIUS POINT 0.8 m BEFORE CUTTING 
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FIGURE 2H. ROUGH STATE OF BLADE TIP AT RADIUS 
.11 m (BEFORE INSTALLATION OF THE RIB) 
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PIilURF ?«. PRILLINCI OP THK TRANSVERSE BOLT HOLE DIAMETER 4S 
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FIGURE 30. FLANGE CONNECTION OF COMPONENTS AND 
COMPONENT TESTING ON THE TEST DEVICE 
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FIGURE 34. LOAD IN THE DEFLECTION DIRECTION 
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FIGURE 35. FRACTURE TEST 
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FIGURES 36 and 37. DETAILS OF THE FRACTURE POINT 
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VIEW OF TOWER HEAD 
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ROTOR HEAD IS LIFTED TO THE TOWER 
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EXTENSIVE PROJECT (CONTROL OF PROJECT ON THE MANUFACTURING 
DEVELOPMENT, BUILDING AND TESTING OF THE GROWIAN ROTOR 
BLADES ( ET 4223 A) 

'Project No. ET H 323 B 
R. Lebmhus 

1. Goal 

The goal of the contract is to develop an ordered sequence for 
the work performed under the project ’’Production development, building 
and testing of the Growlan rotor blade”. 

2. Work program 

Examination of calculations, specification drawings of the rotor 
blades : 

1 . Evaluation and control of testing of design calculation 
drawings . 

2. Production supervision of the rotor blades and testing 
during manufacturing . 

3. Monitoring of the testing program for investigating a test 
blade . 

3 . Status 

Item 1) 

The building specifications, calculation data, construction data 
and production data for the steel spar of the rotor blade prototype 
presented by MAN were examined. An operational strength calculation 
(MAN) was not yet available. However, discussions were carried out 
about the way in which the operating strength demonstration would be 
implemented. 

*^ProJ ect no. ET 4323 B 
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The data submitted up to the present by the MAN firm were pri- 
marily design data. MAN states that it can only present the final 
data during the first half year 1981 which will be the subject for 
evaluations. 

The previously submitted strength demonstrations essentially 
agree with the results of our own FK investigations. Our calculations, 
however, show that in the range R = 35 m (transition from transverse 
stiffening unit of the spar to longitudinal stiffening member) one 
would expect locally limited and impermissible stress excesses. At 
the present time it is assumed that these can be reduced to an unaccept- 
able level by design changes. 

Item 2 ) /133 

The production of the prototype steel spar at the Gustavsburg 
plant of MAN was examined by supervisors of the German Lloyd firm as 
regards material used, individual part manufacturing, assembling includ- 
ing Welding work and evaluation of the results of the welding seam 
inspections, using random samples. 

With one exception, the case of the welding weld at Hetebrueg, no 
important discrepancies in the steel spar were found as concerns the 
execution of the design. /I 3^ 

Item 3) 

We expect the testing of the prototype blade at lABG Munich not 
before the second quarter 1981 according to the present plan of MAN. 

At the present time, we have reviewed the preliminary planning data for 
carrying out tests. 

- Problems which occur 

The development of data which can be used for examination by the 
MAN firm has involved substantial time delays in the advance of the 
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Blade bearing tip fcr test blade at the MAN-Werk Nuernberg. Mr. 
Hetebrueg of Germanischen Lloyd is measuring the upper diameter of 
the several meters long bearing and its deviation from circular 
shape (2050 ±0.5) mm. December 198O 

(Photograph: R. Windheim) 
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production of the prototype blade. The MAN firm was told that possi- 
ble consequences of this may have to be borne by MAN. 

- Deadline situation 

According to our license, the project running time extends from 
November 13* 1979 to April 30, 1981. Since this project represents 
a task which accompanies FE work at the MAN firm (ET 4323 A), we can 
only perform our part of the project on time if the MAN firm holds to 
the schedule. At the end of 1980, it is clear that the MAN firm cannot 
maintain the planned schedule for their task. /135 

DEVELOPMENT OP A MEASUREMENT TEST 
PROGRAM FOR GROWIAN* 

Friedrich Koerber 

CONTENT /I 36 

“ l.v Summary 

2. Classification of measurement program 

3. Previous measurement -programs 

4. Purpose of GROWIAN measurement program 

5. Tasks of the measurement program 

6. Variables of interest 

7. Measurement points 

8. Measurement philosophy 

9. Data recording and processing systems 

10. Data evaluation 

7 figures 
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Iv SUMMARY 


Under the title "Development of a measurement and test program 
for GROWIAN", at the end of 1979 a contract was awarded to the MAN- 
NEUE TECHNOLOGIE firm to develop a program for test operations for the 
GROWIAN wind energy facility in the -planning stage. 

The present report discusses the measurement program and shows 
solutions. 

The main topics of the measurement program include the identifica- 
tion of the facility behavior, verification of the facility concept 
including models and the collection of operational experience with 
large wind energy facilities. . Primarily, we are Interested in the 
action of the environment and the reaction of the facliity based on 
properties. The data selection and processing required for the Work 
is suggested. 

The purpose of the measurement program is to obtain knowledge with 
which technical problems of large wind force facilities and their econ- 
omy can be better evaluated. 

The ET ^ 36 ^ A project Includes a complete measurement program. 

The extent had to be restricted in order to save costs. The hasic 
ideas of the program however remain unaffected. /I38 

For simplicity, we will call the Growian measurement and test 
program simply the measurement program in the following. 

2. CLASSIFICATION OF MEASUREMENT PROGRAM 

Figure 1 shows the situation. 

For the first phase of the overall project of the large Growian 
wind energy facility, there was included the development of "data ready 
for production". 
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At the present time the second phase is being implemented with 
the realization of the facility. In parallel with this, there is the 
development of a measurement and test program for Growian. In the 
next step, the measurement installation has to be integrated during 
construction. The measurement program specifies which measurement 
types are to be carried out and where certain measurements are to be 
performed at various points of the facility. We discuss the recording, 
processing and evaluation of the data. 

In the third section we will discuss a multiyear measurement and 
test period, the purpose of which is to collect •operational experience 
with a large wind energy facility. We will examine assumptions and 
will expand the state of the art in thfe calculation and design of large 
facilities . • 

3. PREVIOUS MEASUREMENT PROGRAMS AND WIND FORCE FACILITIES 


The review of previous work must precede »any new task. 

Pacllitles.'-wlth the size of the a.re only in their 

beginning stages. From the power level point of view, the American 
facility MOD-1 with 2 MW is the closest to the Growian. The first 
machine started in the summer of 1979 in Boone NO. At the present time 
the first measurement results are being published. For the most part 
they confirm the calculations of the design of the wind mill located 
there . Extensive data are available about the MOD-0 facilities (100/ /I 3 9 
200 kW) in the USA which have been tested since 1975. However, their 
power level is less than 105? of that for Growian. Nevertheless, mea- 
surement procedures and results are a point of reference for the planned 
German pro j ect . 

Extensive measurements exist for the 20 year old Gedser wind mill 
(200 kW) . A measurement program started with two other new Danish 
wind mills (650 kW) and results are already available. Also measure- 
ment activity has started at the . TVIND wind mill . 
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Finally, we should mention experience which was selected during 
the 60 ’s on the 100 kW wind generation facility of the Wind Force 
Study Association as well as work in Sweden and Holland. 

Summarizing, we have to point out that detailed investigation 
results are primarily available from smaller facilities and can only 
be transferred to a facility of the size of thfe ^Growian with reserva- 
tion. This is especially true for determining the wind conditions, the 
component -loads and the operational behavior of the wind generation 
facility. The statfe of the art cannot give any information about the 
long time quality of wind energy or the composite operation of wind 
generation facilities for producing electrical power. 

i|. PURPOSE OF THE GROWIAN MEASUREMENT PROGRAM 

It is the purpose of the measurement program to obtain Information 
which will allow the understanding of the technical problems which 
then should lead to better solutions. This is both true for the selected 
Growian design as well as for basic boundary conditions of wind energy 
exploitation. For example. Wind properties or facility types. Together 
with the operational experience with the prototype facility during 
test operation, we wish to obtain information about the economy of the 
wind flow from individual facilities and wind generation parks. /l^O 

A number of steps will have to be followed in order to obtain these 
goals. Figure 2 shows some of them. In the design and construction of 
the Growian, assumptions and models have to be used. The first goal is 
therefore to verify the models. This means that th6 operational facil- 
ity has to be identified so that its realistic properties become known. 
This'then has to be compared with the assumptions and corrections have 
to be introduced. One would expect that, due to lack of information, 
several influences cannot or were not sufficiently modeled* As an 
example, we should mention the only approximately known damping property 
of the facility structure or its effect on loads . Also there is the 
influence of gust increase and decay They are important influences 

on the dynamics and would therefore have to be included in a load analysis 
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IDENTIPICATION OF GROWIAN PROPERTIES 
VERIFICATION OF MODELS 

CLARIFICATION OF ACTION OF THE ENVIRONMENT AND 
REACTION OF THE FACILITY 

TESTING OF THE DESIGN 

TESTING OF CONSTRUCTION DESIGNS 

COLLECTION OF OPERATIONAL DATA 

DEMONSTRATION OF LARGE SCALE USABILITY OP WIND ENERGY 
INPORMATION ON ECONOMY 
IMPROVEMENT OP FUTURE MACHINES 
OBTAINING ADDITIONAL DATA 


FIGURE 2. PURPOSE OP THE MEASUREMENT PROGRAM 



for the installation components. 


Therefore, the model has to be refined so that effects not accur- 
ately represented in the model can be included. For example, this 
includes rotor instabilities or the realistic effect of the tower wind 
shadow. This influences the power yield and, therefore, occurs as an 
exciting element for oscillations of the overall systems. 

The measurement -ppogram is tailored to offer generalization of 
models, a trend which is developing. Comparison of model calculations 
and measurement results Is, therefore, supposed to lead to a non-refined 
model which then will be used to give rough estimates of future machines. 
This model will have to be free of unnecessary ballast and will have a 
low amount of input information. 

5. PURPOSES OF THE MEASUREMENT PROGRAM 


The purpose of the measurement program defines the problems it will 
solve. During the Growian test operation, a number of test series will 
be performed. The measurement program plans and coordinates multiyear 
Investigations by defining tasks and procedures. 

/I4l 

Growian was conceived as a standard unit of a wind energy facility. 

Its function as a measurement carrier during test operation means 
it is necessary to tailor the design to the requirements of measurement 
technology. This results in design changes which have to be traded off 
with the overall design. 

The measurement plans for the Individual parameters are the most 
extensive part of the measurement program: 

- wind conditions 

- facility power output 

- facility kinematics and dynamics 

- loads 

- behavior of electrical system 

- operational control 

- emission and immission^^^^^^^^^ 
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A large niimber of data are obtained during the measurement pro- 
grams, Therefore, it is also necessary to develop a concept of data 
recording, processing, storage and evaluation. Finally, a test and 
measurement catalog will be developed with which the test can be per- 
formed . 

6. PARAMETERS OF INTEREST 

The engineers are interested in the following three quantities of 
a wind energy facility: 

- power - loads - emission 

Once these are known, it is possible to develop the Information 
required from the measurement program. 

/1 4 2 

Due to the action of the environment and system properties, one 
then can develop thte parameters of the Growian reaction. Figure 3 shows 
this relationship and the arrows indicate the directions. 

On the left column wt show the contact points bf Growian and •the 
environment. For the most part we understand by this the action of the 
wind and the action of the network to which Growian supplies power. 

The wind is most easily described by means of its velocity vector. 

Its absolute magnitude and direction change at a fixed point over time. 

In addition, the velocity vector at a fixed time changes greatly as a 
function of position. ^ field changes, therefore, in three 

space dimensions and in time dimension. Work previously carried 
out under the Growian project again always showed that the wind field 
and its time change is a large and uncomfortable unknown which becomes 
very important considering the size of the facility. The measurement 
program, therefore, stresses this point. 

The connection of the wind driven generator to the available medium 
voltage, electrical network is another important point of contact of the 
fLcllity. During normal operation, there are strict requirements for 
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FIGURE 3. ACTION AND REACTION 




the quality of the energy delivered: constant frequency, uniforiu power, 

low higher harmonic component. During disturbances in the net, a state 
is required from the facility which will allow resynchronlsatlon in an 
orderly manner. On the other hand, in the case of frequent net distur- 
bances and, therefore, sudden removal of the load, which especially 
occurs during storms, must not be endangered. The network is very 
important because the usability of the large wind energy converters 
during ;jolnt operation depends greatly on the way in which it is 
connected to the network, especially during critical network states. 

Influences of the environment have to be experienced and this will 
especially include natural influences. 

Let us again refer to Figure 3 and consider the left side of the 
page. This shows a list of system properties of the wind generation 
facility . The absolute dimensions of Growlan, the facility concept and 
the component types, the type of operational modes and the overall pi'o- 
pertles of the system influence the three quantities: power, loads, 

emission. 

7. M EASUREMENT POINTS 

According to the previous chapter, in order to determine the para- 
meters of interest measurement points are selected over the entire wind 
generation facility. This includes transducers for mechanical, thermal 
and electrical measured values. Figure *4 gives a summary about the 
local eorrespondence of the measurement points and the obsei’vation 
points. 


The measurement 

points Include the following processes or states: 

strain 

force 

acceleration 

path 

speed 

temperature 

voltage 

current 

frequency 

pressure 

torque 

rpni 

flow speed 

flow direction 
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global radiation 
acoustic pressure (level) 
surface state 
tremors 

operational measured values 
perception 

/IHH 

The measurement intervals during the measurement program vary 
widely. All of the data whicl:| are involved with the global energy 
supply and energy gain, lifetime and economy, have to be collected 
over long periods of time. The time scale can extend from one day 
(only for energy analysis) up to one or several years. These measure- 
ments are called long time measurements. 

The study of the loads caused by the action of the envirGnment 
will involve shorter time periods . The cyclical return of load , for 
example, due to rotor motion or gusts specifies a measurement duration 
of seconds and minutes. The analysis of the ccnverslon of the flow 
energy of the air into the rotory energy of the rotor and the electrical 
energy delivered by the generator or transformers will involve this 
time scale. The influence of daily rhythms, such as air temperature 
and sound propagation, leads to a measurement duration on the order of 
one day. The prescribed order of measurements can be performed arbi- 
trarily if one waits for the meteorological boundary conditions and if 
the corresponding state of the machines is implemented. These measure- 
ment sequences are, therefore, called measurement campaigns. 

A wind generation facility is subjected to natural forces . Extreme 
cases of wind conditions have important influences. They occur sto- 
chastically and do not fit in with the previously discussed measurement 
processes. Also special operating states of the facility itself can 
occur, even in conj unction with meteorological special cases . These 
special operating states are controlled by the control system within 
the safety monitoring activity. The measurement program is intended to 
give information about the sequence and optimum possibilities of 


flow amount * 
humidity 
switching state 
PS signal 
network state 
optical 

8 . MEASUREMENT PHILOSOPHY 
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controller involvement. In addition to long time measurements and 
measurement campaigns, we therefore introduce a third kind of measure- 
ment process. This is called event measurement and is continuously 
ready because of its task. This process begins when a limiting value 
of a number of exceptional measurement signals J.s exceeded. The record 
Ing duration of relevant measurement values then extends from several 
seconds up to 60 seconds. By the continuous buffering Of all measured 
values, we then also have available the previous history of the event. 
This method of event-triggered initiation of measurement recording will 
make it possible to record the causes, the event itself and thO ;conse- 
quences. 

Figure 5 shows the correspondence between measurement topics and 
measurement modes. A solid dot Indicates the emphasis of the activity, 
the smaller point shows processing of the topic . Figure 6 gives a 
summary of the measurement modes plan. We can read off the following 
characteristics : 

long time measurement: measurement during real operation. 

measurement campaign: measurement of selected states. 

event measurement : measurement of extreme states with 

automatic triggering. 

9 . DATA RECORDING AND PROCESSING SYSTEM 

The design of the wind energy facility GROWIAN and the tasks of 
the measurement program with the measurement modes require a relatively 
close framework for the data acquisition system. Figure 7 gives a 
summary. 

The data are distributed in space over the wind measurement tower 
1, the wind measurement mast 2+3 (measurement grid), the facility 
itself and the operating building. In addition, there is a mobile 
outer station which is especially used for collecting emissions. 

All of the measured values can be divided into groups. These 
groups correspond to the local positions of the measurement points . 

Most of the data Involve the rotor blades, the hub and the rotating 


These 


units of the blade displacement device and the drive train, 
data are collected in a data concentrator 1 and are merged into a data 
single stream. A data concentrator 2 is located in the machine house. /1^6 
There it receives the incoming measured values and combines them to an 
additional data single stream. The same is done by data concentrator 
DK 3 in the tower, DK V in the Operating building, DK 5 in the mobile 
outer station. In the data station the data individual streams are 
collected into a collective stream. 

The number of measurement points (about 500) and their sampling 
rate (l-50Hz, in exceptional cases 160 Hz) suggests the use of a PCM 
facility. This makes it also possible to classify the individual data 
stream into a single total data stream. The subsequent distribution 
then leads the relevant measured values to the evaluation and recording 
units for long time measurements, measurement campaigns, event measure- 
ments and other tasks. 

By dividing up the dedicated computer, we wish to achieve a greater 
operational safety and better handling capability. For a measurement ■ 
and test project of this extent, it makes sense to select a flexible 
concept for data recording and processing because the operational 
experiences certainly will require modifications. This can be better 
satisfied by a building block scheme than with a complex large facility. 

10. DATA EVALUATION 

In order to describe the data evaluation, we again have to dis- 
tinguish by measurement modes. 

The long time measured values are continuously evaluated by the 
program depending on the position. The status of the statistics must 
be callable in a short time. 

The measurement campaigns are evaluated beforehanci in the facility 
for orientation. The detailed analysis of the results with all scienti- 
fic and technical aids will be done by specialists. The results of 
event measurements are treated in a similar way unless a valuation at 
the location is required. 
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MEASUREMENT TOPICS AND MEASUREMENT MODES 














LONG TIME 
MEASUREMENTS 


! measuremen'C under true j 
operational conditions. I 
Time duration: day, month, Year(s; 


action: 

environment (wind, network) 

reaction : 
facility behavior 



wind power gain | 

network compatibility 
operational control 
load collective 
acceptance 
fatigue/ corrosion 


MEASUREMENT 

CAMPAIGN 


measurement at selected states 
time duration: second, minute, da',’ 


[action: 

environment (wind, network 


radio traffic) 


reaction: 

component (facility) 
behavior 



wind structure 

energy conversion 

dynamics 

loads 

control 

emission/acoustics/ 
radio disturbance - 


EVENT MEASUREMENT 

measurement of extreme states 
with automatic triggering 
time duration: seconds. 

action: 



environment (wind, 
extreme) 

network. 

wind structure 
\ dynamics 

[reaction: 


load 

[component (facility) bebavioi 

r control 


FIGURE 6. EXPLANATION OP MEASUREMENT MODES 
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data concentrator 1 



DATA ACQUISITION SYSTEM 




We would like to show this cartoon (publisher’s) about the 
measurement and test program for GROWIAN which th« 
also showed. 
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BLADE TECHNOLOGY PROGRAM 
FOR LARGE WIND ENERGY FACILITIES* 


/155 


Dieter Muser 

Project -goal /156 

Within the blade technology program for large wind energy facil- 
ity tasks, optimum price rotor blades are to be developed. This means 
cost minimization of rotor blade structure as well as an optimization 
of the rotor blade aerodynamics. Interactions between the structure 
and the aerodynamics are to be considered in more detail than before. 

Schedule situation 

The project started at the beginning of the task in October 1979 
and will run to the end of October 1982.> According to schedule, at the 
end of 1980 , the optimizations of the blade structure and the concept 
catalog will be developed. Part 'Of the material investigations were 
done ahead in order to obtain better inputs for the structural calcula- • 
tions. 

Program structure 

The program was divided into four main tasks 

- optimization 

- concept catalog and concept collection 

- concept design 

- test program 

During the optimization, statistical data are collected in order 
to obtain evaluation data. The aerodynamic optimization is still going 
on because no suitable computer models were available for displacing 
the blade sections and for the use of aerodynamic aids. In the /157 
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r WING TECHNOLOGY PROGRAM I 

AND STRUCTURE 



OPTIMIZATION 

- 

POWER STRUCTURE 


- 

STRUCTURE 


- 

COST 

CONCEPT CATALOG 

— 

STEEL 


- 

LIGHT METAL 


- 

COMPOSITE FIBER 


- 

OTHERS 

CONCEPT SELECTION 



TEST PROGRAM 

- 

MATERIALS 


> 

MANUFACTURING 


- 

TESTING 


OPTIMIZATION OF 

ROTOR BLADES 

PERFORMANCE 

- AERODYNAMICS 

- GEOMETRY 

- ADJUSTMENT 

STRUCTURE 

- MATERIALS 

- CONSTRUCTION 
METHOD 

COSTS 

- BLADE 


- FRACTION BLADE 

- POWER-^STRUCTURE 


structural optimization, steel, light metals and composite materials 
were considered which were used with various structural concepts. 
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CONCEPT CATALOG INITIAL DATA 

GROWIAN BLADE DESIGN - GEOMETRY 

- PROFILING 

- LOAD CASES 

GROWIAN - MAX. BENDING 3 m 

(LP 2.$ negative gust) 

- EIGEN FREQUENCY > 0.8 Hz 

TECHNOLOGY PROGRAM - BLADE BEGINNING R = 10.85 m 

- SEPARATION POINT R = 2 7. 50 m 

- SUITABILITY FOR THE SERIES 


Cost •determination based on single element costs was incorporated 
in the calculation program in order to obtain the influences of various 
structures . 

/158 

As initial values for the solution catalog and for predimensioning, 
values of Growian were taken over and our own data was added. 

With these data the following were calculated for the various basic 
structural concepts 

- mass 

- bending 

- eigen frequencies 
; - cost 

During the preliminary construction work, various details were pro- 
duced at a scale of 1:1 in order to find the economical possibilities, 
for example, the spar — end box connection or the end beam. These 
points differ because of the strong profile changes over the blade 
length and are only conditionally suitable for rational design. 
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***silicx>n (?) reinforced plastic 


The calculation of the blade masses results In clear advantages 
for the fiber composite rotor blades. 

/160 

This mass difference also seems to trickle down to the cost. Even 
though the composite fiber materials are much more expensive both for 
procurement and machining, drastic cost reductions can be expected by 
using machines in mass production. More exact data can only be obtained 
after conclusion of the first component test. 

From all of the concept variations, according to the evaluation, 
two similar designs evolve for the further work. 

The D spar up to of the profile chord will use winding techno- 
logy and is attached with the end box using hose winding building tech- 
niques and can be separated from the spar. 

In the case of the box spar, between 5 and ^05t of the profile chord 
will use prepeg technology for manufacturing with nonremovable and glued 


ORIGINAL PAGE JS 
Of POOR QUALITY 


163 



























on sandwich end box. 

Both concepts have a sandwich structure in order to stiffen it 
against bulges in the spar region. By using protruding profiles, it 
becomes possible to reduce manufacturing times and premanufacturing 
of the elements with machines will bring about a uniform high quality. 

/I6l 

By using these high value materials and by exact dimensioning, 
aojltlonal mass savings could be achieved. For the D spar we found a 
mass reduction of ^800 kg for the preliminary design without connection 
elements at the separation points. 


16k 



ORIGINAL PAGE 13 
OF POOR QUALITY 



* GFK = fiber reinforced plastic; ** GST< - glass fiber reinforced plastic 
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Additional work includes detailed structural calculations for 
the separation points and the laminate structure as well as aerody- 
namic calculations of blade design. The manufacturing of the compo- 
nents in the form of partial shells will start in the second quarter 
of 1981. Apparently, the entire project will be carried out within 
the framework and scheduling constraints. 
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Status of pro.lect ET ^342 A. construction and operation of GROWIAN T 


/164 

Lecture during the seminar on building, testing and development of 
large wind energy facilities within the framework of projects supported 
by BMFT. 

MAN Neue Technologie, Munich, March 23-24, 198I. Author: H. Witt. 

Contents 1. Introduction 

2 . Status of project tasks 

3. status of proj ect 

4. Change methods 

5. Location 

6. Official permits 

7. Scheduling situation 

1. Introduction 

Based on the ’’production data”, which were developed during the 
first development phase of the Growian I proJ ect , at the beginning of 
1979 its realization was started. This second development phase was 
structured in several parallel research and development projects in 
order to reduce risks and because of the schedule. Among these we have 
the ’’construction and operation of Growian I” project. 

The "Large Wind Energy Construction and Operating Association 
m b H" has taken on this task. The following are participating in 
this : The Hamburg Electrical Company AG, the SCHLESWAG AG and the 

Rhine Westphalia Electrical Company AG. The company was founded in 
January 198O . 

2. Status of the proj ect 

In January 1980, the firm MAN in Nuermberg was given the contract 
about engineering work on the project based on the production data. 

In August of 1980 it was given the contract for delivery and assembly 
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of the system components. This includes: The rotor hub, the machine 

hubs, the control equipment and the electrical technology. 

/165 

The firm Maurer Sons, Munich, was given the contract for the 
steel tower delivery. 

The operating buildings will be built using the space element con- 
struction method and are being delivered by the firm Klrschner, Duelmen. 
These elements are being manufactured. The HS switching rooms, the 
transformer boxes and the emergency power diesel unit are now being 
awarded to contractors. These use conventional construction methods. 

At the beginning of April 198 I we expect to let contracts for the 
foundation work for the tower and the guide wire foundations. 

At the present time bids are being received for conventional 
electrical technology and the two 170 meter measurement masts, the 
building of Which belongs to this project. 

3 . Status of proj ect 

3.1 Rotor hub 

- The pendulum frame and traverse construction group has been modi- 
fied to correspond to the chain rotor blade bearing concept and 
blade adjustment device. Manufacturing heu stai’ted. 

- The design data of the rotor bearing are available to the AG 

for examination. The manufacturing can be started shortly. 

- The design detailed processing of the blade displacement device 

with hydraulic emergency shut down and the blade retraction ^ 

3 . 2 Machine house 

calculation and construction of the supporting elements 
has been concluded >except for the final design of the connection 
to the tower head. 
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- The following are in process: 

Concept and configuration of the oscillation damper/damplng_ 
system in the front part of the machine house, 

concept for auxiliary construction for the assembly equalization 
ballast, 

concept of the wind direction following device. 

Installation of the hydraulic lifting device 

3 . 3 Control technology and electrical technical equipment 

The concept and clarification of the automatic systems (opera- 
tional control) of the measurement installation and the control for the 
emergency shut down system as well as the facility behavior during power 
shut down has been started according to the design stage. /I 6 6 

3.4 Tower construction 

The detailed construction of the steel tower component elements 
has been started with consideration of thfe devices required for the 
lifting process and the support . 

3.5 Operational buildings 

The planning has been coneluded. The building work at the con- 
struction site have started with the digging of the foundation. 

4 . Change methods 

In order to be able to cover design changes, for example, based 
on technical developments or from calculations and to be able to deter- 
mine their influences on the design of other components and thfeir 
effects on schedules and costs, a special change procedure has been 
developed. 

5 . Location 

The authorization for building the facility at the Kaiser Wilhelm 

170 ■ , 



Koog was obtained formally In 1979 and is now concluded. Considering 
the large area for the support cables, approximately a 6 ha area had 
to be leased. 

With consideration of the present ground conditions, it was found 
that a low foundation of the buildings was the most economical solution. 

The tower and its cable support have local concrete driven in posts 
with driven in base plates which are arranged vertically and obliquely. 

The ramming depth is about l8 m. The measurement mass foundations 
will have profile seal posts. The access road was started in the fall 
of 1980. The asphalt cover will be applied in March of this year. 

6. Official authorizations 

The permit to build the access road was given on October 21, 198O. /I67 

The permit for the entire installatloh was submitted on December 
22, 1980 at the authorities. First partial authorizations are expected 
Shortly. In agreement with the district building authority at Heide, 
the lowest authority, the German Lloyd firm will be used as an official 
authority for testing the safety aspects of components. | 

I 

7 . Schedule situation I 

. ■« 

'■I 

s, 

Today it looks as though no difficulties will occur with maintain- | 

ing the schedule which plans for delivery on November 30, 1982. Delays I 

in the projects for parts of the azimuth bearing can be equalized by I 

measures for accelerating the production schedule. 
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Full sccile GRCWIAN wind facility at the Kaiser Wilhelm-Koog, Elbe Delta. View in 
tte train wind direction towards the Elbe. The 3 wind itBasurement nasts are to be 
bailt during the measurement and test program (see page 135) . Performance of the 
facility 3000 kW after 12 n\/sec. Wind speed, rotor diameter; 100.4 m, tcwer height 
100 m. Planned completion, fall of 1982. 
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INVESTIGATION OP NEW CONCEPTS 
FOR LARGE WIND ENERGY FACILITIES IN THE CLASS 
OP UP TO 135 m ROTOR DIAMETER* 

(ACCORDING TO 5 ^Wg^^g^^p . AT A WIND SPEED OF 11 m/sec.) 

Building and testing of a demonstration facility on a 1:3 scale 
Development of design data for a prototype full scale facility 

R. Meggle 


1. Purpose 

The Investigation of new concepts for large wind energy facilities 
In the class of 135 m rotbrdlameter (according to 5 MW delivered elec- 
trical power at a wind speed of about 11 m/sec). Testing and building 
of a demonstration facility on a model scale (about 1/3). Development 
of design data for a prototype full scale facility. 

2. Working program 

a) Preliminary concept definition for a two- and a one-blade system. 
System design. Investigations of the oscillation and stability 
behavior, Determination of loads. Selection of preferred 
concepts . 

b) System optimization of full scale facility - definition of the 
demonstration facility and production of design data. 

c) Hardware definition of full scale facility - component manu- 
facturing and partial ilntegration of the demonstration facility. 

d) Hardware development (building and testing of selected test 

^ ^ - assembly and operation 

of demonstration facility. Test operation. 

e) Preparation of design data for the full scale facility - dis- 
assembly of the demonstration facility and recultivation. 

*'proJect: No, ET 
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3. 


Status (end 1980) 


1 ) 


- Work on the full scale facility/ (5 
Design studies for rotor blade and head 

- Work on the demonstration facility (370 

- Generation of design data 

- Advertisement of the main component groups. Contracts to the 
most' important contractors 

- Beginning of manufacturing for rotor blade, rotor head, gondola, 
generator system and measurement and data systems, inclusion of 
building interrogation for the Bremerhaven location. 

/17Q 

GROWIAN II Program summary 


-Goal 


Inysstlgations of new concepts for large wind energy facilities 
for the power class >5 MW^^ for the rotor class >135 M0. Building 
and operation of a demonstration facility for a model scale of 1:3. 


- Sequence 


program beginning 7/78 

concept selection 3/79 

design data demonstration facility 12/80 

production/assembly of demonstration 
facility 12/81 

sample operation demonstration afterl2/8l 
facility 

design data for full scale facility OPEN 


- Design goals (for full scale facility) 

specific investment costs according to 
power production 

Design for wind conditions in the North German coast area. 
>20 years lifetime 
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completely automatic operation 
small operating and maintenance costs 

- Location of demonstration facility ; Bremerhaven-Weddewarden 

171 

2) SELECTED CONCEPTS * 

- Single blade system^ free wheeling with horizontal axis 

- Static imbalance equalization by counterweight^^ 

- Rotor blade with composite fiber design: Shell design with 

support foam 

- Pendulum hub 

— Blade angle adjustment (redundant) of the entire blade Including 
counterweight 

- Azimuth following by redundant actuator units 

- RPM yielding force transmission/generator concept 

- Stiff shaft 

- Transmission gear 

- Slipping cup link, rotor brake. Jointed shaft 

- Synchronous generator 

- Direct current Intermediate circuit, static frequency rectifier 

- Computer control operating control system with separate 
emergency off logic (fail safeO 

- Soft (over critical) tower concepts: cylindrical steel tower 
with 3 guy ropes 

- External lightning protection and protection for the electrical 

and electronic systems ■ 


Detailed descriptions can be found in the status report on wind 
energy, VDI Publishing House, Duesseldorf 1980, publishers remarks . 
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Key: 1- — 1st elastic deflection eigen frequency (rotor); 2-^lst 

elastic flapping eigen frequency (rotor); 3— 2nd tower bending; 

^ — alternate; 5— "2nd tower bending; 6 — 1st tower torsion; 7“~support 
rope; 8~lst tower bending; 9 — rpm range, 

10 — rptor rotating frequency 
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KEY TO FIGURE PREVIOUS PAGE: 

1 — •Independent power source; 2— switching box Space cell; 3 — tower; 

4 — elevator; 5 — switching box gondola; 7 — space cells; 8-.-wind measurement 
tower; 9— crane; 10 — grease lubrication pump; 11 — hydraulic brake; 

12 — ventilation/heating; 13 — battery; 14 — facility; 15 — supply; 

16^ — supply; 17— wind measurement installation; l8 — azimuth drive; 

19— blade angle drive; 20 — motor auxiliary drive; 21 — oil pump gear; 
22~disturbance generator system; 23 — power supply operation; 24 — 
braking power operation; 25— network supply station; 26 — synchronous 
frequency converter; 27--slip. rings; 28 — operating control system; 

29— measurement and data system; 30 — rectifier; 31 — alternating 
rectifier; 32— static net switching device; 33-'-rectifier generator 
control 
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l8 — space cell; 19 — wind measurement facility; 20—10 analog inputs; 
21— test and monitoring device; 22 — •disk storage 
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OPTIMIZATION OP LARGE WIND ENERGY FACILITIES j 

John H. Argyrls and Kurt A* Braun ’ : 

t 

Suinmary I 

r 

The project -was carried out jointly by the University of Stuttgart I 

I 

and the Statics and Dynamics Research Association (m.b, H.) with the I 

supervision of the Institute for Statics and Dynamics of the Aviation I 

and Space Plight Design Department. The currently running phase I of I 

the project started on November 11, 1979 and was extended by six months | 

up to December 31 » 1981 within the financial constraints. | 

i 

I 

During phase I, by developing new computer methods and by build- j 

Ing the rotor hub including the data selection systems, conditions | 

will be generated for optimization of large horizontal axis wind energy j 

facilities. I 

During phase II which will not be discussed in the present report, | 

i ; 

various rotor concepts will be optimized and compared using the calcul- ' |- 
at ion and test tools developed. Figure 12 gives a sketch of the open | 

air facility from the present day point of view. /190 | 

. I ' 

[ 

Working points j 

The OPTIWA project can be structured along three lines: 

1. Development of computer programs for the static and dynamics 
investigation of various rotor .concepts. 

2. Calculation methods for a new tower concept and design of a 1 

I 

model tower. r 

■ I 

3. Design and building of a computer control functional model | 

with special rotor hub for simulation of various design and | 

control concepts. I 


1. Development of computer programs 


Here computer programs are developed by the computer engineer 
which allow static and dynamic Investigation of various rotor concepts 
(rigid rotor, pendulum hub. Individual flapping and/or deflecting with 
and without blade feedback, etc.). Precursors for these programs were 
already developed within the research project on ’'investigation of 
rotor stressing and smoothness of operation of large scale and wind 
energy conversion system" which was supported by BMFT and the lEA , 

(International Energy Agency) (see [1-3]). The programs were tailored 
to the requirements of the design engineer In order to allow a fast and 
simple change of characteristic blade data on the monitor (see Figure 1). 
These are based on finite element packages ASKA [4,5]. Since we will 
discuss the calculation methods In [6,7], we will not discuss them here. 

2. Calculation methods for a new tower concept 

Under 2, calculation methods are developed and a model tower Is 
calculated for a new kind of tower concept (Stuttgart tower) which Is 
based on an Idea of Professor Wortmann [8] . Based on the preliminary 
Investigations In [9], the kinematics of the tower shown In Figure 2 
are Investigated and a static design is carried out using the nonlinear 
finite element package LARSTRAN [10]. The dynamic Investigations, /191 

resonance and responses are also carried out using LARSTRAN (see Figure 
3). 

The following disadvantages and advantages can be discerned com- 
pared with conventional tower concepts: 

Advantages : 

- the ground resonance problem is made less important with the 
rotors under consideration because the thrust' — dependent eigen 
frequency of the lateral oscillation (azimuth oscillation here) 
is far below the rotating rate and the eigen frequency of the 
vertical oscillation (pitch oscillation here) can be maintained 
far below the rotating rate. 
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- The tower orientation relative to the wind is always the same 
which allows a static, dynamic and aerodynamic matching of 
the tower cross sectional geometry. 

- The tower will align itself in the direction of the wind with- 
out azimuth drive but may need a yaw damping. 

- The rotor free wheeling condition is good and the rotor axis 
can probably be placed horizontally (tower inclination about 
8.3°). 

- The tower-rotor blade distance is larger which is favorable 
for the tower shadow problem. 

Disadvantages: 

- The construction complexity is probably higher than for a con- 
ventional tower. 

The mast length of the investigated model tower is 15 itl and the 
radius of the base hexagon is 6.5 m. The tower was considered as a 
space gear with rigid connection rods between spherical joints for the 
kinematic investigation (see Figure 4). The actual joint at the tower 
base allows a rotation of the tower around the z axis (yaw) and pitch 
motions of the tower in the plane defined by it and the z axis. This 
means that the rotor axis will always be in the direction of the wind. 

A ring bearing is located at the tower tip which can rotate around the 
mast, to which ropes are attached. However, neither is important for 
the kinematic analysis. For a suitable selection of the rope geometry, 
the mast tip can carry out a circular trajectory around the z axis and 
the rope will remain tight. Since the tower should also adjust in the 
direction of the wind for low wind speed (shear forces), the trajectory 
of the mast tip must be as close as possible to a horizontal circle. 

With consideration of the required rotor free wheeling condition, the 
rope geometry was optimized so that the deviation from the 'circular tra- 
jectory with a radius of 2.18 m was only ±0.1/5, which, corresponds to a 
waviness of the trajector;^ of ±0.32 mm in the vertical direction. When 
installing the tower on terrain which is Inclined by 15°, by slightly 
modifying the rope length, it is possible to maintain the deviation from 
the circular trajectory with a radius of 2,l8 m to within ±0.12^, which 
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corresponds to a waviness of the trajectory of ±0.39 mm in the verti- 
cal direction. 

When the anchor point was set at 11 cm, the trajectory curve was 
again smoothed by extending the base rope by 7 •5 cm so that the devia- 
tion from an average radius of 2.19 m was only ±0.55^ or ±1.77 mm in 
the vertical direction. 

All of the cases mentioned can be looked upon as uncritical, that 
is, the tower tip runs along a circular trajectory which is perfectly 
horizontal for practical applications. 

However, this is different whfen a rope fails. For example, if a 
circumferential rope falls, then the tower has a preferred position 
provided that no other ropes are torn. The preferred position is oppo- 
site the position of the fractured rope. 

We still have to Investigate the dynamic properties of the tower. 

3 . Computer control functional model 

Under 3 a computer control functional model of a horizontal model 
axis wind generation facility (first without a tower) is designed and 
built which later on will be used in open air tests. The hub will 
allow the simulation and testing of design concepts, control concepts 
and regulation concepts without construction changes because for the 
optimization task a large number of parameters have to be varied. This, 
can be achieved by the fact that the blades can move arbitrarily and 
the force control can be achieved as an arbitrary function of measured 
deformation and load values. /193 

Deflection, flapping and <mgle of attack are varied independently 
using hydraulic cylinders. Measurement transducers for angle position, 
blade root moments and accelerations allow the valuation of the opera- 
tional state and, therefore, optimum control with a minicomputer. The 
minicomputer integrated in the rotating system provides nominal values 
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for the analog hydraulic control system. 


The supply of the minicomputer with programs and test parameters 
is done from the main computer using an optical transmission system 
without disturbances. Operational data of the rotor are transmitted 
to the same system to the main computer for processing. 

In the laboratory the hub is driven by a controllable hydraulic 
motor which is also used for power removal as a generator considering 
phase II of the project. 

The work on the computer control functional model can be divided 
into three points. 

3.1 Mechanics and hydraulics 

3.2 Process data processing 

3.3 Communication between the stationary and rotating system 

3»1 Mechanics and hydraulics 

Figure 5 gives the structure of the functional modelj the indivi- 
dual components are shown in Figure 6 , and the complete hub is given in 
Figure 10 . 

The rotor shaft is suspended with two ballbearings in a frame. 

The hydraulic disk brake, the mechanical drive for the gear and the 
flapping arms with the deflection arms are attached to it.- The entire 
gondola is suspended with four rubber bearings to the test band. The 
rotor is driven by the hydraulic motor in the present phase 1 through 
the gear. During phase II, it -can be used as a pump with a power up- 
take of up to 45 kW. The supply of the hydraulic flapping, the deflec- 
tion and angle of attack variation in the rotating system is done with 
an axial hydraulic coupling. 

3.2 Process data processing 
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The diagram of the facility including connections to the external 
computer are shown in Figure 7. Figure 8 gives the data flow diagram 
for the rotor. 

Two types of signals have to be distinguished. The measurement 
data are only transmitted from the rotating system to the stationary 
system. The control information is exchanged between the minicomputer 
LSI 11/23 and the main computer PDF 11/70. The transmission system is 
a further development of the system given in [11]. 

The analog measurement signals are partially amplified and given 
over to the data collection system. There they are digitalized and are 
transmitted bit-parallel through optical transmission systems to an 
XB terminal (see [12]). From here, using tbfe the X-bus system> they 
are transmitted to the main computer or an interactive graphic system. 
Here they can be stored on a disk or in the main computer, for example, 
for Controlling the facility or for data reduction and they can be 
graphically represented. 

The measurement data can be also transmitted to the minicomputer 
LSI 11/23 in digital form and can be used very rapidly using programs 
for controlling the hydraulic cylinder. 

The minicomputer and the main computer can communicate through the 
optoring system when the facility is rotating and the test engineer can 
intervene into programs of the minicomputer through the graphic system. 

3.3 Communication between the stationary and rotating system 

The supply with hydraulic oil is done through an axial hydraulic 
coupling, the power supply occurs through sliprlngs and rotating trans- 
formers . The data transmission is done without contact through so-called 
optorlngs (see Figure 9). This is a set of plexiglass rings for trans- 
mission away from the rotating system and into the system. The indivi- 
dual values are transmitted using infrared diodes bit-parallel. The 
transmission rate is about 30,000 words/seconda. . The measurement data 
transmission is already operating satisfactorily on the test stand 
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because due to the high scanning rate, any errors can easily be detected 
with software. However, the required reliability has not yet been /195 

achieved for communications between the minicomputer and the main com- 
puter and the graphic system. 
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FIGURE 1. BLOCK DIAGRAM FOR LINEAR ROTOR BLADE CALCULATION 
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KEY TO FIGURE 1 m PRECEDING PAGE: 

i ; 

l—aerodynamic data, wing geometry; 2 — input and change of spar data; 
3“Calculation of equivalent beam values; 4 — linear static PE calcula- 
tion of beam models; 5 — loads for design; 6--forces, stresses, dis- 
placements; 7 — eigen modes and eigen frequencies (at rest, rotating 
without aerodynamics); 8 — stability data; 9— ‘eigen modes and eigen 
frequencies with aerodynamic coefficient matrices; 10— aerodynamic 
matrices; 11— time response to dynamic loads; 12 — time dependent forces, 
stresses and displacements; 13 — aerodynamic and other time dependent 
loads 
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FIGURE 3. BLOCK DIAGRAM FOR CALCULATING THE STUTTGART TOWER 
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FIGURE 5. STRUCTURAL OF FUNCTIONAL MODEL (LEEWARD VIEW) 

































FIGURE 7. SCHEMATIC VIEW OP FACILITY 
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KEY TO FIGURE 8. ON PRECEDING PAGE; 

1~ analog data; 2— digital data; 3 — force; 4— hydraulic cylinders 
5 — force; 6 — acceleration; 7 — amplifier; 8— displacement; 9 — data 
acquisition system; 10 — ^channels; 11— hydraulic control unit; 

12 — clock; 13--channels; l4 — memory; 15 — words; l6 — rotating; 

17— optorings; l8 — opt®rings; 19 — main computer; 20— Interactive 
graphic system; 21— at rest 
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FIGURE 9. OPTORING TRANSMISSION UNIT 
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FIGURE 10. FACILITY WITHOUT BLADES ON ASSEMBLY DEVICE 
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DEVELOPMENT OP A GUST GENERATOR FOR WIND TURBIi^ES 
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Project No. ET 4372 A 
G. A. Walter, P. X. Wortmann 


1. Goals 

A gust generator is used to produce reproducible gusts in time or 
in space with variable frequency, amplitude and spectral distribution 
in order to follow the responsive wind turbine models of suitable size. 

2. Working program 

First of all, experimental investlgiUdns were performed with a 
model tunnel in order to not excessively load the drive blower in spite 
of gust excitation. Special attention has to be given to the avoidance 
of possible tunnel resonances . The development of production data and 
contract awards occurs after the Information has been collected. An 
investigation program in parallel with the manufacturing process is used 
to clarify the behavior of wind turbine models under gust loads and 
comparisons with theory are made. 

3. Status (end 1980) ; 

A model has been built with which experiments on the mechanism of 
gusc production and the reduction of the unsteady additional loads for 
the blower can be performed. The Icnoviledge obtained allowG the devel- 
opment of planning and design data for defining the project. After 
tenuous and extensive discussions with competing firms, the contract 
was let for building the gust tunnel. At the present time the final 
documents, position documents foundation documents ’.and detail plans are 
being finished for guaranteeing the scheduled beginning of construction 
in the spring of 1981. 
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In parallel with this, computer programs are being developed for 
the numerical treatment of gust problems of wind turbines, in order to 
make comparisons with the subsequent investigation programs and their 
measurement data. 

/210 

Abstract: 

The task definition for the gust generator now being built at the 
Institute for Aerodynamics and Gas Dynamics of the University of Stutt- 
gart, consisting of the gust wind tunnel and the gust grid. Includes the 
production of reproducible gusts viitb variable frequency, both in space 
and time. with variable amplitude and spectral distribution.. This is in 
order to determine the dynamic behavior of v?ind turbine models of suit- 
able size. 

A direct comparison between the measured values and results from 
the numerical treatment of gust problems then becomes possible. We 
then have an Important instrument for testing computer programs, theories 
and assumptions and we can obtain direct information and functional con- 
trol of various control concepts for the model rotor with consideration 
of all of the degrees of freedom. 

The gust tunnel is designed as open Eiffel type and is determined 
by important model law requirements as far as dimensions are concerned. 

The gusts are produced by a gust grid attached to the beginning of the 
measurement chamber, using actuation wings which can be set at different 
settings according to gust type. An external cylinder ring chamber over 
the measurement chamber allows one to generate flow around obstacles 
with a large streamline extension. 

/211 

THE GUST GENERATOR AT THE INSTITUTE FOR AERODYNAMieS AND GAS DYNAMICS 
OP THE UNIVERSITY OP STUTTGART 

1. Introduction 


Experience in the past has shown that at the present time it is 
not possible anywhere in the world to build large wind turbines which 
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are characterized by high operational safety and high lifetime as well 
as by simplicity and low Investment costs. 

If one follows the development of the research of these very 
complex machines, then one finds that often one proceeds directly to 
tests on full scale designs and this is related to the associated 
technical and financial risks. The results are not yet convincing 
because these machines are subjected not only to the enormous destruc- 
tion potential of storms, but also to the continuously changing loads 
due to turbulence and gusts of a natural atmospheric wind. This 
requires detailed research under control and defined conditions in 
suitable wind tunnels. 

For this project,- it would seem that the use of low speed tunnels 
of conventional design would be suitable. Model rotors have very exact 
requirements in terms of design size and because of their high drag 
under loads they represent a large blockage to the tunnel. 

In addition to the required measurement cross sections, substantial 
large changes are to be superimposed on the air stream in order to 
simulate typical gusts. At the present time, no such facility is avail- /212 
able . 

In the following we give a brief description of the gust generator. 

2. Task and eoals of the eust venerator 


In order to investigate the complex behavior of andividvMl systems 
and the entire system of wind turbine models which allows transfer to 
full scale, designs, a gust generator is being built at the Institute 
for Aerodynamics and Gas Dynamics of the University of Stuttgart. 


The facility consists of a gust tunnel and a gust grid and is used 
to produce various kinds of reproducible gusts. These can be specified 
in terms of a space structure and a time structure •with variable fre- 
quency, amplitude and spectral distrl’ '.tion. 



Direct comparisons of various rotor concepts and their operational 
strategies and Important comparisons between measured values and 
results from numerical analysis of gust problems become possible. 

This represents a high powered Instrument for testing computer pro- 
grams, theories and assumpblons. In contrast to calculations, it Is 
possible to consider all of the degrees of freedom of the dynamic sys- 
tem. Parameter Investigations, functional controls and testing of 
different control concepts of the rotor model can be carried out with 
minimum risks. 

mi 

3 • PeslKn of the gust generator 

The Important data about the design of the gust tunnel, for example, 
measurement chamber diameter , speed and, therefore, power requirements, 
can be specified from the Investigations of important model requirements. 
In order to have a meaningful transfer of measured results and exper- 
iences of the model motor to full scale designs, it is Important to 
consider the characteristic coefficientsv, A detailed analysis is 
found In [1], [2]. An analysis of the special case of pure flapping 
motion of rotor blades la given in [31. 

- Reynolds number 

The Reynolds number is an Important parameter for describing dyna- 
mic similarity of flows. If the same blade tip speed is maintained, 
the Re number referred to the profile chord should not be less than 
2’10*^. A 4 m model rotor should achieve values which are about 3 •10^. 

This seems to be all right as far as the flow behavior is concerned. 

- Mach number 

The Mach number at the blade tip reaches values of about 0.3* 
Tberefore, Its importance is negligible. 

- Froude number 
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This describes the influence of the gravity on the rotor. If 
other similarity laws are not to be violated, then equal Froude numbers 
of the model and the full scale version do not have to be maintained. 

- Lock number 

This parameter establishes the relationship between the moments 
applied to the rotor blade by the aerodynamic forces and centrifugal 
forces. This number in general can be maintained without great diffi- 
culty and does not give a minimum lower limit to the size of the model 
rotor. 

/2li\ 

- Aeroelastic similarity 

The aeroelastlc similarity of the full scale version and the 
model rotor requires equal Strouhal and Cauchy numbers. Maintaining 
this requirement requires a great deal of complexity, but should be 
possible with model rotors with a diameter of 4 m. 

- Free wheeling number 

The ratio of blade tip speed and wind speed should be maintained 
for a true model gust ■excitation. This then leads to the derivation of 
relatively low speeds of about 15-20 m/s. 

- Gust excitation 

The controlling quantity here is the dimensionless time t = • t. 

It should be maintained constantly. However, if the same blade tip 
speed of the model rotor and the full scale version is required, then 
for a decreasing size model this means an increase in the rotor rotation 
frequency Q and, therefore, a reduction in the gust-specific time t 
for example. For the gust grid, however, one very rapidly reaches the 
technical limit at the rotational frequency of the actuating fins. This 
implies a lower model rotor size. 
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Technloal desorlptelon of the faollliiy 


- Location 

Tbe university area at the Institute for Aerodynainlcs and Qas 
Dynamics in the area of the laminar wind tunnel already in existence 
has been selected as the place to set up the gust generator. The posi- 
tion is protected against natural wind by buildings and forests. 

- Facility components 

The following main components characterise the gust generator and 
can be seen in Figure 1 with the most important dimension: 

- measurement chamber with installations 

- gust grid 

- blower unit 

- diffuser 

The project is divided into two Individual projects in terms of 
organisation and execution. The gust tunnel consisting of the measure- 
ment chamber, blower unit and diffuser were let as a contract to a firm. 

The construction of the gust grid is being carried out within the 
facilities of the Institute for Aerodynamics and Gas Dynamics. 

t.l Measurement chamber with installations 

Air is sucked in from the free atmosphere through a noasle into the 
measurement chamber. In the measurement chamber the flow has the possi- 
bility of flowing around an obstacle with large streamline expansion, 
by evading into the ring space between the slotted measurement chamber 
inner wall and the external housing wall. The permeability of this 
internal wall can be varied widely by displacing U profiles on tracks. 
The cross sectional area of the outer cylinder ring space is almost of 
the same slge as the cross sectional ax'ea of the internal measurement 
chamber space. 


The flow is achieved which is as close as possible to reality. 

If one compares the relatively reduced power requirement, an average 
flow speed of 17 m/s is achieved. 

Protection grids are planned at the nozzle inlet and at the head 
of the lower wheel. 

The model rotor is supported on a separate foundation socket in 
order to uncouple the oscillations from the measurement chamber. 

. 2 Gust grid 

* 

The gusts are produced by a gust grid Installed at the begliiinlng 
of the measurement chamber and according to gust type, it has different 
configurations and controlled actuating wings. The flow behavior for 
the closed grid for unloaded wind turbines is shown schematically in 
Figure 1 . How the flow enters the cylindrical ring space for open gr3.d 
wing position and highly loaded model rotors is shown in Figure 2. 

The conditions can be easily described with a measurement chamber 
cross section perpendicular to the channel axis as can be seen in Figure 
3. The cross sectional area available to the flow is seen and the 
cross sectional area of the evasion ring space is about the same as 
the cross sectional area of an Internal measurement chamber space. The 
area Influenced by the gust grid is about 1/4 of the total cross section- 
al area. 

It is natural to ask the question as to which goals of artificial 
gust production are to be achieved by a gust mechanism. For present 
problems in the case of wind turbines, there are close relationships 
between aerodynamics, structural dynamics and control. The exact simu- /217 
lation of the Wide spectrum of the natural wind, its gustiness and the 
related stochastic properties are secondary for the most part. Instead, 
a reproduGlble and deterministic gust specification with defined enei’'gy 
density is very important, that is, with defined amplitude, frequency, 
structure, spectral density, etc. 


The following grid configurations are suggested for thevEealiza- 
tion of characteristic gust shape in the gust tunnel : 

- Longitudinal global and local gusts 

Radially arranged actuating wings with symmetric circular arc 
cross section which are maintained by a wing structure can cause single 
or permanent longitudinal velocity fluctuations by changing the angle 
of attack (see Figure 4). The drive of the actuated wing is done 
within the wing ring. 

- Lateral global and local gusts 

Parallel arranged actuated wings with profile cross section which 
are also held by a ring wing structure produce single or permanent vel- 
ocity fluctuations with one component -perpendicular to the tunnel 
axis by changing the angle of attack. Compare Figure 5 - The drive 
for the actuated wings is most simply installed in the middle stiffen- 
ing rib. 

Both configurations have an Increase in the blocking degree of the 
gust grid with the possibility of flow in the wing ring. At the pre- 
sent time about 50^ of the measurement chamber internal cross sectional 
area is used as the flow ring area in order to maintain the unsteady 
load due to pressure fluctuations as small as possible on the lower 
runner wheel. Also the required minimum model size must be mairstained. 

The air mass to be accelerated and delayed for this gust generator 
design has to be held to about 25^ of the air mass in the gust tunnel. / 2l8 

It is apparent -that it is very important to have good flow around 
the ring wing structure. Therefore, tests with a grid model were made 
in a water tunnel in order to clarify questions such as profile select- 
ions, profile angle of attack, maximum permissible blockage and flow 
separation in the steady case. 

Irt order to simulate the rotationally symmetric case for pure 
longitudinal global gusts, a NACA 0012 profile with 6 actuated wings 


in the opposite direction were used in the model as the wing ring 
structure. The position angle of 0® corresponded to the open position. 

A position angle of 90° amounted to practically 100 % blockage of the 
wing ring internal cross sectional area. 

The incompressible flow processes are made visible using colored 

U ^ U 

probes. The achievable Re numbers were in the range 2*10 - b'lO . 

Figures 6-9 give an impressive of the steady case for a ring wing 
position of = 0°, referred to the undisturbed flow. For angles of 
“st ^ 20°, the flow separated. 

4 . 3 Blower unit- 

The runner wheel supported on both sides is driven by a bydromotor 
without intermediate gear. The hydraulic units and the 315 kW 3-phase 
motor are Installed underneath the blower unit in the foundation socket. 
The blower running wheel has a diameter of 5.6 m and consists of 8 
individual running blades made of high quality spherical casting. 

They can be adjusted when standing still. It achieves a maximum circum- 
ferential speed of about 56 m/s. This means that the noise development 
can be looked upon as being minimum. 

Both for reasons of costs and because of testing, we decided for 
rpm control of the runner wheel using hydraulics. This design has sub- 
stantial hydraulic disadvantages. However, in addition to avoiding 
network feedback, for example, which occurs by operating with 3-pbase 
current and a rectifier, hydraulic control has the important advantage 
for testing of a higher switching possibility and a faster braking of 
the runner wheel. 
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4.4 Diffusor 

Downstream of the blower there is pressure recovery in a diffusor. ' 
For cost reasons, the diffusor length had to be limited to 8.7 m for an 
area ratio of inlet/outlet -area of 1:1.75* A rolling door is the shutoff 
at the end of the diffusor. 
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5 . 


Status of the work 


The planning, position plan, foundation plan and detail plans of 
the gust tunnel are now ready. 

The blower and the facility elements are already being manufac- 
tured. 

At the present time, contract work on the foundations is being 
advertised. We count on the beginning of assembly in September 1981 
and start up of the gust tunnel in the spring of 1982. 
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Summary about the most important technical data 


Summary about the most important technical data / 222 

- Data on the measurement chamber 

- internal diameter of the measurement chamber inner space 6300 mm 

- cross sectional area of measurement chamber inner space 31.2 m^ 

- general diameter of measurement chamber housing 8700 mm 

- cross sectional area of measurement chamber housing about 59.5 m^ 
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- measurement path length 


6500 mm 

- inlet nozzle length 


1500 mm 

- perforation of measurement chamber internal 

wall about 

15-2055 

- external diameter of gust grid 

about 

i |200 mm 

- cross sectional area of gust grid 


13.8 m'^ 

- frequency of gust grid up to about 


10 Hz 

- average flow speed for empty measurement chamber 

17 m/s 

Data on the blower and the motor 



- running wheel diameter 


5600 mm 

- hub diameter 


2660 mm 

- max. amount delivered 


530 mVs 

- max. total pressure p ® 1.2 (kg/m^) 

about 

300 Pa 

- total blower efficiency 

about 

8055 

- power at blower shaft 

about 

230 kW 

- rotation rate of blower 

about 

190 rpm 

- power of 3 -pbase motor 


315 kW 

- rpm of 3 -pbasG motor 

about 

1500 rpm 

Data on diffusor 



- inlet diameter 


5600 mm 

» net cross sectional area at inlet 

about 

19.0 m^ 

- output diameter 


6500 mm 

- cross sectional area at outlet 

about 

33-2 m^ 

- diffusor length 


8700 mm 

- geometric opening angle 

about 

6 . 11 ® 
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FIGURE 1. Intersection along the channel axis of the gust generator. The gust generator 
Is in the closed position. 
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Figure 3. Measurement chamber cross section perpendicular to channel I 
axis with viewing; direction from blower to inlet. The following: is | 
represented in sequence: Model rotor, gust grid ring wing with support#* 

perforated measuremprit chamber in a wall, external measurement chamber | 
wall 








FIGURE 5. GUST GRID WITH PARALLEL CONFIGURATION 
OF ADJUSTABLE WINGS. 
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FIGURE 6. FLOW FOR OPEN GUST GR7D. ANGLE OF ATTACK 


FIGURE 7. FLOW OVER GUST GRID RING WING WITH A 
SETTING ANGLE OF =10° 
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FIGURE 8. FLOW OVER GUST GRID RING WING FOR A 
SETTING ANGLE OF ot,,. » l8® 


FIGURE 9. FLOW SEPARATION OVER RING WING FOR AN 
ADJUSTMENT ANGLE OF a„. _ 21® 
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OPERATIONAL BEHAVIOR OP WIND ENERGY FACILITIES 
Project No. ET 4362 A 

W. Leonhard and W. Kleinkauf 

1 . Goals : 

The operational behavior of wind energy facilities (horizontal 
axis machines) of various sizes is to be investigated by using suit- 
able control methods and suitable mechanical-electrical energy con- 
verters. 

2. Working program : 

Determination of the dynamic behavior (transfer function) of the 
components of a wind energy facility — analysis of various electro- 
mechanical energy converters — development of control concept and 
Investigation of the entire facility behavior (network and island 
operation) — verification of theoretical knowledge with practical 
investigations. 

3. Status ( end 1980) : 

The principal control structure of a wind energy facility was 
developed. For simple systems with blade adjustments equipped with 
asynchronous or synchronous generators and used in a network or indi- 
vidually, we designed control concepts. Using digital simulations, 
we obtained information about the dynamic behavior. We started with 
measurement technical work for covering the operational behavior of 
small facilities (10-20 kW) . 

As a first application, we investigated the use of a slipring 
asynchronous generator. The machine works in a stand -within the three 
phase network of constant, frequency. The 'rotor is applied by a direct 
converter so that sub synchronous and oversynchronous operation are 
possible. A simple variation to this is the oversynchronous current 
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circuit where the slip power is delivered to the network through a 
direct current intermediate circuit. The static and dynamic processes 
were analyzed theoretically uslnt; digital simulation and also with a 
22-kW laboratory setup. The control used for field coordinates allows 
uncoupled specification of the effective power and the blind power. 

A 6-pulse circuit with 3-phase current output is used as direct con- 
verter. Control is done with a microcomputer. 

General remarks 


The project is being performed in conjurctlon with the Kassel 
Technical High School, Electrotechnical Division (Energy Technology 
office) and the Institute for Control Technology of the Braunschweig 
Technical University. 


- Work division between Gh-Kassel and TU Braunschweig 


The simplified structure of a wind energy facility in the figure 
shows the interface for the work. 
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The work division allows relatively great independence, espe- 
cially for the practical Investigatioi’s and the simulation calcula- 
tions, However, a great deal of cooperation is required so that common 
evaluation and presentation of knowledge is assured. 
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Gh Kassel? 

Control and dynamic behavior of wind energy facility I ; 

Table of contents 

1. Introduction 

2. Control structure of a wind energy facility 

2.1 General remarks 

2.2 Mechanical control circuit components 

2.3 Simulation of facility behavior 

3. Comments on the control of small wind energy facilities 

4. Practical work 

M.l Measurement investigations of components of wind energy 
facilities 

4.2 Control of wind energy facilities with Joint operation 

5. Other procedures 
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1. Introduction 

According to the detailed project descriptions of the contract, 
the investigations are concerned with wind energy facilities which 
have a horizontal axis design where the blade adjustment angle can be 
changed and which are intended to convert wind energy into electrical 
energy. With this restriction, we are assured that the investigations 
will concentrate on facilities now being developed and will, therefore, 
have a high degree of application. 

The several basic considerations about the topic will now be pre- 
sented for explanation. 
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Turbine 



Figure 1. Wind energy facility, energy flow and 
control . 


2 . Control structure of a wind energy facility 
2 . 1 General remarks 


When a wind energy facility Is used for producing electrical 
power, the facility function has to be tuned to the user requirements. 
Detailed knowledge about the energy flow from wind energy up to elec- 
trical supply network for individual users Is, therefore, very Impor- 
tant. Figure 1 shows this chain. 

The primary energy Is specified by the fluctuating speed of the 
air flow. Control manipulations must, therefore, assure that excessive 
component loads, excessive power and rpm variations are avoided. The 
adjustment of the wind wheel blades around their longitudinal axis (see 
Figure 2) is the most often used possibility of influencing the energy 
absorption from the wind. 

Control of energy absorption by operating the wind direction fol- 
lowing system (turning of the towerhead) can only occur very slowly 
because of blade loads (gyroscopic forces) so that in the following 
they are not used for power limitations. 
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a ... blade adjustment angle 
(a = 0® - flag position) 

^act ^^act^ supplied 

by the actuator (or torque) 

... torsion moment applied 
by aerodynamic forces, etc. 


Figure 2. Blade adjustment device 

Figure 3 gives a simplified structure of the wind energy facility 
from the control point of view. 

The entire installation is monitored by the operating control 
system. The control paths are characterized in their behavior by the 
bold face locks. 


This includes the following: 

1. Behavior of the blade (or blades) for twisting around the long- 
itudinal axis, blade adjustment path, where the actual value 

of the adjustment angle is the output. 

2. Formation of the drive moment from both of the main influenc 

ing parameters wind velocity v and adjustment angle a. 

w 

3. Influence of the moments of inertia of rotating masses (essen- 
tially the wind wheel and the generator runner) on the actual 
rpm of the generator n^^ and 

conversion of mechanical energy into electrical energy by the 
electrical generator with formation of the electrical resist- 
ance moment and the energy transportation to the network 
and users. 


2 . 2 Mechanical control circuit members 

Numerous additional influencing parameters which influence the 
control circuit of the facilluV have to be considered in a closer anal- 
ysis. The behavior of the blade and its influencing by torsional moment 
is especially complicated. The following facts have to be considered. 




ORIGINAL PAGE IS 
OF POOR QUALITY 


external Input 

r-U 


rpm 


‘wind speed 



Figure 3. Simplified structure of wind energy facility . 

- moments due to lift forces over the blade 

- moments due to bending of the blade and related mass 
displacements and lift displacements 

- moments caused by the pendulum motion of the rotor 

- propeller moments 

- moments due to blade Inertia, which for example, change 
with the blade bending (possibly additional consideration 
of the acceleration of air mass particles) 

- moments by aerodynamic damping 

- friction moments at the blade bearings 

- moments due to spring and damping properties of the blade 
adjustment facility, etc. 
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Depending on the type, these moments depend on changeable varia- 
bles such as the rpm of the wind wheel n , blade adjustment angle a as 
well as its derivatives, blade position during rotation and wind 

speed V . 

^ w 

The calculation of the lift moments from the wind speed v^ with 
consideration of the adjustment angle a, wind wheel rpm n^ and blade 
position (for example. Influence of tower shadow) requires exact 

knowledge about the power characteristic field of the wind turbine. 
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KEY TO FIGURE 4 ON PRECEDING PAGE: 

1— blade adjustment; 2 — drive moment; 3 — wind wheel gear generator 
(mechanic); ^ — network or user (island operation); 5 — electrical 
generator; 6 — operational control nominal value specifications; 

7 — controller installation; 8 — actuator for blade adjustment; 

9 — blade adjustment; 10 — blade angle adjustment; 11 — drive torque; 
12 — wind wheel hub; 13-— gear coupling; l4 — electr.; 15--network or 


user 


Lambd a , R 1 ph a 
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Lambda 

Figure 5 . Performance coefficient as a function ~of free wheeling coefficient X with 
adjustment angle a as parameter. This is a characteristic line field stored for stimu- 
lation of a two-hlade wind energy facility (10 kW range). 











Among the additional dynamic influences due to mechanical compo- 
nents which participate in the energy transfer, we have the torsional 
stiffness and the damping behavior of the element as well as the iner- 
tia of the rotating masses. 

Figure 4 gives a detail of Figure 3 and shows the structure of a 
wind energy facility which makes it possible to consider the Influences 
mentioned. Starting with this structure, we can make simplifications 
depending on the Investigation goal and the structure of the wind energy 
facility. 

Exact analysis also requires the actuator for blade adjustment 
which is the power part of control. Special investigations at the 
present time are continuing for often used units with hydraulic oil 
and cog wheel pumps, control valves, adjustable pistons, etc, 

2.3 Simulation of the facility behavior /24l 

Wind speed variations have to be specified for the digital simula- 
tion of the facility behavior. The reactions of the coupled facility 
components are calculated and printed out for the wind speeds with 
inclusion of the time dependences. The basis of the facility control 
is a control system tailored to the properties of the components and 
the purpose of the facility. The important parts of the simulation pro- 
gram include: 

- specification of the wind speed variation 

- control and monitoring system 

- facility structure (according to Figure 4) with subprograms 
for calculating the drive moment, blade behavior, etc. 

In designing the program, we were especially interested in its 
general validity. We wanted to have the possibility* of inputting the 
wind speed variations and to store extensive rotor characteristic line 
fields (Cjjj, Cp = f(X)). 
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With the software developed, it is possible to include their 
arbitrary wind speed conditions as well as directly measured condi- 
tions as excitation functions for the simulation of the facility 
behavior. Specification of the wind data can be done using a digital 
device using a satellite con^uter or directly from the wind measurement 
device through a telephone modem to the computer center (mass storage). 

In this way, we have the possibility of comparing the simulation re- 
sults later with the actual behavior of the facility. 

The most exact method at the present time for determining the 
drive moment is by using the characteristic line fields of the power 
coefficient (Cp) or the rotation rate coefficient which is usually 

used by the blade designer for quasi- stationary states and requires a 
great deal of effort (it is usually represented as a function of the 
free wheeling coefficient \ with the blade angle a as a parameter, see 
Figure 5). In addition to the high accuracy in the representation of 
the characteristic line fields, for a dynamic evaluation there has to be 
a good and rapid extrapolation in the intermediate range in order to 
avoid long computer calculation times. The evaluation of the stored 
characteristic lines was done for the simulation using a two-step method 
with Interactive transfer from one spline function and was developed 
into a linear Interpolation program which satisfies the requirements. 

In particular, critical and rapidly dropping characteristic field ranges /243 
(see Figure 5) can be included for stability investigations. 

At the present time, test runs of the previously designed simulatldn 
program are being made with a control system for island operation of 
small wind energy facilities (10 kW range). 

3. Principal aspects of the control of small wind energy facilities 

When using wind energy converters for generating electrical energy, 
three kinds of operations can be distinguished. 

- island operation 

- network operation 

- island and network operations (selectable) . 
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Figure 6. Configuration for Island operation of wind 
energy converters 
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Figure 7. Configuration for network operation of 
wind energy converters with 
generator. 


asynchronous and synchronous 


The following discuss^.^n3 refer to small wind energy facilities 
(up to 10 kW electrical power output) which are equipped with three- 
phase current generators. The use of special generators and current 
rectifier units which would allow variable rotor rpm would hardly be 
worthwhile because of the cost. 

3 . 1 Island operation 


2h6 












The wind energy facility Is not connected to the open public 
electrical supply network. For the mechanical-elecstrical energy con- 
version, a synchronous generator (SO) is suitable. The use of a asyn- 
chronous generator would mean that a controlled excitation blind power 
would have to be made available in order to achieve a constant voltage. 

The complexity for a controlled compensation device is improportion- 
ately high in the normal case. In order to check the output frequency, 
the rpm of the wind energy I'aclllty has to be controlled (through the 
blade adjustment device). If the users require approximately a con- 
stant voltage, then when thfere is insufficient wind speed, in order to 
match power and rpm, less important energy users will have to be turned 
off through a user control system according to a priority method. In 
order to assure the supply to Important users, depending on the user 
requirements, devices will have to be used which operate independent of 
the wind (Figure 6), 

/2M4 

3 • 2 Network operation 

Now the wind energy facility will operate with the public supply 
network (see Figure 7) with approximately constant voltage and fre- 
quency. Asynchronous generators (ASG) or synchronous generators (SG) 
can be used as generators. Since ^he frequency is specified by the 
network (i.e., the rpm of the synchronous generator specified and the 
rpm of the asynchronous generator can be varied by slip values), then 
in order to avoid overload the power uptake of the facility has to be 
controlled (through the blade adjustment device). In order to switch 
these facilities equipped with an asycbronous or synchronous genera- /245 
tor to the network, in addition a corresponding switching and synchro- 
nation device is required. 

In both facility types, .the power control would have to be asso- 
ciated with an rpm control system so that during startup, shutdown and 
during synchronization network failure, etc. , the control of the facil- 
ity would be maintained by it. 
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Figure 8. Configuration for network and Island operation 
(selectable) of wind energy converters. 


3 . 3 Network and Island operation (selectable) 

In the normal case the facility will operate on the network, but 
when the network falls. It Is used also for electrical power supply. 
This operational condition requires the maintenance of the conditions 
given above for network and Island operation. The use of a synchro- 
nous generator Is probably more suitable for satisfying the require- 
ments of an Individual facility (see Figure 8). 

Practical work 


^ . 1 Measurement technical Investigations of components of wind 
energy facilities 

First of all. Individual components of the wind energy converters 
such as the rotor blades and the transfer elements (for example, gears) 
are most important for the analysis. The examination of the operating 
behavior of the entire facility will be done later on. 

When selecting measurement procedures, we are especially Ir.tent on 
using those which will not require extensive facility modifications. 

The adjustable rotor blades represent the only actuator components 
In the power train of the wind energy facility. Its dynamic properties 
during the adjustment along the longitudinal axis are, therefore, very 




important. The investigation of the torsional behavior of the rotating /2M? 
blades requires a high degree of complexity. 


There with other methods the photogrammetric determination of | 

the blade behavior can be done with relatively a small amount of effort | 

without construction changes to the blades. Therefore, we planned to | 

illuminate marks on the rotor blades stroboscopically and to then mea- | 

sure the images of two geometrically separated cameras. The optical | 

I. 

evaluation gives the object coordinates in conjunction with a computer | 

facility as well as the deformation of the rotor blades. The determina- | 

I 

tion of the damping properties of the fins, etc., requires additional I 

calculations. An estimation of the accuracy and Information content of | 

I 

the measured values is only possible after measurements have been I 

carried out. Since exact analysis has shown that in particular the | 

image evaluation requires greater costs than expected, we are attempt- I 

ing to use another measurement method whose information con1:snt can be { 

better estimated. For example, we are considering model analysis. | 

I 

A structural dynamic analysis using accelerometers on the facility | 

k 

components to be investigated requires a high degree of complexity, but | 

allows a simpler evaluation of measurements. The use of a structural I 

dynamic analyzer also allows the Investigation of other facility com- I 

ponents such as the hub, gears, coupling, generator, machine house, I 

tower and guy ropes, and foundations, etc. We can also obtain inf or- | 

mation about the dynamic behavior of the entire Installation. At the | 

present time we are investigating this possibility as an alternative i 

I 

for phototometry. ! 

I 

I 

In addition to the behavior of the rotor blade, it is very import— | 

ant to know the time dependence of the drive moment in the gear and | 

the generator for other investigations. Instantaneous fluctuations 
due to wind speed changes, tower shadow Influences, etc., are espe- 
cially Important here. In collaboration with the DPVLR Stuttgart, 
measurements of the modulus facility in Schnlttlingen are being per- 
formed on the Swabian Alp. For these measurements, the facility was 
equipped with a torque measurement shaft by the DFVLR on the wind wheel i 
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side. Since this measurement method implies rather high equipment 
costs and is not easy to attach to already built and operating facil- 
ities, we are Investigating a method which can also be used for other 
facilities . 

/m 


Generator 


■planetcury geetr 


chain gear 



torque at the chain gear — modulus facility UO kV/) of the DFVLR, 


We are investigating whether a rotation angle difference mea- 
surement between the input stage and the output stage on the gears 
would allow one to calculate the torques and the transmission ratio. 
Static preliminary Investigations on the evaluation of dynamic measure- 
ments were already made on the installed transmission elements of the 
DFVLR modulus facility. As Figure 8 shows, the rotation angle on the 
chain gear was measured for static moment application. The same mea- 
surements were also carried out on the planetary gear and through both 
transmission stages. The measurements allow one to determine the gear 
stiffness and can be used for calibration. 

Additional investigations of Individual components and the exam- 
ination of control concepts will probably be performed soon using the 
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Aeroman (10 kW power class) which is an internal MAN facility. 
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4.2 Control of wind energy facilities with Joint operation 

In addition to control concepts and the Investigation of opera- 
tional behavior of individual facilities, recently we have dealt with 
problems of the control of the wind energy joint facility. Since the 
funds for the practical Investigation of a joint facility were not 
available from the project, using engineering excess funds we have 
started to build four controllable sets of machines using our own 
funds which are suitable for simulating joint operation (alternating 
current generator coupling). The generator power of the machines is 
between 5 to 15 kW. 

The Individual facilities will be controlled by a computer (in 
the form of fluctuating drive moments to the generator caused by varia- 
tinns in wind speed). The -behavior of the facility and of the joint 
facility can later on be evaluated by a fast measurement collection 
system now being built (multiprogrammer for computer coupling, meas- 
urement value Intermediate storage, control of processes, etc.). 

5 . Further procedures 

The important points of our procedure correspond to the work plan 
of the project and its steps. 

We are especially Interested in the measurement of the blade 
behavior and the final specification of the investigation method (photo- 
grammetry or modal analysis method). 

In addition to completing the program by including additional 
details (also the Investigation results determined in Braunschweig) 
for the facility simulation, we are also pursuing another path which 
Will lead to a great deal of program simplification. By comparison 
with the results of the extensive program, we are sure that there will 
be no Important loss of information. Possibly we can use practical 
results obtained from the DPVLR and the MAN facility. 
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In addition, we have started with the building of a model con- 
trolled by microprocessors which will represent the properties of 
various wind energy facilities in a lucid way by using specified facil- 
ity data. 

Abbreviations used 
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P 


M 


M 


AG 


M 


AW 


M. 


BG 


power coefficient of the wind wheel 
drive moment 

drive moment of the generator 
drive moment of the wind wheel 
acceleration moment of the generator 
acceleration moment of the wind wheel 
electrical resistance moment of the generator 
actuator moment on rotor blade 
M torsion moment on the blade 

T 

’^G* *^i generator (actual value) 


^BW 

Mel 


WG 


M 


St 


n. 




W 


w 

“i 

^B1 

\ 


nominal rpm 
wind wheel rpm 

power output of generator (actual value) 
nominal power value 
wind speed 

blade adjustment angle (nominal value) 
blade position of wind wheel 
free wheeling coefficient- 
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1 . Introduction 

Large wind energy facilities (MW rang eV) where the dynamic in- 
fluences on the operating bebavior are especially interesting, in most 
cases, are now only in the design stage. This is also true for facil- 
ities which are being developed in West Germany. 

In order to be able to examine calculation methods and theoretical 
results as early as possible, therefore, we have to apply them to exist- 
ing facilities in the power range of about 10 kW electrical power output. 

In the second report (I), in addition to general remarks and des- 
criptions of the work performed during the project, we also present 
measurement results which for the first time allow a comparison with 
theoretical results. 

2. Operational behavior of small wind energy facilities 
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2.1 Network operation 


2.1.1 Control concept 

As discussed in the first report (I) for wind energy facilities 
having asynchronous generators or synchronous generators, the power 
uptake from the wind has to be controlled or limited in order to avoid 
overloads. In order to allow startup and shutdown and switching of 
the generator to the network, additional controls are required which 
influence the rpm of the facility. 

For the most often discussed cases of the use of 3-phase asynchro-^ 
nous generators with short circuit runners, we developed the control 
concept shown in Figure 1. We are considering a facility where by 
changing the blade adjustment angle, it is possible to vary the power 
uptake from the wind. 

/253 

The power controller and rpm controller are arranged in such a 
manner that the rpm controller has a limiting effect on the output of 
the power controller, that is, the nominal value of the blade adjust- 
ment angle a . A reversal of the cont3?oller configuration would also 
s 

be possible. For normal operation, that is, for network operation, the 
nominal rpm value n is set a few per cent higher than the corresponding’ 

S' 

value of the network frequency. Except for turnon and shutdown pro- 
cesses, the rpm control system will only intervene during disturbances, 

for example, network failure when n. > n . For turning the system on 

i s 

or off and for providing network switching conditions, the nominal rpm 
value n_ can be changed (for example, turning on through a maximum value 
transducer and similar concepts). 

The variable limitations at the output of the rpm controller (see 
figure) allow intervention in the rpm control without delay for _ > 0 
and I fractions as well as unencumbered operation of the power controller. 

The dynamic behavior of a wind energy facility depends importantly 
on the type of operation. The selected configuration makes it possible 
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to tune the rpm controller to free wheeling condition and the power 
controller to network operation. 

The two control circuits mentioned have also a control of a blade 
adjustment angle a and the blade adjus'tment angle rate .^—in order to 
improve stability and dynamics. We did not have an additional loop for | 

controlling the blade angle acceleration ^9" and, therefore, the act- I 

uating moments applied to the blade because of the fact that th© fric- I 

i 

tlonal moments sometimes change discontinuously in the blade bearings. | 

i 

For small facilities, the actual value determination (determina- I 

tion of the blade adjustment angle and the corresponding rate a^) 1 

will often be too complex. Using simulations, however, it was possible | 

to show that when the two Internal control loops were dropped (the H 

power controller then directly affects the blade adjustment device) I 

there will be stronger loads on thb blades and the adjustment devices, f 

I 

but stable control behavior can be achieved nevertbeless . | 

• / 25 ^ ! 


2.1.2 Installation behavior and measurement results 


When calculating the facility reaction to various wind conditions, 
one has to determine the behavior of the blade in the longitudinal 
direction, its influence caused by the blade adjustment device, various 
moments applied to the blade as well as the formation of the drive power 
and the drive moment. 


There are no extensive theoretical methods for determining the blade 
behavior. Numerous Individual effects (aerodynamic torsion moments, 
moments due to blade bending, propeller moments, etc.) were incorporated 
in the simulation structure. At the present time we were not able to 
measure the blade dynamics because certain modifications to the converter 
are required. The preparation for carrying out a structural dynamic 
analysis was continued as mentioned in the first intermediate report. 

We plan to perform such measurements during the summer months at the 
DPVLR modulus facility in Schnittlingen. 






Figure 2 . Electrical power P. of a 10 kW wind energy facility with 
asynchronous generator for network operation (constant 
blade adjustment angle) wind speed. 

The calculation of the drive moment from the wind speed variation 
using the power coefficient/free wheeling coefficient characteristic 
fields was improved considering the fact that it could be evaluated 
on a digital computer. Figure 2 shows the critical determination of 
the output power of a 10 kW facility compared with actual measurement 
results (MAN facility AEROMAN). Since the wind speed was almost always 
below the nominal value, the power controller did not operate and the 
facility operated with a constant blade adjustment angle which was 
not completely optimum. 

In the determination of the variatlGns, we proceeded as follows: 

Using fast reacting measurement devices, we measured wind speed and 
electrical output power at the wind energy facility. This was recorded 
and stored on tape. The wind speed variation was later on specified 
through the input device to a digital computer. Using the simulation 
model, we were able to calculate the electrical output power in addi- 
tion to other parameters. Complete agreement between measured and /256 
simulation results would mean that one bad used the wind speed varia- 
tion which determines the power uptake for calculating the facility 
behavior. 

The calculated power variation is based on the wind results of one 
anemometer shown in Figure 2a. A wind speed variation of this type has 
only local applicabili '-y and can only be representative of the facility 
wind -speed under certain conditions. As would be expected, details 
occur even though there is good agreement between the calculation and 
the measurement. One can see that the calculated power variation has 
a well defined contour and, therefore, confirms the expected trend that 
the measurement with an anemometer can only follow the different wind 
speeds of a rotor surface under certain conditions. There will be better 
agreement between theory and practice if several anemometers are used 
for determining the decisive Wind speeds distributed over the rotor area 
(see Chapter 3* 3). 
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2.2 Island operation 




2.2.1 Control concept 

Based on the numerous applications, smaller facilities in the kW 
range are especially interesting for island operation. Two applica- 
tion regions have to be distinguished. 

1. Supply of simple electrical users which have no great 
requirements for constant generator frequency and voltage, 
for example, heating resistances. 

2. Supply of electrical users which require that the frequency 
and the voltage of the generator may only vary over a small 
fixed range.- 

The mechanical-electrical energy conversion is suitable for a 3- 
pbase synchronous generator. Starting with the control requirements 
discussed in the first preliminary report, we developed the control 
scheme shown in Figure 3* 

■ 7258 

An rpm controller influences the power drawn from the wind by the 
variation of the blade adjustment angle in such a manner that for suffi- 
ciently high wind speed, the rpm is maintained approximately constant. 

As already mentioned for the control of the facility with an asynchron- 
ous generator (Chapter 2.1.1), the two inner control circuits (a and a) 
are not needed if one disregards the influence of the dynamics during 
blade adjustment processes. 

For the supply of less demanding users, it makes sense to use a 
self-excited synchronous generator whose output voltage is a function 
of the rpm so that below the nominal wind speed, there will be an equil- 
ibrium between the power offered by the wind and the power generated by 
the generator due to the fact that the voltage drops with rpm. The rpm 
controller can often be made in the simple form of a centrifugal force 
controller. 
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If frequency and voltage are only allowed to vary over a small 
range, then it >is necessary to better tune the rpm controller to the 
behavior Qf the converter so that a mechanical controller usually does 
not satisfy the requirements. In order to prevent a dropping of the 
converter rpm over the partial load range, a user control which brings 
about turnoff of users of low priority has to provide for maintaining 
the delivered power below the maximum power which can be taken from 
the wind. A synchronous generator with electronically controlled vol- 
tage can be used as a generator. 

2.2.2 RPM control and measurement results 

For calculating the facility behavior, we proceeded just like in 
Chapter 2.1.2. The simulation calculations showed that when one uses 
rapidly operating blade adjustment devices (for example, hydraulic sys- 
tems with electrically controllable valves), it is possible even for 
highly fluctuating wind speeds to maintain the rpm fluctuations to 
less than ± 1 %. In order to test this informati^^n, we performed appro- 
priate measurements on the MAN converter AEROMAN (10-kW facility ) . 

They are shown in Figures ^la and 4b. Figure 4a is based on a reX^^vively 
slow working blade adjustment device with an electrical motor. 

/26t) 

Highly fluctuating and large wind speeds were present and the con- 
verter delivered the nominal load. From the measurement results, we 
find that the maximum rpm and frequency fluctuations are ±4^. 

In the measurements of Figure 4b on the other hand, we used a rapidly 
operating electromagnet ically controlled hydraulic valve. The converter 
operates at wind speed below the nominal value under the- free wheeling 
conditions which is a very difficult operating state for the control 
system. As a maximum rpm fluctuation, we find that ±1^* This value 
was not exceeded even for higher wind speeds which fluctuated even more 
and for higher loads. If one considers that we did not use any opti- 
mized controllers and actuators in the facility, then as a result of the 
investigations we find that the rpm fluctuations can be maintained small- 
er than ±1^ even for smaller wind energy facilities. Therefore, it is 
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possible to satisfy the requirements of demanding electrical users 
with regard to voltage and frequency fluctuations. 

3. Practical work 

In addition to the measurements made, the activity also included 
a unit for determining wind data, building of a microprocessor — 
controlled operating control system and the building of laboratory 
devices for investigating a Joint facility. 

3.1 Wind measurement 

For testing calculations of the facility behavior, accurate mea- 
surements of the Instantaneous wind speed are required. Almost no 
indications for this were found in the literature. We developed a 
shell cross anemometer which operates faster than commercial units and 
in addition to the wind speed and its average values, can also be used 
tn measure the wind acceleration. The measured data can be used for /26l 
processing and transferred directly to tape. In order to determine 
wind conditions typical for a wind energy facility, we started with 
measurement evaluations where the results were taken from several 
measurement points. Figure 5 shows the measured value v^^ and anemo- 
meters which were set up outside on a 10 meter high mast separated by 
2 meters, and their connection line was perpendicular to the incoming 
wind. In spite of the small distances, there are wind speed, fluctua- 
tions of more than 2 m/s as can be seen from the differences (v^^ ” ^2^ 
which were determined with the computer. For a larger separation of 
the measurement points (for example, for a distribution of the entire 
rotor area) correspondingly larger deviations should be expected. In 
order to obtain wind speed variations which determine the reaction of 
the wind converter, the values of several measurement points had to be 
processed. A simplified method without wel^ting of the wind speed effect 
(p A, v^) would, for example, be forming the geometric average (see Fig- 
ure, (v^+V 2 )/ 2 ). 
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Figure ^a . Electrically driven Figure 4b . Electrically con- 

blade adjustment nominal load trolled hydraulic valve In the 

blade adjustment direction. 

Free running 

Figure 4 . RPM control of a 10 kW wind energy facility with synchro- 
nous generator for Island operation. Excerpt from direct measurement 
value recording (multichannel recorder). wind speed, n converter 

rpm, f generator frequency. 




Figure 5. Wind speeds of 2 adjacent measurement points. 
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Figure 6. Laboratory facility for Investigating a wind energy 
converter complex. GM DC current machine (15 kW drive) 

DM 3-phase machine (10 kV; generator) 
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power part with measurement value collection 



















The wind speed determination will be very Important In the other 
measurement investigations of dynamic processes in wind energy facil- 
ities. 

3*2 Model of a wind energy facility 

Within the development of a facility model for demonstration pur- 
poses, we started on an operational system which would have a wide 
validity (both in terms of facility siae and operating range) . At the 
present time, we are working on a microprocessor — controlled unit, 
which in particular satisfies the requirements of large facilities 
operating on a network. In the next step we plan to extend this to 
several operating conditions. 

3 • 3 Joint operation of facilities 


The further continuation of the Joint operation system discussed 
in the first report (1), financed by university funds, was continued. 
Out of the block diagram of the planned overall structure shown in 
Figure 6, up to the present three motor generator sets were produced 
including the controlling regulation systems. Fart of the measurement 
determination and measurement processing system are also available and 
at the present time are being complemented by interface units* 
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1. Comparison of various generator types for a wind power plant 

For the selection of a suitable wind power plant generator, the 
following decision criteria have to be mentioned: 


1. Operation for Joint operation or Island operation 

2. Power range 

This first report Is restricted to high powered facilities (>'100 
which are to supply a Joint network. The operator of the public power 
supply company (EVU) will require a supply with the constant effective 
power and idle power from a power generator source. In this case, it 
is necessary to store the fluctuating wind energy supply In a rotating 
mass. This Is only possible when the rotor rpm of the generator f 2 
is not directly related to the network frequency f_ (see Figure 1.1). 


wind 

speed 








gear 


Turbine 


trans- 

former 


blind 


mech • 
elect, 
converter 

Figure 1.1 Power fraction of a wind energy facility 


network 

(user) 


Synchronous power plants usually used for power production do not 
have this rpm flexibility. Therefore, either a variable gear or a 
dynamically high quality blade adjustment device is required In order 
to use this which will compensate for a variable; wind energy supply. 
For Gonstant power delivery to the network, primarily a synchronous 
generator should be considered, if the stator frequency f *2 is uncoup- 
led from the network frequency f^ to a direct current intermediate 
circuit. This concept has advantages for parallel operation of 
several wind turbines which together supply a direct current collect- 
ion bus. By forming the sum, the local power fluctuations are smooth 
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network 



Figure 1.2 Wind energy facility with asynchronous generator and 
slip power recovery 

in addition. However, this configuration does not allow the produc- 
tion of capacitive idle power for direct switching of the synchronous 
machine and the network. As an alternative to the method above, we 
can consider a slip ring runner asynchronous machine which is oper- 
ated in an oversynchronous current rectification cascctde circuit 
(Figure 1.2). For the oversynchronous operation, the rotor slip power 
is supplied through a direct current intermediate circuit to the network. 
Compared with the synchronous machine with the direct current inter- 
riedlate circuit, this variation becomes more favorable the smaller the 
maximum rpm deflection is, because the rectifier only has to be designed 
for the slip power. The network feedbacks are also smaller because 
only the idle power is smooth through the current rectifier. This 
facility also requires condensers for producing idle current. 

In order to eliminate this disadvantage, it is possible to use a 
direct converter in the rotor circuit instead of the direct current 
intermediate circuit (Figure 1.3) with which an independent rotory 
field is superimposed onto the rotor. In the case of load fluctuations, 
it then becomes possible to deviate into oversynchronous and under- 
synchronous operations and capacitive idle current can be produced. 
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These advantages, however, must be contrasted with the relatively 
great complexity for the direct converter and the more complicated 
control. 

2. Double supply asynchronous machine as a wind power plant generator 

2 . 1 Principle of the double supply ASM 

Compared with the synchronous machine, the double supply ASM 
has the advantage that the flow direction in the rotor can be speci- 
fied by means of a rotory field system which is imposed from the out- 
side. For the synchronous condition, in spite of different load 
states, the machine is driven at a constant pole wheel angle using a 
control system. In addition, there is the possibility to influence 
the idle power requirement of the machine for a constant effective 
power delivery by means of rotation of the flow direction. Therefore, 
it is possible to have uncoupled effective power control and idle power 
control . 


If the rotor is supplied with a low frequency rotory current, 
then the machine will operate in the oversynchronous or undersynchro- 
nous condition depending on the direction of the rotary field. Also 
an energy flow through the rotor in both directions is related to this 
(Figure 2.1). The rotor power is calculated to be = s • 


network 



network 



Figure 2.1 
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This means that in spite of a variable generator rpm, it is 
possible to maintain the delivered stator power constant. The runner 
flow has to be controlled in terms of frequency and phase by the rpm 
and angular position of the rotor so that the rotor and the stator 
rotary field are always synchronous. The machine has an operating 
behavior similar to an asynchronous machine. 

2.2 Theoretical discussion of the double supply asynchronous 
machine 


We will briefly outline a control method which will lead to a 
far reaching uncoupling of the magnetic processes along the longitudi- 
nal and transverse axis. The mathematical model of an asynchronous /269 
machine consists of two complex differential equations and two equa- 
tions for the mechanical motion. A simplification of the relationships 
can be achieved if one assumes that the rotor currents controlled 
through the current rectifiers are assumed to be superimposed. This 
means that the voltage equation of the rotor circuit has no influence 
on the dynamics of the machine control. After conversion to the same 
wndlng number in the stator and the rotor, we have the following equa- 
tions ; 


io + L- 
-S S dt 


dt 





Wei 


+ ni 


m 



Lh 

n in 




( 2 ) 




For the rotor current assumed to be superimposed, we assumed that 

/aa 

* (4) 


-iR * ^ "" IIrI ® 




In the following we will assume a two phase rotor winding. Using 
abbreviation 


U; “ (1 ♦ Oj) • 
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equation (1) can be written in the following form: 

»s 


•Is*^ ■«[(> < »s) * in - Hn ■ 


( 6 ) 


The expression in the square brackets is defined as an extended mag- 
netization vector 


-^is 


= (1 * Oj) Is + e’<^ = |1^,| e 


.3y 


du 

- — = 0 ) 
dt ms 


(7) 

( 8 ) 


and this is used for eliminating the stator current. Using' Equation 
(2.6), we find that 
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3 777^ -Ml - J . (10) 

By multiplication with e”'^^ equation (9) is then transformed into a 
moving coordinate system which is defined by the magnetization current 
vector i . After introducing the polar coordinates for i' i d 
and U_jj, we then find the following equations by splitting the above 
into real part and Imaginary parts 


Rs dt 




* tins' - - “n + 


+ JIrI cos(eH‘y)» 


( 11 ) 




ms '-ms 
S 


ii 3^ 

t»Lmc * 


Rj — — “n ♦ 


+ Ur I cos(e+5-u) . 


( 12 ) 


”e) • - I f'-— liis'ljR! sin(«5-u) 


(13) 


Figure 2.2 shows the various angular relationships, e + & - . 

U = <S corresponds to the load angle of the machine. The quantities 
[i_p|co5 6 and |i_p| sin 6 incorporate the longitudinal and transverse 
components of the rotor current vector referred to the expanded 
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magnetization current vector i_jjjg* The longitudinal component 
controls the magnitude of the magnetization current. The transverse 
component controls the angular speed of the magnetization current 
vector and the torque. 


If these two rotor current components are changed independent 
of one another, then we have an independent uncoupled control for the 
idle power and the effective power of the stator or for the voltage 
and the torque of the machine. 


After several transformations and normalization with the nominal 
values of the machine, we finally obtain the following equation system 


di 


ms 


■S dt- * - 0 ♦ "s) 


1 j 1 + a 

_• 


[f^ ‘^Nd I *^0 ^Rdj* 


The stator time constant is 

Lc 




Rs 


i-x 

•■s “o ■ 


In addition, we have 


= 3 COS 
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2 Ro «d 
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(U) 

(15) 


(16) 


(17) 


Here is the normalized synchronous reactance 
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and is the ratio of the stator short circuit current for nominal 
excitation divided by the nominal current. 

These relationships are shown in Figure 2.3 in the form of a 
block circuit diagram. The normalized input quantities of this mach- 
ine model are the two rotor currents, the network voltage and the 
drive torque. The nominal current values which vary sinusoidally 
for oversynchronous and undersynchronous operation have to be made 
available by the controlling control circuit. The quantities re- 
quired for control are available by measurement for a doubly supplied /272 
machined^* w, c)or can be easily derived (^njs* • m)» 

cosfe ^ w)). 

A microcomputer can be used for this task. 

Figure 2.3 gives an example how the field oriented control by 
a moment controller which affects i^^ would be designed and a voltage 
controller which involves ip^j* 

The^ limiting dynamics and the restricted actuator path of the 
rectifier require that the rotor voltage induced for s ^ 0 be switched 
in to improve the dynamics of current control which can be looked upon 
as a perturbation quantity. Figure 2.4 shows this and also we indi- 
cate how the residual delay of current control can be reduced by a 
linear delay unit in the nominal value input. 

2.3 Design of the power part of the rotor supply 

For supplying the rotor of the double supply asynchronous machine, 
we have selected a direct rectifier circuit according to Figure 2.5. 

This circuit consists of three antiparallel current rectifier bridges 
which are connected with a network through a transformer. The current 
rectifier bridges are switched in a star. The direct converter can 
allow power flow in both directions, depending on the phase sequence 
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(over or under synchronous operation). The star point of the rotor 
winding is not connected. 


The transformer has the task of reducing the network voltage to 
the maximum required rotor voltage for control. It also prevents the 
formation of equalization currents between the three current rectifier 
bridges. 

The three phases of the direct converter are commercial ^ quadrant 
current rectifier component groups which consist of a power part, a 
monitoring part and a nominal value determination device as well as a 
control logic system. The control of the individual current rectifier 
bridges in the present design is done with a microcomputer which from 
the actual and nominal values of the enter control circuit calculates 
three actuator signals which are displaced by 120®. The current rec- 
tifier bridges and the transformer are selected with power reserve in 
order to provide flexible testing. In this way, we wish to bring 
about that even in the most unfavorable operating case, we could exam- 
ine how a moment shock of twice the nominal value could be controlled 
with a slip of 1555. In order to be able to investigate slip values of 
> 1555 and not to limit the controller dynamics by actuator limits, a 
transformer output voltage was designed for 21 times the rotor direct 
current for the synchronous condition (the factor of 21 was taken from 
a simulation). In order to determine the requireu limiting values, 
an actuator transformer is switched in ahead of the primary winding 
with which the maximum voltage can be continuously brought down. 

2.4 Realization of a controller concept using microcomputers 

2.4.1 Hardware summary 

At the Institute we have developed a computer card with a micro- 
processor Intel 8085 according to Pigure 2.7 for controlling of the 
double supply ASK (Figure 2.7). This design allows simple operation 
in monitoring of all control circuits. In order to have a high scann- 
ing frequency, in the design of the computer system we have made sure 
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that time intensive operations will be carried out if possible through 
hardware. For example, the multiplications are performed by a fast 
8x8 bit hardware multiplier and the phase deductions and phase divi- 
sions are realized by simple analog circuits. 

The actual value determination of the rotor current and the 
stator current is done through Hall generators whose signals are 
amplified in order to obtain a larger perturbation voltage separation. 

For the position and rpm determination, an angle decoder was attached 
to the motor axis which allows a maximum resolution of 0.1°. These 
actual values are prepared according to the requirements of the con- 
trol algorithm, filtered and are supplied to the computer through I/O 
ports. The interface between the analog part and the digital part is 
selected so that the minimum number of quantities have to be directed 
through the relatively time consuming AD/DA converter. 

: i 

In order to allow extensive valuation of the dynamic behavior of 
the control circuit and its variables, the quantities of lnter«=^st can 
be plotted if required. An analog output has been built for this. 

/274 

2.4.2 Software summary 

Shortly a program will have been finished which has the structure 
shown in Figure 2.8. The microprocessor takes over the calculation of 
the actual values for the electrical moment and the voltage. The cal- 
culation of the induced rotor voltage (perturbation voltage) and its 
switching as well as the transformation between the individual coordi- 
nate systems is also done in this way. The PI controller required 
for control can be performed, by software with simple control algorithms. 

The entire programming sequence will be determined by the required min- 
imum scanning rates of the individual partial control circuit and 
great calculation accuracy is required. The nominal rotor voltages 
in the a,b system are the output quantities which are directly given to 
the phase splitting device and the following direct converter. 
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Figure 3.1 Wind energy facility with oversynchronous 
current rectifier cascade 

3 . Oversynchronous current rectification cascade 

As a generator a slip ring runner ASM is used which operates 

only in the oversynchronous rpm range > 0 . The slip power ^’^mech 

"o. 

is supplied through a direct current intermediate circuit to the net- 
work, that is, in the rotor circuit it is possible to have only a 
single power direction in contrast to the facility discussed in section 
2. The runner rotary field can be Influenced in terms of current ampli- 
tude (by controlling the Intermediate circuit current 1 through the 

o 

rectifier). The phase and frequency, however, are determined by the 
stator rotary field and the rotor rpm, that is, only the effective power 
can be controlled. The machine as well as the rectifier take idle power 
from the network inductively which may be controllable by means of con- 
trollable condensers. The current rectifier complexity is small (6 
diodes, 6 thyristors) and a circuit of the runner circuit to resistances 
is possible as overload and over rpm protection. 

The rpm range is specified so that the rotating masses are capable 
for the maximum deviation from the nominal point to give off or take up 
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Figure 3.2 Characteristic line field of the over- 
synchronous current rectifier cascade 


the same kinetic energy. For the normal operating range, the facility 
will, therefore, be along characteristic i of Figure 3.2. If the 
limiting rpm is achieved, then the power delivered to the network has 
to be changed (characteristics 2 and 2*). For facilities with a larger 
slip, it may be necessary to hold the entire power ^rotor^ 

constant. Then it would be necessary to reduce the moment according 
to characteristic 3 with increasing rpm. This makes the stability con- 
ditions more difficult and this will be tested by means of a digital 
simulation program. 

■! 31 § 

. Present status and furthei rocedure 

For the double supply asynchronous machine, a computer program 
has been developed which considers the unsteady operation of the current 
rectifiers and simulates the dynamic behavior of the control generator. 

The results will be compared with practical laboratory tests in order 
to allow generalization. The computer structure shown in Figure 2.7 
with peripherals is now finished and at the present time is being tested. 
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The first testing will be done with a small 1.7 KW- ASM and will be 
paid for from: our own funds and this is supplied with our own trans- 
istor rectifier in the rotor circuit .i This intermediate stage has been 
introduced in order to overcome the long delivery times for the direct 
converter. Therefore, it becohies possible to build up the operating 
program in the following report year and to test it in steps. We plan 
to transfer to the larger machine (22 kW) and there will be additional 
problems because of the direct converter which can be isolated by pro- 
ceeding in this way. 

For oversynchronous current rectifier cascades, at the present 
time a computer program is being designed which will be usOd to invest- 
igate the dynamic behavior in parallel with the test set up. For the 
double supply ASM as well as for the current rectifier cascades, this 
same machine will be used^ therefore, we have a direct means of com- 
paring both methods. 

Instead of the wind turbine, a direct current machine is used for 
the drive for which a moment variation will be implemented corresponding 
to that of the turbine. First test functions (jump, ramps) can be 
specified. Later on the torque variations determined by 'GHK will be 
used as nominal values. 
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FIGURE 2.7 Hardware summary: Control of a double supply ASM with 

microcomputer 8o35 
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Block diagram of a microcomputer generator control system 
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ANALYSIS OF VARIOUS ENERGY CONVERTERS 

1. Introduction 

In the first report we investigated generators and their controls 
and compared them. The theory and laboratory setup of the double 
supply rotary field machine were discussed in detail. 

In this report. (II), we first discuss the idlgital simulation 
results of the double supply rotary field machine. After this we dis- 
cuss the measurement protocols obtained with the laboratory setup. In 
addition, the report has in it extensive description of the theory and 
the simulation results for oversyncbr'Onous current rectification cas- 
cades. 

2 . Simulation of a double supply three phase generator 

The differential equations of the double supply three phase machine 
and the block diagram of the entire control concept can be taken from 
the first intermediate report. Pgiure 2.1 gives a simplified circuit 
diagram. 

2.1 Simulation of the closed control circuit 


For the best possible representation of the dynamic properties of 
the generator, the differential equations of the three phase machine 
and of the network are related with the properties of the direct recti- 
fier. 

In order to limit the calculations, the following simplifications 
are Introduced: 

1. The output voltage of the direct reetlfler consists of segments 
of a sinusoidal function, depending on the control angle a, and the 
commutation processes are ignored. 

2 . The blockage time for zero crossing of the current is not 
considered. 
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network 










The entire block diagram of the controller concept is shown in 
Figure 2.3 of the first intermediate report. 

2.2 Simulation of the wind energy supply and the mass inertia 
of the generator unit 

In order to be able to approximately simulate the operation of a /287 
real wind power plant, the dynamics of the wind and the mass Inertia 
of the real wind energy converter have to be simulated for the labor- 
atory machine. This is done by a corresponding control of the direct 
current motor which drives the three phase machine and in this way 
gives the fluctuating energy supply of the wind to the generator. 

Figure 2.1 shows a simplification of the control of the direct current 
machine. The moments of inertia of the real facility can be simulated 

by a suitable selection of the integration time constant T . 

m 

2.3 Design of the controller 

All of the controllers (see Figure 2.3 of the 1st intermediate 
report ) are determined for the synchronous idle running condition with 
the assumption of complete uncoupling. Since the paths are of higher 
order, the theoretical design can only give indications by means of 
replacement time constants. The final dimensioning is done experi- 
mentally during a simulation. 

3. Simulation results ^ 


Figures 3 T la and 3vlb show the simulation results calculated from 
the test facility and normalized for the nominal quantities and the 
control circuit was subjected to several disturbances. 

3 • 1 Jump response of the generator for synchronous conditions 


It is assumed that the generator is in the synchronous idling con- 
dition for t < 0, i*e., the rotor current and the rotor voltages are 0. 
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At the time t * 0 the nominal quantities of the generator are intro- 
duced as nominal quantities instantaneously in the form of a Jump to 
the voltage controller and the moment controller. At the same time, 
the electrical moment of the generator is specified in the same magni- 
tude. The rpm of the machine remains approximately constant. 

Diagrams a, b in Figure 3* la show the development of the two-phase 
rotor currents. The real three-^phase rotor currents are shown in dia- 
gram c. Two of the three-phase rotor voltages are taken from d, e. 

The electrical moment m^j^ and voltage = ^nora “ ^act have 

achieved a stationary state after about 100 ms (diagrams f, g). Dia- 
gram h shows the slip s of the machine . The dasb and dot line gives 
the slip of the energy converter. The difference between the electrical 
moment and the moment of the wind wheel which results during the trans- 
ient process of the controller circuit is not sufficient to accelerate 
the machine noticeably. This is not completely achieved with the labor- 
atory facility (see diagram h, curve but the dynamic behavior of 

the simulated wind energy converter is quite close to the real facility. 

/290 

3.2 Braking of the generator from the synchronous condition to 
minimum rpm 

The generator is braked with 25 times the nominal moment within 
175 ms to its minimum rpnn Diagrams f and g show that these perturba- 
tion quantities are quickly controlled by the controllers and that the 
moment and voltage controls are uncoupled. By means of the control 
system a slip-frequency sinusoidal rotary current is superimposed on 
the rotor (c). 

3.3 Voltage actual value and nominal value Jump for asynchronous 
operation 

The generator is in the under synchronous condition (s = 15/S). For 
t > 450 ms the network voltage and the nominal value of the voltage con- 
troller are reduced by 10J5. Diagram g shows that the stator voltage 
has reached the new nominal value after about 70 ms. A comparison 
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Figure 3 -la: Dynamic behavior of a double supply three-phase 

machine (simulation) 
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Figure 3 ‘lb: Continuation of Figure 3 •la 


between the electrical moment and the stator voltage shows that the 
two control circuits have the desired uncoupling for this perturbation 
case. 
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A detailed description of the microcomputer system and the power 
part was given in the first report. In the meantime, the individual 
components have been operated and have been switched together accord- 
ing to the structure shown in Figure 2.1. The structure and operation 
of the direct converter resulted in some difficulty because of the fact 
that the star point was not connected. 

An assembler program was built for controlling the double supply 
three field machine and was tested with a laboratory setup. A scanning 
time of 2.5 ms can be achieved for the two current controllers and a 
scanning time of 5 ms can be achieved for the moment and flow controllers. I 

The first results with this controlled three phase machine were 
presented in October at the INTERKAMA 80 fair. Figure 4.1 shows the ; 

fair stand of the Institute for Control Techniques during the building \ 

phase. The figure shows the 22 kW machine, the controller unit and on /291j 
the left side there is the rectifier console for the direct current | 

machine. The direct rectifier cannot be seen on the figure. The poster 
(see Figure 4.2) in the background gives a summary about the purpose, 
theory and control method. 

: 

!■ 

4.2 Measurement results ! 

i 

I 

In the following we will give a brief discussion of measurement ; 

protocols of the test set up. Since the facility is not yet optimum, | 

these are only intermediate results. For example, there are still diff- | 

icultles because of the fact that the star point of the rotor windings \ 

is not acceptable and because of the switching delays of the partial I 

bridges of the direct rectifier. In addition, the voltage switching I 

, I- 

of the Induced rotor voltage was not yet operable. However, we already /294 I 

can see that the results obtained with the digital simulation will be i 

I 

confirmed in practice. 

Figure 4.3 shows the variation of the normalized two-phase rotor 
current of the machine. The controlled direct current machine drives 
the three phase generator and transfers it from oversynchronous condi- | 

tion to the undersynchronous condition. The phase position of the rotor | 
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KEY TO FIGURE 4.2 ON PRECEDING PAGE; 

1-^replacement diagram for superimposed rotor current; 2— current 
vectors, longitudinal component of rotor current '..a-uii* flux, 
transverse component of rotor current 'n*'*****- ; 5— behavior of 

the generator for jump-like moment and voltage changes (simulation); 

7 — the power delivery of a wind power plant depends greatly on vrind 
speed. Short time fluctuations can be equalized by using rotating 
masses (wind wheel, gearing, rotor) as fly wheels. This implies that 
the generator operating on a fixed network can operate with variable 
rpm. In this project, the operating behavior of various generators 
is Investigated for wind power plants. The figures show th$ method 
adapted for Growlan of a double supply 3-phase machine; 6 — uncoupled 
effective and idle power control of a double supply 3-phase generator; 

8 — properties of a double supply 3-phase generator. -Over and under 
synchronous operation -Uncoupled effective and idle power control 
-capacitive or inductive operation; 9 — simulation of wind v/heel; 

10 — controlled 3-phase generator; 13 — -laboratory set up. In the 
laboratory the generator is driven by a controlled direct current 
machine instead of the wind wheel.- Stored moment variations are used 
as controlling variables of the direct current drive. This means gusts 
can be simulated. Supported by BMPT. l4— wind power plant with double 
supply 3-phase machine ; 15—line ; 16 — control; 17— actual value detarmination 
18 — direct rectifier; 19 — gear; 20- — wind wheel; 21— subsynchronous opera- 
tion; 22— over-synchronous operation. 



Figure Measurement stand of the double supply 3-phase machine 

Interkama 80 


currents changes during transition through the synchronous condition. 

Figure 4.4 shows the variation of the magnitude of the magnetiza- 
tion vector Ijjjg for constant rpm of the generator referred to the mag- 
netization current (for an Idling machine). At the time t^^ the 
nominal value Is changed In jump-llke fashion for the magnetization 
current (flow controller). For the maximum voltage of the direct con- 
verter selected here, the actual value adjusts to the new nominal value 
after about 120 ms. The corresponding variation of the a component of 
the two-phase rotor current Is shown In Figure 4.4b. 

5. Further procedure ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 

As the protocol shows, the switching of the direct converter and 
the discrete operation of the microcomputer lead to higher harmonics. 
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10 t/lsl 


Figure ^.3 Variations of 2-phase rotor currents for 
transfer from over synchronous to under synchronous 
operation. 



0 2,5 5 7,5 10 t/ls) 


Figure 4.4a. Magnitude of magnetization vector for a 
Jump-like change in the nominal value. 
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Figure a — component of the 2-phase rotor current 


t/(sj 


In order to further reduce the Influence of the scanning raster of 
the microcomputer on the higher harmonics, at the present time we are 
building a new controller unit with a much higher performance l6 bit 
microcomputer. This means that the phase transformations are also 
performed digitally. The scanning times expected for the controlled 
circuit are between 1 ms and 1.5 ms. In addition, we will Investigate 
whether there will be an improved switching and control of the direct 
converter with regard to the higher harmonics. 

6. Oversynchronous current rectifier cascade 

6 . 1 Method of operatio n 

Figure 6.1 shows the circuit diagram of a current rectifier cas- 
cade. The rotor current and, therefore, the rotor power of the slip 
ring runner machine which Is connected to the network can be controlled 
through the control angle of the rectifier on the network side. Other 
details about the advantages and disadvantages of this concept are con- 
tained in the first report (I). 
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Figure 6.1. Wind energy facility with rectifier cascade 


The asynchronous generator for the oversynchronous current rec- 
tification cascades circuit is simulated with consideration of the 
differential equations of the machine and of the current intermediate 
circuit. The results obtained will be compared in another part of the 
work with measured results of a test set up and they will be interpre- 
ted. 


Calculations 


stationary 


)eration 


The nominal slip of the current rectifier cascade is specified so 


that the rotating masses a**e capable for permissible rpm deviations, 
300 
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to take up or five off the same kinetic energy. For the simulation 
we set s_ 4 ^ = 0 and s„ » 0 .I 5 . Therefore, we find s * -0.283* 

Later on we will discuss how the minimum slip cannot exceed a limiting 
value because of the expected commutation difficulties. 

Since we already have simulations of a double supply asynchronous 
machine, in order to have a simple comparison, these machines data are 
used for the cascade as well for comparing the systems. 

7 . Simulation of the oversynchronous current rectification cascade 

The differential equations of the asynchronous machine have to be 
related to those of the direct current intermediate circuit for the 
mathematical description of the oversynchronous current rectification 
cascade. In the following we will outline the procedure. 

7 • 1 Differential equations of the direct current Intermediate 
circuit 

If the voltage of the rectifier is assiuned as an imposed voltage, 
then depending on the state of the rectifier, we can write down the 
following simple differential equations for the direct current inter- 
mediate circuit (Figures 7*1 and 7*2). 



For the simulation the method of a direct coupling in conjunction with 
a cyclical exchange is used for the fixed pulse rectifier. The poss- 
ible 12 states of the rectifier, 6 commutation and 6 noncommutation 
phases, are related to the equations of the rotor with the correct 
Initial conditions and in the correct sequence. One then obtains a 
set of 12 differential equation systems. 

7.2 Differential equations of the rotor circuit 

The two phase description of the electrical quantities is selected 
for describing the asynchronous machine in the a, b system moving with 
the rotor. 
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Figure 7.1. Replacement diagram of the current rectifier 
cascade for establishing the 12 differential equation systems 
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Prom Figure 7.1» we can read off the equations of the rotor 
voltages. 


df, 


R1 


" “R1 


d^Rp 

Rr iR2(t) ♦ = Up2 

d^Rs 

•*R ’R3^^^ nr * “R3 


(7.2) 
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If the rotor fluxes are expressed by means of the inductivity of the 
machine and the corresponding transformed electrical quantities, then 
the rotor voltages can be calculated as follows after a few conversions: 


“ri 


’Ra 


P’Ral 

dt 


’Sa 


"*R1 


“Sa 

“R2 

>3 

“ A 


+ B 

"3t" 

+ C 

.’Sb 

+ r,^ 

’R2 

_’R3 

+ D 

“Sb 


The elements of the matrices A, B, C and D consist of the inductions 
of the machine and the numerical factors resulting from the transfer- • 
mations . 



7 • 3 Relationship between the differential equations of the 
rotor and those of the intermediate circuit 


Any state of the direct current bridge can be represented as a 
combination of the commutation current ij^ and the intermediate circuit 
i . The following connection relationships can be specified 

o 


R1 


’R2 


R3 


* S 


(7.4) 
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(7.5) 


■<Ra 


V 


= V 




s. 


T 

S, S and V are switching matrices which depend on 
rectifier. By combining equations (7-1) to (7.6), 
12 equation systems of the form 


(7.6) 

the state of the 
we finally obtain 



A V 



- B V 




C 


'Sa 




- 0 


*^Sa 

“Sb 

W' 


(7.7) 


7-^ Stator equation and mechanical motion equation 

The differential equations of the stator currents and of the elec 
trical torque are obtained if one expresses the rotor currents using 
1 and 1^ In the known equations for the asynchronous machine. 

8. Dimensioning of the current controller 

A simplified block diagram can be designed of the current recti- 
fier cascade for dimensioning the current controller. The following 
structure can be given for the closed control circuit: 
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Figure 8.1. 
cascade. 
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Simplified block diagram of the current rectifier 


The A quantities have the following values: 




■ ^Gfj 

*“UR 

“ “WR 

” “gn 

“el 

If 



m 

-«N 


“o * ® 


The factors T^ and V are determined using the stimulation models. The 

factor -s^ considers the fact that for idling the nomlnail: rpm is 

to/o) = 1 - s^. 

' o o 


Using this block diagram the PI controller of the current control /301 
circuit can be designed according to the known dimensioning specifica- 
tions. 
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9. Simulation results 


i 



The first simulation results are available. The numerical Inte- 
gration of the differential equations Is done using the Runge-Kutta 
method. The rectifier voltage u^ consists of sections of a sine func- 
tion, depending on the control angle a. The controller Is dimensioned 
according to Chapter 8. 

9.1 Operation of the cascade at the nominal point 


Figures 9.1a and 9.1b show the most Important electrical and mech- 
anical variables of the cascade at the nominal point. The Intermediate 
circuit current 1 and the moment m contain the upper harmonics of the 
6 X rotor frequency as well as the hlgber harmonics of the 6 x network 
frequency. They are produced by the uncontrolled rectifier and the 
rectifier In the rotor circuit controlled by the network. 


These voltage jumps caused by the rectifier also affect the commu- | 
tatlon of the rectifier (see Figures 9.1a, 9.1b). At the time t^, | 
the switching condition of the bridge Is satisfied but the voltage | 
time area Is not sufficient In order to build up a sufficiently large | 
commutation current. The current can only be commutated to the new | 
phase after t 2 » I 


9.2 


Jump response for a change In the nominal current value 1^^ 


We consider simulating a jump In the nominal current value from 

» , » 

1.13 to 0.63 Iqjj* Figures 9.2a and 9.2b give the result. The 
actual current value adjusts to the new nominal value with a slight 
undershot . According to the power division 

^Rotor = ” ® ^Stator (9.1) 

the stator power becomes smaller. 

9 . 3 Jump response for a change In the supplied mechanical power 
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Starting at the nominal point, the generator Is accelerated up to 
maximum rpm with an assumed torque within 0.1 sec . Figures 9.3a and 
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Figure 9.1a. Nominal point of the cascade 
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9.3b show that the rotor current and, therefore, the moment of the 
machine can be maintained constant by control. For uniform power deli- 
very into the network, it is necessary to reduce the nominal value of 
the electrical moment of the machine according to a function which 
depends on rpm. 


10 . Further procedure 

The simulation results discussed above allow a first evaluation 
of the oversynchronous current rectifier cascade. In the second step, 
the icf-thematical model will be Improved and especially the higher har- 
monic spectrum and the operating points important for commutation 
will be investigated. In parallel with the theoretical results, we 
will build up the test stand and bring about microprocessor control. 


Abbreviations used 
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s 

E 

0) 

m_ 


m 


m 


el 


U. 


g 

■g 


u 

V. 

p 

u. 

R. 

R. 


R 


R 


WR 

R 


slip 

mechanical rotation angle 
rpm 

mechanical moment 
electrical moment 

ohmic resistance of the direct current intermediate circuit 
inductivity 

intermediate circuit direct voltage 

intermediate circuit direct current 

rotor current 

stator current 

rotdr voltage 

stator voltage 

magnetic rotor flux 

electrical power 

commutation current 

rectifier voltage 

rotor resistance 

stator resistance 
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mechanical moment of inertia 
mechanical time constant 
magnetization current 
idling magnetization current 


CONCENTRATION OF WIND ENERGY IN VORTEX FIELDS AND THEIR /3H 

EXPLOITATION FOR PRODUCING ENERGY (PHASE I) 

Project No. ET 4252 A* 

H. Staufenbiel and H. dry 


1 . Purpose /312 

Investigation of the physical foundations of wind energy conver- 
sion which exploits the concentration of wind energy within vortex 
fields. According to measurements with wind tunnel and water tunnel 
models accompanied by theoretical analysis, an economical system is to 
be defined and the design of an outside model facility is planned. 

2. Working program 

- Development of the physical foundations of suitable vortices: 
production, structure and stability for disturbances. 

- Estimation of the wave power of various designs for the exploit- 
ation of vortex systems as wind energy concentrators. 

- Variation of the parameters for the optimization of the aerody- 
namics with the purpose of producing high energy vortices. 

- Critical analysis of concepts with consideration of costs of 
construction and problems of structural realization. 

- Concept selection and definition of wind tunnel models. 

- Investigation of selected concepts for wind tunnel testing and 
calculation methods. 

- Economy considerations compared with conventional wind ■mills. 


The project was introduced to the participants of the Seminar for 
Information purposes because it is not directly related to the 
Seminar. In this report, therefore, brief summaries are given of the 
Institutes for Aviation and Space Flight and for Plight Construction 
as well as two photographs. A more extensive final report will be 
published in the middle of 1982 (BMFT PB series). During Phase II, 
a 10 kW-WINK prototype will be tested. 

R. Windheim 


315 


3. Status (end 1980); 


The work on the topic wind energy concentration during the report 
year were designed to investigate the properties of various delta sur- 
faces with realistic installation conditions of the surfaces. 

The natural solution for a delta surface on a tower has substan- 
tial aerodynamic and structural disadvantages. Support walls below the 
wing which are more favorable from the structural point of view, could 
even increase the vortex intensity for some developments. Regarding 
the use of a concentrator under natural terrain conditions with over 
speeds, such as would be found over cliffs, for example, seems practical 
by using a perpendicular half delta surface and seems to be economic. 

The power brake was improved in accuracy and the measured power 
coefficients show an increase by a factor of 10 in the power delivery 
compared to a rotor with unreinforced Incident flow. In addition, the 
Influence of a wind direction change is demonstrated. Up to a slip 
angle of 6°, no tracking is required. When the azimuth actuator motor 
has a tracking speed of 6°/min, then no bursting of the vortices is to 
be feared. 

The work on the YEN vortex tower demonstrates the functional capa- 
bility of the principle. First measurements give a power intensifica- 
tion by a factor of 6. Similar values were also achieved with delta 
surfaces. The power Increase with reduced turbine diameter as given 
by YEN does not occur and towers with large height— diameter ratios 
"block” the outlet . 

/313 

The Institute for Aerodynamics and Space Flight (ILR) and the Inst- 
itute for Light Construction of the RWTH AAchen since 1978 have been 
developing concepts for a wind power generation facility using new 
technology sponsored by BMFT and under the direction of the KPA Juelich. 

This is a so-called wind energy concentrator (wink) . 

The small density of the natural wind energy makes it natural to 
corisider a concentration for technical and economic reasons. The 
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concentration can be done using natural terrain forms and suitable wind 
direction surfaces on designs. One special possibility of concentrating 
the energy is the production of vortices at suitable, delta — shaped 
guide surfaces. 



an angular momentum flow, then 
each blade element requires 
substantially greater twists 
than usual as the velocity 
triangle shows. 


A vortex is created at the leading 
edges of the delta wing. 

In the vortex the speeds are substan- 
ially increased compared with the wind 
speed 


The speeds in the vortex have approx- 
imately distribution over the radius. 

A substantial throughput Increase 
results in a high rotsttlon energy and 
a concentration effect. 
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An optimum distribution 
of the blade adjustment 
angle Is shown above, 

1. Rotor hub; 2. Blade 


For the same design power as a convent- / 
ional horizontal turbine 

- the wink requires a much smaller rotor 

- the wink will start at smaller wind speeds 

- the wink can control the vortex intensity 
by adjusting the angle of attack and 
therefore equalize wind fluctuations 


tip. 

Therefore, the wink results in a greater yearly energy yield. 

With consideration of the power characteristic line shown below and 
the frequency distribution of wind speeds prevailing at the location 
List/Sylt, at a 100 m hub height we find the following sum curve of 
energy to be gained per year. The reference speed is the measured wind 
speed at 12 meter height. 




1 — concentrator power plant 

(wink); 2 — conventional 1 — year; 2 — concentrator power plant 

horizontal turbine; 3— design; 3 — conventional horizontal turbine; 

4 — hub height; 5 — power char- 4— design; 5 — yearly energy gain at 

acteristic of a wink location list/sylt 


Of PO* - Q’ - r-T' 


A WINK COULD LOOK LIKE THIS 



a high power version (5-10 MW) 
on a stand 

1 — bearing and generator housing; 

2 — deflecting arm; 3 — vortex axis; 
4 — vortex generator; 5 — tower 



Structurally simpler version with 
half wing, perpendicular near the 
ground (<5MW, see title figure) 
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rr.odel-visuallzatlon v;ith smoke 


Vortex structure of half 
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Wind energy concentration in the wind tunnel. The stable vortex 
made visible by condensed air humidity in the center of the picture 
is produced by the black delta structural part and drives the small 
white turbine on the left (arrow). A power brake which can be 
aligned and has 5° of freedom is coupled with the rotor in one unit. 
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CALCULATION OP LARGE ROTOR BLADES PROM THE CALCULATION 
.■" ENGINEERING I POINT OP VIEW 

# K Hit 

John H, Argyrls , Kurt A. Braun , Warner Lang 

SUMMARY 

Extensive calculations are required for designing rotor blades of 
large wind energy facilities. In the present report we will discuss a 
part of finite element calculations. We are discussing the cross sec- 
tional data for rotor blades which are produced automatically for 
specified profile contour, spar geometry, materials and wall thickness 
distribution, and are produced for further processing In an PE system. 

The Input and change of typical parameters can be done Interactively 
by the engineer on a screen. An Iterative matching of the structure 
Is, therefore, easily possible and can be done rapidly which saves 
cost and time for the design of the blades. 
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SCHEMATIC SEQUENCE OF AN FE-CALCULATION 


Before the finite element calculation, there Is an aerodynamic 
design of the rotor blades. With consideration of the wind conditions 
for a blade of given length, profile type, chord distribution, tblckness 
distribution and twist are specified. This aerodynamlcally specified 
wind geometry Is available to the engineer In order to house the sup- 
porting components In a rotor blade, for example. In a box spar. The 
geometry of the supporting components and the materials have bo ge speci- 
fied, The Input of this Information Is done Interactively on the screen 
for various wing sections. From these spar data, the equivalent beam 
values are calculated which are required as Input for the following PE 
calculation. This complex and difficult computation step was performed 
up to the present by the engineer and the beam cross-sections were Input 
for the various blade- positions Into the PE program. 


University of Stuttgart, Institute for Statics and Dynamics of Avia- 
tion and Space Plight Designs 

**Statlcs and Dynamics Research Association mbH, Stuttgart 
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Using static PE analysis of the rotor blade idealized as a beam 
model, it is possible to determine forces, stresses and displacements 
due to various design loads (for example, eigen weight). The calcula- 
tion engineer changes his spar data until the reaction of the structure 
are satisfactory. As a rule, this makes it necessary to change the 
spar very frequently and determine the equivalent beam values . The 
automation of the last step connected with a substantial facilitation 
of the input and change process of spar data means a substantial savings 
in time and cost. 

If the static design of the blade has been Included satisfactorily, 
then a dynamic investigation is also necessary which gives stability 
Information and time responses of the structure. Here again, as a 
rule, it will be necessary to change the spar and structural data sev- 
eral times until satisfactory results are achieved. 

To interpret the PE results, it is necessary for the beam ideali- 
zation used to recalculate the generalized beam forces into stresses 
in the spar. In order to unload the engineer, this is also done auto- /319 
matically. Using clear graphical representations prepared by the com- 
puter of the stresses and displacements, a fast interpretation of the 
results becomes possible. 

In the following, we will only discuss the input and change of 
spar data. 

INPUT OP SPAR DATA 

In the Idealization of a rotor blade as a beam model using finite 
elements, for every node point and end point of the individual beam 
element, one requires the characteristic geometric data. In the pre- 
viously used Investigations we use the beam element BECOS and BETACX 
of the PE system ASKA / 1 / which are the following: 

- transverse area 

- moment of Inertia 
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- St. Venant torsion constant 

- reciprocal shear surfaces 

- coordinates of the thrust center 
(see Table 1 and Figures 2 and 3) 

In addition, the mass and inertia variables are determined for 
the dynamic analysis: 

- coordinates of the mass c.g. 

- mass of the blade cross section 

- mass inertia around c.g* 

and the latter two quantities are referred to the longitudinal dimen- 
sion of the blade (for example, dm / dr). 

The cross section of the rotor blade is idealized as a multicell 
closed shell. It can have walls, ridges and spars made of different 
materials so that a composite wing made of many materials can be treated 

The program for the interactive input of spar data requests the 
following data required on the screen, for example 

- number of blade cross sections 

- are the cross sections represented graphically? 

- number of thrust steps 

- number of flanges on the upper and lower spar side 

- coordinates of the shell corner points 

- configuration of walls, thrust steps and flanges 

- material values for individual parts 

- wall thicknesses for walls and thrust steps 

For simplifying the input, the program asks whether certain special 
cases prevail, for example 

- spar coordinates defined by a known profile? 

- perpendicular thrust steps? 

- spar walls parallel to x axis? 

- constant wall thickness in the cross section" 

- constant E modulus over cross section? 

- pure flange idealization? 
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and if a special case occurs, it Jumps over the unnecessary 'Questions 
to the engineer. 

REPRESENTATION OP THE PROFILE AND THE SPAR GEOMETRY 


On request, the profile and/or spar contour can be automatically 
drawn for the blade sections. The profile geometry can be specified 
with 

- subprogram (up to the present NACA 00 series) with parameters 

- data filed with value pairs (previously PX 77-W-153) 
therefore, it is very simple to introduce new profiles. In Figure 4, 
for example, the FX 77-W-153 profile including the spar selected by the 
engineer is shown. In this case, only the positions of the perpendicu- 
lar thrust steps are given and determined that the spar corner points Ite 
along the profile contour. In Figure 5 the spar is shown alone. The 
NAGA 0052 was represented in Figure 6 including the fitting spar. 

The graphical representation of the spar in the profile contributes 
substantially to the clarity of the design and clearly. shows possible 
errors before calculations are started. 

CALSULATION OF THE BEAM CROSS SECTION DATA 

The beam cross section data mentioned above are calculated and made 
available to the FE program (in our case, ASKA) . 

REPRESENTATION OF RESULTS 

The stresses calculated from the generalized beam forces as 

the result of the PE calculation, are automatically plotted in a figure 
similar to Figure 6 and allow fast evaluation. This part is still being 
worked on. 

REFERENCES 

[1] - , ASKA UM 202, ASKA Part I, Linear Static Analysis, User’s Refer- 
ence Manual, ISD-Report No. 73, S^dttgart, 1971, Revision F, 1979. 
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unci Vtrschiebungen 

I 

PIGURE 1. Block diagram of the linear rotor blade calculation 
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KEY TO FIGURE 1 ON PRECEDING PAGE: 


1 — aerodynamic specifications wing geometry; 2 — Input and change of 
spar data; 3— calculation of equivalent beam values; 4 — linear static 
FE calculation of beam model; 5 — forces stresses displacements; 

6 — load for design; 7 — eigen modes and eigen frequencies (standing still 
rotating without aerodynamics); 8 — stability data; 9 — eigen modes and 
eigen frequencies with aerodynamic coefficient matrices; 10 — aerodynamic 
matrices; ll~time dependent forces, stresses and displacements; 12 — 
time response to dynamic loads; 13 — aerodynamic and other time dependent 
loads 


BECO 


A 


OF POOR QUALITY 



i 

Solid bcom oF constant crost-stcHon, St.Venant torsion, no shear effects in bending 

Number of nodes 2, locoted on centroid of end cross-sections 

Degrees of freedom o, v, w, # 4^ ot eoch node 

X y X 

Figure 2. The BECOS element 



Closed thin<*walted tube with lineor toper connected to eccentric nodal points, 
St . Venant torsion, sheer effects in bending 
Number of nodes 2 

Degrees of freedom u,v,w, 0 , 4> , <t> ot each node 

X Z 

Figure 3. The BETACX element';- 






(illegible) 

Figure 4. Profile PX-77-W-153 including spar 




(il legible) ISO STUTTGART 

(il 1 egibl e) - ‘ ■ 

Figure 5. Three cell spar for profile PX-77-W-153 



(illegible) ■ — 

(illegible) 

“ Figure 6. Profile NACA 0052 including spar 



WING TECHNOLOGY PROGRAM FOR LARGE 
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In the following we discuss several problems which resulted 
during the first phase of a wing technology program which obstructed 
a fast and exact optimization of large rotor blades. The individual 
problems of various importance have various influence on the rotor 
blade design. 

DEVELOPMENT OF AN OPTIMIZED ROTOR BLADE 


The main purpose of the wing technology program is the develop- 
ment of low cost rotor blades. The rotor blade can only partly be con- 
sidered as an isolated system because there are substantial feedbacks 
of the rotor blade on the entire system. Assuming equal total facility 
concepts and the same rotor concepts, for example, Growian, it is possi- 
ble to make a simplified analysis for an Isolated rotor blade. 

The optimization path of the blade leads through aerodynamics and 
structural designs up to the rotor blade design. The energy yield ob- 
tained and the structural costs together give the energy cost. If these 
reach a minimum, then we have a cost optimized blade. 

INFLUENCES ON THE ROTOR BLADE 

For the rotor blade, the effects of the profile, geometry and 
adjustment philosophy as well as the interaction of the supporting and 
covering structure combined with the ma.terials can be computed separat- 
ely. Other influences on the energy yield and blade costs are also 
important for the entire installation. 
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ORIGINAL PAGS IS 
OF POOR QUALITY 



For a limited wing technology program, only a few combinations 
can be treated out of the entire matrix of solutions. According to 
experience with Growian, it is Important to investigate new aspects of 
the profiles and adjustments in conjunction with the structure. However, 
already during a comparison of calculated data .with those of other facil- 
ities, the influence of the calculation, the location, the total facil- 
ity and the number of units will very clearly lead to distortions. At /332 
the same time various operating states can only be estimated rougblj^ as 
to their effects. The influences of contamination (increase of the 
roughness by a factor of 10 for insects) or the erosion effect on loads, 
power and lifetime are for the most part unknown. 

Also not much information is available about gust distribution 
and it is difficult to estimate the expected nonrunning times of the 
facility. 

POWER COMPARISON 

Prom the cost estimates, one cannot derive certain information 
about the expected energy costs. At the same time, when comparing 

"" ■■■■■■■■■■ .r.: 333 , ' 



KEY TO FIGURE BELOW: 


ORIGINAL PAGE |S 
OF POOR QUALITY 


1— thicknesses; 2 — twist; 3 — chord; — tape ring; b'— slenderness; 

6 — blade; 7 — tower height; 8 — rotor position; 9-^control; 10 — wind 
following; 11 — wind profile; 12 — frequency; 13“gusts; 14 — calcula- 
tion methods; 15 — profile; 16 — geometry; 17“adjustment; 18- — facility 
19 — location; 20 — calculation; 21 — aerodynamics; 22 — structure; 

23— rotor blade; 2^1 — energy yield; 25 — blade costs; 26 — energy costs; 
27 — optimised blade; 28 — supported; 29 — covered; 30 — materials; 

31— manufacturing; 32— number of pieces; 33- operations; 33a“J^ spar; 
34 — box spar; 35 — shell; 36 — Individual shells; 37 — entire shell; 

38— steel; 39 — light metal; 4o — composite; 4l— mechanical; 

42--manual; 43 — investments; 44 — delivery state; 45 — facility; 

46— total; 47 — yearly; 48 — assembly; 49 — maintenance; 
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various facility sizes and types, it should be realized that the caj.- 
culation results are only partly supported by wind distribution, tower 
shadowing, wind direction fluctuations and wind gradient and the result- 
ing power losses can lead to different results. 

/333 

Changes in the power yield can amount to the following: 
wind distribution ± 5 percent 

tower shadowing - 2 to -7 ” 

oblique incident flow -0.5 to -1 " 

and wind gradient ±3 ” 

In the superposition case, differences can occur, depending on 
the assumptions. This can affect the energy yield and the facility cost. 

COST DEVELOPMENT 

The prediction of cost development is just as difficult as the 
prediction of rotor blade manufacturing costs. If one considers the 
optimized cost estimates of hydrocarbon manufacturers with a minimum 
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price of 30 DM/kg for the year 1990, then it may be possible to devel- 
op a rotor blade whose cost will remain the same up to the year 1990, 
in spite of increase in resin, glass fiber and labor costs. Under these 
assumptions, the development of a hydrocarbon blade would have addi- 
tional advantages. 

One condition for this is the development of rational manufactur- /33^ 
ing techniques and nondestructive testing methods for composite fiber 
components. The present state of processing speeds of composite fibers 
varies widely. Per work place or machine, there are values between 1 
kg/hrs and 200 kg/hrs. GPK (glass fiber plastic) is assumed and a simple 
recalculation to other materials can be done using the density of the 
various fibers. The differences and the difficulty of processing are 
not important. From this we find that for a suitable blade design, the 
winding of rovings with a hi.gh single thread number seems to be the most 
favorable manufacturing method at the present time. This is confirmed /335 
by the American GFK rotor blades of Kaman and Hamilton Standard. 
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The conversion of these possible processing times to mass pro- 
duction of rotor blades becomes problematical. Comparisons with numbers 
for 1000 rotor blades known from the United States over 10 years could 
be made. The resulting 100 blades per year, however, would not necess- 
arily have to lead to the lowest price. For large components with high 
investment costs of the production facilities, the manufacturing costs 
can again increase with increasing numbers when the facilities are not 
completely used. At the present state of the optimization work for 
large wind rotor blades, therefore, one has to design the number of 
pieces made according to optimum manufacturing methods. At the present 
time for a light metal blade, we find an optimum piece number of 38 
pieces per year. Similar results (32 GPK blades) are known from 




Hamilton Standard. Price comparisons with other facilities can be per- 
formed but are difficult to evaluate. 
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The cost of light metal and the composite fiber version for mass 
production reached the level approximately for the MOD-2 model. This 
corresponds to the price of about 0.4 pfennig/kilowatt hour. In order 
to achieve these low costs, which are by a factor of 4 less than the 
present Growlan blade cost, there has to l^e a substantial manufacturing 
and material development . | 

i : 

i ■ 

. t-. ■; - • : 

MATERIAL AND MANUFACTURING DEVELOPMENT I 

Material development does not mean a determination of simple cbar- 
acterlsties but the development of favorable materials for wind energy 
applications regarding weathering and loads. Various developments are /337 
required for the various materials which increase the manufacturing 
speed and rbduce the manufacturing costs, and are also Intended to Ymprove 
the lifetime of the components. Here especially we are interested in 
the composite material sector, the most recent group of constructlonV 
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o LONS STORABILITY 
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O STRAND PRESSED PROFIIEl 5 
o GOOD CORROSION RESISTA CE 

• HIC2I DYNAMIC CHARACTER iST ICS 
o STRAND DRAWITC PROFIIEI! 

® SEMIHARD F.EMFITS 


materials. It Is important to use the possibilities of the material 
and to reduce overall costs. 

The influence of various concepts on the total facility have to 
be considered. Also there may be required design changes. For two 
load cases with the example of Growian, one considers the bending 
moment variations over the blade length . Then for a lifetime calcula- 
tion of a hub, for steel of LF 1.3 (nominal operation), the hub is not 
as important as for a composite fiber blade. Using the assumed mass of 
16300 kg for a steel blade and 4800 kg for the GFK/CFK blade, LF I .3 
becomes dominating because it is on the same order as 2.5 (negative 
gust). It is important for the hub design to consider the change load 
conditions. 
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For the dimensioning of these components with the load cycle num- 
8 Q 

bers between 10 and 10^, there is not much data available. This is 
also true for the type of lifetime calculation. 
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SUMMARY 


Considering the various aspects, we can consider the following 
partial areas for additional wind energy research: 

- nominal load list 

- performance calculation 

- material development 

- manufacturing development 

- testing 
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With increasing knowledge of these subjects, the further develop- 
ment of large rotor blades will become easier and will lead to lower 
energy recovery costs faster. 
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THE ANALYSIS OF THE OVERALL DYNAMICS OP QROWIAN I 
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K. Kehl* 


SUMMARY 

Computer programs for the analysis of the overall dynamic beha- 
vior of Growian I are presented. Problems and solution possibilities 
in dynamics are discussed with some examples. 

73^2 

1. Introduction 


Large wind power facilities with a horizontal rotor axis can per- 
form oscillations when standing still and especially during operation. 
Because of operation under load, static loads occur which can be ampli- 
fied by excitations caused by ground tremors, network disturbances, 
control influences, facility imbalances and wind disturbances. This 
affects the dimensioning of the facility, both for the extreme low case 
and for the creep strength. Therefore, certain oscillation deflections 
and accelerations must not be exceeded. 

Large wind power generation facilities have to be dimensioned 
for operational reasons and cost reasons. In contrast to earlier wind 
mills, substantial supporting components have to have an almost rigid 
behavior instead of an elastic behavior. Because of the design yields 
(for example, tower, rotor blades, drive train, etc.), dynamic insta- 
bility problems occui?, for example, whirl — flutter, bending/torsion 
and flapping/defleciion flutter of the rotor blades. There are also 
oscillation responses due to system disturbances. 

2. Computer programs for the analysis of the overall dynamics 

In general, a computer model for determining the overall dynamic 
behavior of large wind power generation facilities must contain both 

¥~ ^ 

MAN - Neue Technologic 

3^13 


the structural responses as well as the couplings to the control sys- 
tem with the electrical components and also the process control system 
which is present in many systems. All of the structural parts are 
modeled with an elastic-yielding method and their relevant eigen fre- 
quencies are of comparable magnitude. Such a universal model has the 
advantage that practically all of the interesting total dynamic problems 
can be treated in a design. However, one disadvantage of this is that 
the development times of such an extensive program system are very 
long and the user costs are very high because of the enormous memory 
requirement and the long calculation times. 

73^3 

In the case of Growlan I, there is a possibility of exploiting 
certain design properties in order to solve certain individual prob- 
lems. This has the consequence that some of the problems mentioned 
above are avoided. Figure 1 shows this. A criterion which Justifies 
an uncoupling of the control system from the structure are the small 
time constants which have very short control delays as a consequence. 

This means that the control technology can be considered alone as part 
of the operational control and the rotor blade is contained as a simple 
control loop. Using the program SIMBR, simulations of the operating 
behavior overtime were performed, whei*e the wind is introduced as a 
stochastic quantity and the blade adj ustment angle , rotation rate and 
power are calculated in real time as responses. 

Structural dynamic problems of the coupled rotor/tower system 
can be simplified if the lowest flapping eigen frequency is much larger 
than the first tower bending eigen frequency (with the rotor blade as 
a point mass) with respect to the deformation in the direction of the 
wind. Symmetric rotor deformations are assumed. If In addition the 
antimetric rotor loads are uncoupled by a pendulum hub from the tower, 
then the rotor can be modeled with a simple tower Idealization. On 
the other hand, the elastic tower can be modeled with a rigid rotor per- 
forming pendulum motion (program RMOD 1). 

To a good degree of approximation, the rotor modeling can be sim- 
plified to the problem of single axis bending when the principle axes 
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of inertia of the rotor blades are untwisted and the external loads 
during operation primarily act in one of the principal directions of 
inertia. A model which contains this idealization is the program RBM 
with which one can calculate dynamic rotor blade interface force magni- 
tudes and deformations due to the various disturbances. 

Mi 

Design reductions in the rotor flapping eigen frequencies have led 
to the requirement for a refined structural modeling. The rotor is 
considered to be yielding like the elastically yielding tower and the 
flapping and deflection bending are considered (program RMOD 2). 

A coupling of the structure with the control unit and the electri- 
cal network as an additional disturbing quantity allows one to simulate 
extreme operational situations, such as network collapse for extremely 
unfavorable wind conditions (program SIM) . 

1 ) 

3 . Loads of the rotor due to tower shadow excitation ■ 


The rotor load model designed for the single axis bending 'RBM) 
on the one hand consists of a finite element part with which the eigen 
oscillation variables of the conservative system are determined and is 
shown in Figure 2. Using the subsequent modal analysis, all of the 
eigen oscillation variables of the system damped by aerodynamic forces 
can be calculated. 

The Runge-Kutta method is used to integrate the equations of motion 
and provide the dynamic responses over time. By direct calculation of 
the complex frequency changes after corresponding transformation into 
the frequency range, the dynamic response spectra are determined. As 
right Size for the motlvnn equations (disturbances), one can input con- 
stant wind, the linearized boundary layer, symmeti’ic and antisymmetric 
individual gusts, tower wind shadow or wake systems, eigen weight 
excitations and centrifugal force excitations due to coning angle and 
oblique rotational axis. The model also considers the pendulum motion 

^^with the collaboration of Prof. E. Giencke - TU Berlin, 
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of the rotor which has a blade angle feedback control. As an example. 

Figure 3 shows the variations of the blade tip deflection and the 

clamping bending moment due to tower wake disturbances for the stable 
state. The figure applies for quasistationary aerodynamic forces for 
a symmetric incident flow constant and time over the rotor area. For 
the bending moment, we have two main maxima displaced by l80® of about 
the same magnitude. The first maximum is the direct response of the 
rotor blade when passing through the tower wake displaced by about 20® 
in phase. The tip is displaced by ir and this is the response of blade 

A to the tower shadow excitation of blade B. This means that in spite 

of the antimetric disturbance of the rotor from the tower shadow, we 
have essentially symmetric loads on the rotor blade (that is, 2 £2). 

In the case of the blade tip deflection^, this is not so clear which can 
be explained by the pendulum motion. 

Figure 4 shows the influence of time dependent aerodynamic forces. 
Here we show an example of a single gust which has a trapezoidal shape 
according to the sketch. For the case of quasistationary aerodynamic 
forces, there are only unessential deviations compared with Figure 3- 
This shows that the disturbance of the tower wake Impulse has already 
decayed after one revolution and that the' rise time of the previous gust 
is so large that Only unessential additional dynamic amplifications 
occur. It is interesting to consider the unsteady condition for the 
aerodynamic force trial solution (using the Theodorsen and Schwarz func- 
tions). This leads to maximum bending moments which are about 30 % 
smaller which are displaced by about 15° with respect to the stationary 
case in phase. From this it becomes clear that the unsteady aerodynamic 
forces have a substantial Influence on the dynamic responses because of 
the high frequency perturbation components. 

If the results of the tower shadow calculation with dynamic res- 
ponses from the stationary flow (static), eigen weight and centrifugal 
force excitation as well as an Idealized boundary layer are superimposed 
then for the nominal operating case one finds that the tower shadow is 
not so Important with respect to maximum loads as shown in Figure 5. 
These loads are much more Important, however, for the operating strength 


design of the facility. 


/ 3^6 


^ . Overall dynamic structural behavior 

Calculations were performed with the program RMOD 1 for analyzing 
the dynamic stability and response behavior of the coupled tower/rotor 
system. The computer program assumes an elastic tower with a lumped 
rotor mass which is described by its 3 ° of freedom. 

= 1st tower bending in wind direction 
q2 = 1st tower bending perpendicular 
q^ = 2nd tower bending in wind direction 
qij = 2nd tower bending tranverse 
q^ = 1st torsion 

as well as an elastically rigid and zero mass 2~blade rotor performing 
pendulum motions having biade angle feedback control shown in Figure 6. 

The elastically rigid machine bouse can be aligned with the wind. The 
rotor rotation axis is inclined by an input angle with respect to the 
horizontal. 

The computer model I shown in Figure 7 was discussed in [1]. Ref- 
erence [2] gives detailed results. The modal tower motion equations 
and the discrete rotor equation motions are coupled to the entire sys^ 
tem using a coordinate transformation with the unknown constraint 
forces. The stability analysis is performed for 2-blade rotors which 
perform parameter-forced oscillations using the Floquet theory. The 
State matrix is integrated using a Runge-Kutt a method of ^th order with 
Gill coefficients starting with unit initial conditions over the rotor 
circulation; the so-called transition matrix is found. The complex 
eigen values of this matrix are called characteristic multipliers and 
can be used for the stability analysis. In the program, as a reference 
point, at the same time the eigen value analysis of the state matrix is 
performed using coefficients which are averaged over 1 revolution. The 
eigen value calculation is done using the QR algorithm and the matrix 
is first transformed to the Hessenberg form. 
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If the system is dynamically stable, then response calculations 
of the steady state can be performed for various excitations, such as 
rotor imbalances, wind gradients, wind direction tracking and oblique 
flow . 


4*1 Dynamic stability calculation 

Figure 8 gives the result of a stability analysis. It can be seen 
that the nominal operational point is in a stable region. A doubling 
of the nominal speed, wliicb is permissible according to Growlan oper- 
ation, can cross over the Instability boundary for excessive rotation 
rates with zero loads. Operation with power generation is practically 
always in the stable region. One eigen oscillation mode has been found 
to be a critical degree of freedom which contains primarily tower tor- 
sion and to a great extent, the first tower bending mode in a transverse 
direction. Prom the lower half of Table 1, which contains the eigen 
vectors normalized to 1, we can see that the other tower degrees of 
freedom are only negligible in the 5th eigen vector. This is different 
for the calculation with constant coefficients. Here there are substan- 
tial couplings with all of the tower degrees of freedom which leads one 
to the suspicion that the stability property will also be much different 

Figure 9 shows this assumption. The Lehr damping measure is 
plotted over the nondimens ional rpm which is referred to the damping 
for the nominal operating state. In the representation "without addi- 
tional system" (lower curve), one can clearly see the dropoff to the 
unstable region for excessively high rpm. As a reference, we show the 
input value of the structural damping which here was selected at a value 
of 0.255 which is low for steel design (corresponding to a logarithmic 
decrement of about 0.0125). The calculation with constant coefficients 
gives results which are both quantitatively (more than a power of 10) 
and qualitatively completely different from the calculation with period- 
ic coefficients. For time variant coefficients, the system damping 
decreases with rpm,i for constant coefficient calculations, it increases 
practically linearly. 


The danger of dynamic instability led to the requirement of damp- /3^8 
Ing the critical degree of freedom with additional measures. A dynamic 
additional system which acts in a transverse direction consisting of an 
oscillating mass and a spring and damping member was designed for the 
front machine house shown in Figure 10. With this configuration, on 
the one hand the critical torsion degree of freedom in the first tower | 

transverse bending mode is improve in terms of damping properties. On I 

the other hand, the dynamic responses in the transverse direction which I 

lead to large oscillation deflections are reduced. Figure 10 gives 
results from calculations of the optimization of the additional system. 

From the restrictions above, we made a compromise which, however, is 
not represented by any local maximum. 

With a selected additional system, the damping properties of the j 

critical degree of freedom are substantially Improved (Figure 9 — upper | 

curve), so that we are quite certain that the danger of the unstab illty | 

has now been removed. I 

k . 2 Dynamic response ealoulation ] 

In addition to the deflections in the transverse direction which 
are defined for the hub point, there are also pendulum deflections of 
rotor motion which are substantial dynamic responses of the modeled 
total system. 

Figure 11 shows nondimens ional maximum transverse deflections of 

8 

the hub point plotted against the dimensionless rotation rate for the f 

steady state. Two resonance phenomena can be established. The first 
peak occurs at about 0.7 n^ahd is the result of an excitation of 
twice the rpm harmonic. Boundary layer influences occur as disturbances 
here and the oblique attitude of the rotational axis with respect to 
the horizontal. Oblique flows and wind direction tracking Increase 
these values in some cases . The second resonance peak excited by the 73^9 

first rpm harmonic is the result of rotor imbalances and is about pro- 
portional to the amount of imbalance which here is defined as an ideal- 
ized additional mass at the blade tip. For excessive rpm and incident i 


flow speeds above the nominal speed, the instabilities occur without 
the additional system which lead to ever increasing deflections for 
linear systems. 

Figure 11 shows the deflections with and without the additional 
system in the machine house. Comparison of the two curves clearly 
shows that with the damped additional system the resonance phenomena 
are reduced by almost an order of magnitude. | 

Figure 12 in the top part shows the maximum and minimum transverse 
deflections which are not the same in every operating state. One rea- 
son for this is that for the steady case, in general, no harmonic oscil- 
lation occurs but instead, there is a superposition of two harmonics | 

at and 2 n which are weighted differently depending on the load state 
(small sketches in the top of Figure 12). Especially for power genera- I 

tion (which is not shown here), there are substantial differences 
between both extreme values in the operating range between 0.85 and | 

1.3 n^. i 

In addition to the representation of the transverse deflection of 
the hub point. Figure 12 also shows the relative motion of the additional 
system. In the resonance region the magnitudes of the damper deflection j 

are about twice as high as the hub point deflection, whereas in the 
operating state range they are only about half as large. This means 
that the additional system during normal operation will practically I 

stand still with consideration of the friction forces ignored here and \ 

the designed restrictions of the dampers. This has primarily a positive | 

effect on the creep strength of the additional system. | 
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The lower figure of Figure 12 shows the phase relationships of the ■ 

extreme deflections of both motion quantities. They are displaced by 
90° in the first resonance region and by l80° in the second region. 

In the transition regions, there are intermediate phase relationships, 
that is, within the operating rpm range . 

i 

The maximum pendulum deflections of the rotor as a function of the [ 
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incident flow speed and rpm are shown in Figure 13 without considera- 
tion of friction forces. One can see that between the turnon speed 
(V = 6 ra/s) and an extreme speed of 30 m/s which can be withstood for 
a short time, there are no noticeable differences in the pendulum 
deflections. This confirms the assumption that in the case of Growian 
I, the stationary speed which acts in a constant manner only has an 
indirect influence on the pendulum motion. The reason for this is 
because a constant speed distribution does not occur over the rotor 
area. The ground boundary layer leads to a height gradient which is 
idealized here as a linear speed drop and for the two stationary wind 
speeds, they are varied in three steps as a quantity to be superimposed. 

From the computer model I we find that the perturbation speed of the 
ground boundary layer has a linear Influence on the pendulum deflections. 

The variation of the pendulum amplitudes with rpm has a hyperbolic char- 
acter which explains that when standing still, there is a neutral equil- 
ibrium state. In order to avoid the large pendulum amplitudes which 
occur for low rpm, additional pendulum dampers were installed in Growian 
I which allow a free pendulum motion over a defined angular range. 
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5. Summary 

The design features (pendulum hub, stiffness relationships, etc.) 
allow uncoupling of individual problems from the total dynamic analysis 
so that for their treatment, one can use simplified computer models 
which, however, do not have a universal character. For calculating the 
structural dynamic behavior of Growian I, two computer models are pre- 
sented and results are discussed. As for the rotor blade loads, the 
excitation by tower shadow which occurs as a perturbation of the Incident 
flow conditions for systems operating leeward Is discussed. The 
dynamics of the coupled tower/rotor system Is deserlbed In terms of 
stability and oscillation response for various operational loads. Mea- 
sures for increasing the modal damping properties and for reducing the 
transverse oscillations of the machine house are discussed. 
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KEY TO FIGURE 2 ON PRECEDING PAGE; 


1 — rotor load model; 2 — simplification: single axis bending; 
3~untwisted Inertial principle axes of rotor blades; 4 — finite 
elements calculation (beam model); 5 — n - nose; 6— 2n-degrees of 
freedom; 7— eigen frequencies, eigen nodes, genes?allzed masses; 

8— conservative system; 9 — modal condensation; lO^-m - modal degrees 
of freedom; 11 — eigen value determination vector iteratlon/QR 
algorithm; 12 — generalized damping; generalized excitation forces; 
13— modal analysis; l4 — Integration according to Runge/Kutta; 

15 — dynamic response calculation; 16 — disturbances (right sides); 
deterministic— constant wind — boundary layer— symmetric and anti- 
symmetric Individual gusts — tower shadow eigen weight centrifugal 
forces. Indeterministic; input power spectra; 17— outputs in the time 
and frequency range — tranverse forces — flapping bending moments — 
flapping deformation — tower spectra of dynamic responses for indeter- 
minate excitation — over oscillation coefficients 
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Figure Influence of unsteady aerodynamic forces on tne flapping 
Cfendtng moment at the blade clamping point due to the tower wake flow 
due to a symmetric individual gust. Wake duration T, = T /20.5 
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Figure 5. Time variation of the hub moment (x = 11 m) and blade tip deflection 
after superposition of the eigen weight loa<^ and wind load (including the tower 
wake) for nominal operation . static states 
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horizontal axis wind turbine (HEWT) 



Figure 6. Model assumptions of the overall dynamics calculation 
of large wind tunnel facilities 
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Figure 7. Computer model I (RMODl) 








nondimens ional rpm 

Figure 8. Stability map of Growian — without additional damper 
structural damping D = 0.2?^ 


nondlmenslonal damping D/D 



unstable* I ** stable 
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Phase relationship minimum deflection nondimensional maximum deflection 



Fi?ure 12. Dynamic response for 100 kg Imbalance in » blade tip 
without power delivery, v = 24 m/s wind speed. 
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Figure 13- Maximum pendulum deflection as a function of 
RPM (n/n^) — parameter: 

— parameter: ground boundary -layer ( AV^,) incident flow 

speed (V) ^ 

— power => thrust = 0 
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CONTROL OF OSCILLATIONS OP LARGE FACILITIES 736? 

GROWIAN II: Single blade mode? 

H, STKEHLOW and G. SEITZ 
Dynamic Design of the WEA-DEMO Facility 

- 1 blade rotor 

- drive system 

- gondola tower 

Oscillation excitation for wind energy facilities with horizon- /368 
tal rotor axis: 


main excitation to 
blade (rotating system) 


fi : rotor rpm frequency 

- measures which contribute to the control of oscillation /369 

problems of large wind energy facilities with a horizontal rotor axis: 

- fast runners with one or two rotor blades of hlKh stiffness 

- reduction of dynamic weight loads for stiff deflection blade 

- avoidance of aeroelastic rotor stability problems for bending 
and torsion stiff blades 

- rotor with flapping Joint and/or pendulum hub 

- reduction of flapping bending moments 

- free rotor rotation for (doubly supplied) asynchronous generator 

- reduction of dynamic deflection bending moments 

- tower tuning for overcrltlcal rotor operation 


gravity • 

shear flow ' 1R » 2G- 

tower shadow s d ^ 12)G- 

inclination of rotor axis : “Ui » 2G- 

incldent oblique flow : lii » 2J2- 


.tti 


*ET 4240 A 
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- reduction of tower load for the operating r pm range 

- traverse problem 

Dynamic features of one and two blade rotors: 

- missing rotor symmetry leads to a system with time dependent 
(periodic) parameters 

- with a central Jointed rotor beating, the influence of the 
time bearing parameters on the total system is greatly 
reduced 
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Wind 



Demonstration facility concept: 

— levered 1-blade rotor with Impact Joint 
— free rotor rotation (rotating Joint) with double 
supoly asynchronous machine 

— soft bending tower with rope suooort for* overcritlcal 
operation 
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count- erwelght for mass 
Vjfc/ equlllzatlon 



a 



Dynamic main component of the demonstration facility 
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. Eigen frequency of the system components of the DEMO wind 
installation 

- rotor blade 

- tower 
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“ 0.4 - 0.2 0.0 0.2 0.4 

Bending eigen inodes of the demonstration rotor (blade (aounterv/elght ) 



rotor blade demonstratlch 



bearing by. double beam mo 























operating 






Influence of an additional mass 
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XI. Specification of the eigen frequencies of the DEMO wind 
facility 

(collection of results) 
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For the specification of the eigen frequencies for the main com- 
ponents of the DEMO facility, the following special items should be 7380 

noted: 


- The eigen values and the eigen modes which apply for the total 
WEA system as a rule have multiple frequencies even without the influence 
of aerodynamic force. 

- These multifrequency "eigen oscillations" which are typical for 
these systems with periodical coefficients in general (without aero- 
dynamic forces) consist of three oscillation components which depend 
on one another: 

( (o-Ji ) - component 
u) - component 
(w+S2) - component 


where B: rotor rotation rate 

to: system "eigen frequency" 

- The WEA system "eigen frequencies" approximately coincide with 
the component eigen frequencies (rotor, drive, tower) as long as there 
is no critical resonance case. 


The following resonance cases have to be considered: / 38 I 

(1) Principal parameter resonance for (2v-1)0 v « 1, 2 ... 

important for stability v = 1 0,5ft resonance 

(2) Harmonic resonance with vft v » 1, 2 ... 

(second parameter resonance) 

important for stability: v * 1 1ft resonance 
important for response : v ;< -► (1-8)ft resonance 

In the evaluation of the resonance points, one has to consider /382 
that for the selected one-blade rotor concept for the tower dynamics. 
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the force transmission from the rotating system to the fixed system 
is important, whereas for the drive system (drive train), the torque 
transmission plays a role. 

We have flapping bending Wg oscillation components important 

for tower (beriding) 

deflection bending (u)^ oscillation components important 

for drive (torsion) 

(w5±S2) oscillation component important for tower 
bending (rolling) 
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3rd tower rolling 


“/ — 3rd tower pitching 
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w 2nd flapping 
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Is deflection 
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*> 1st flapping 
■ 1 

2nd tower rolllr^ 
2nd tower pitching 


1st tower torsion 


1st tower pitching 
rolling 


4 1 nom. '■) 


K) ilr» 


— ra<l/s 

Eioen frequencies in the fixed system of the demonstration facility 
(iJ.. = 4,5n4 rart/s) 


Loads for operation of the DEMO wind facility 
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thrust force distribution dr ] tanpiential force distribution 
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ir rotation (wt) 

rotation in rotor center (bl2K3e side) 





demonstration facility 






• -tower rough (c^ = o. 7) demon stration faci lity 

tower smooth (c = o, 4 ) 
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Flapping bending moment variation 








Radius r/R 






Plapping/def lection angle and elastic blade tip deflections of demonstration 
facility in operat ion 


















IV. Basic discussion of the passage problem in the DEMO facility /391 | 
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Variation of forces over rotor main bearing for the basic case 
with 10 m/s wind speed for demonstration facility 

F - vertical force 

' X 

- horizontal force 

F_ - longitudinal force in mast direction 
tower excitation) 
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■passage problem for demonstration tower: Steady and unsteady amplification functions 


INVESTIGATIONS ON THE DYNAMIC BEHAVIOR OP WIND ENERGY FACILITIES 
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» 

John H. Argyris 
Bertold Kirchgassner 

S UMMARY 


Selected results of dynamic analysis are presented for a special 
rotor blade design with flapping Joint, deflection Joint and a coupling 
between angle of attack and the flapping angle. This includes the 
stability of the rotor blade for nominal conditions and the displace- 
ments and stresses for stationary operation, for a gust condition and 
in the tower wake. In addition, continuing work in the area of the 
dynamics of the total facility are presented. 
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INTRODUCTION 

The dynamic investigation of wind energy converters Includes two 
different problems} on the one hand, self-excited oscillations, dynamic 
Instabilities and on the other hand, foreign excited oscillations, the 
dynamic response of the system to external excitation forces. Dynamic 
Intabilities, for example, flutter and ground resonance, have to be 
avoided at all costs because, due to the rapidly increasing amplitudes, 
the support structure is destroyed in a short time if large amplitudes 
are not damped by nonlinear effects. Using a stability investigation, 
we therefore have to prevent the occurrence of self-excited oscillations. 

Par and excited oscillations cannot be avoided for wind energy 
converters because even in the ideal case of constant incident flow, 
there will be at least cyclical weight forces over the rotor plane for 
horizontal axis machines. For vertical axis machines there will be at 
least cyclical aerodynamic forces. Using a favorable re latlonsbip 
between the eigen frequencies and the excitation frequencies, resonances 
have to be avoided in order to make the deformation amplitude and 

University of Stuttgart, Institute for Statics and Dynamics of 
Aviation and Space Plight Design 
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therefore the alternating loads on the structure small. 


ROTOR BLADE DESIGN WITH FLAPPING AND DEFLECTION JOINT 
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If one considers rigid body motions in a dynamic system whose 
stiffness is only the result of external forces, for example, centri- 
fugal forces, then the alternating loads which occur in the structure 
can be reduced among other things by allowing the periodic or stochas- 
tic excitation forces to produce substantial oscillations. However, 
they are allowed only to contain small elastic components and therefore 
only load the structure slightly. 

During the research project entitled "Investigation of rotor 
stressing and smoothness of operation of large-scale wind energy con- 
version systems", which was supported by lEA with the participation of 
BMFT, a rotor blade design with flapping joint and deflection joint [1] 
^[4] was investigated in which the stationary and unsteady loads are to 
be absorbed for the most part by the centrifugal stiffness. In addi- 
tion, the rotor blade bad a flapping joint on the rotation axis, a 
deflection joint at about 10^ of the blade length and a coupling 
between the blade adjustment angle and the flapping angle for stabiliza- 
tion of the cone angle, the so-called blade angle feedback control. 
Figure 1 shows the blade model and the FE Idealization. 

In addition, a mass distribution was selected for the rotor blade 
for which the flapping bending moments for the most part vanish for 
nominal operation because the resultant of the corresponding centrifugal 
force and aerodynamic force components points in the direction of the 
blade longitudinal axis. For operating points outside of the nominal 
state , of course, changed cone angles occur. 

When a horizontal axis machine rotates, a cyclical gravity force 
applies to the rotor blade and its Importance increases with increasing 
facility size. For the rotor concept considered here, in this way a 
deflection oscillation is excited which runs in a counterpbase to the 
gravity force moment because the low«:St deflection eigen frequency is 
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far below the excitation frequency. For the most part, this is a rigid 
body motion. 

DYNAMIC ANALYSIS OF THE BLADE MODEL 7398 

For the linear investigation of the blade design a number of assump- 
tions are made: 

- rotation axis clamped 

- rotation rate constant 

- quasistationary, linearized aerodynamics 

This resulted in a system of linear differential equations of 
motion with constant coefficients which was set up using the FE system 
ASKA and was condensed using the eigen modes of the system without 
aerodynamic and Coriolis forces C3]. Figure 2 shows these eigen modes 
for the blade model with free flapping motion or free deflection motion 
and blade angle feedback control. In the condensed system, the stabi- 
lity of the system was determined by calculating the complex* eigen 
frequencies and the motion of the system by direct integration. The 
important data of the blade model is given in Table 1. 

A stability analysis of the corresponding rigid body system with 
flapping, deflection and blade angle feedback control [4] gave an 
instability for the rigid body defleetlon eigen frequency. In addition, 
estimates for the gravity^ — excited deflection motion led to the predic- 
tion of excessively high deflection amplitudes and, therefore, the 
principle of free deflection was given up and a deflection joint was 
planned vjlth a viscous damping and a spring stiffness. The lowest 
deflection eigen frequency also remained under the rotation frequency. 

7399 

STABILITY OF THE ROTOR BLADE FOR NOMINAL OPERATION 

For the blade model modified in this way, the stability for nominal 
operation was investigated as a function of the blade feedback control 
factor and the ratio of angle of 'attack change and flapping angle change 
[3]. Figure 3 shows the frequencies and Figure 4 the damping values as 
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a function of blade feedback control factor. The first deflection 
eigen frequency is very slightly damped in this case as well. For a 
continuing investigation of these oscillation modes, probably non- 
linear aerodynamics would be n^hessary in order to determine the aero- 
dynamic damping more accurately. The second eigen frequency belongs 
to the lowest flapping eigen mode which for n = 0 essentially repre- 
sents a rigid body motion. For increasing end the oscillation mode has 
increasing parts of higher flapping eigen modes. At the same time, the 
damping of this flapping motion rapidly decreases above n = 4; the 
motion remains stable in the investigated range up to n - 7. 

DYNAMIC RESPONSE OP THE BLADE MODEL 

Also the response of the blade model was Investigated for nominal 
operation with cyclical gravity force excitation for which there was a 
global gust from 10 to 20 in/s within one second and with a simulated 
tower wake with a wind speed reduction of 45?^ C33. The results for 
nominal operation and tower wake are acceptable within the framework of 
the linear theory. The results for gust excitation have a more quali- 
tative character, first of all, because of the large displacements in 
the deflection direction and also because of the nonlinearities of the 
aerodynamics which occur for the substantial displacements which are 
Ignored in the linearised model discussed here. 

Figure 5 shows the displaceinents of the blade model for gust loads. 

^ is the deflection angle, 1 / is the flapping angle, Ug and u^|j are the /400 
displacements in the flapping direction at the node 8 in the center of 
the blade or at node l4 at the blade tip. Wg and w^|j are the corres- 
ponding displacements in the deflection direction. The gust only 
starts after one rotation so that up to t - 3* 77 s* we only see the 
gravity -excited oscillation of nominal operation. • The response after 
the beginning of the gust consists of the gravity-excited oseillation 
with rotation frequency and the superimposed gust response whose most 
Important component is a deflection oscillation of high amplitude with 
a slightly damped first deflection eigen frequency. Figure 6 shows the 
corresponding stresses at various points of the blade. Here the 
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gravity-excited oscillation with rotational frequency dominates, 
whereas in the outer ranges of the blade, higher frequency components 
also can be found after the gust starts. 

Figure 7 shows the displacements of the rotor blade for nominal 
operation with the tower wake and Figure 8 shows the corresponding 
stresses. There is an additional load on the outer wing because of the 
excited higher eigen modes. 

In 'Figures 9 and 10, we finally show the flapping and deflection 
angles as a function of blade angle feedback control factor n. There 
is hardly any dependence in the deflection direction, but the flapping 
angle amplitude is drastically reduced for Increase in n. 

/ 1{01 

PROGRAM SYSTEM FOR THE STABILITY ANALYSIS OF THE TOTAL SYSTEM 

One of the most complicated problems in the dynamic investigation 
of the wind energy converter is the stability analysis of the overall 
system. During the project ET A '’OPTIWA”, at the present time a 

program system is being developed toward this . The displacements of 
the partial systems tower and rotor are described by the eigen modes 
of the partial systems, the kinetic and potential energies of the entire 
system are set up and then the linearized equations of motion are 
derived using the Lagrange equation. In the general case, the motion 
equations contain cyclical coefficients and the stability analysis is 
done using the Floquet theory. The number of generalized degrees of 
freedom of the entire system is first restricted to about 30. 

NONLINEAR INVESTIGATION OF THE ” STUTTGART TOWER” 

In [I] a new tower design was suggested through v;hich the dynamic 
behavior of the entire system can be substantially Improved so that in 
this way a soft and overcrltical support of the rotor shaft can be 
brought ab out . The tower has a j olnt bearing at the lower end and the 
tower head is guided along approximately a circular trajectory using a 
corresponding rope support. The frequency of the yaw oscillation of 
the tower head is, therefore, proportional to the square root of the 
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rotor thrust and for nominal operation is far below the rotational fre- 
quency. 


Investigations on the kinematics of this tower concept started in 

[ 5 ] were continued during the project ET 4406 A "OPTIWA”. Also there 

is a nonlinear static and dynamic analysis using a nonlinear PE program 

system LARSTRAN. 
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Figure 1 . Rotor blade and finite element idealization 
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position of turning iQ’lnt at r = 

6.0 

m 

position of CG at r = • 

30,657 

m 

blade mast 

12,989^^ 

kg 

moment of inertia around turning joint 

1,0991'*'^ 

m^kg 

moment of inertia around blade longltudlna: 
axis 

1,2580^ 

m^kg 

maximum wing chord 

6,0 

m 

rpm 

1,667 

rad/s 

cone angle for nominal operation 

12 
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Table 1. The most Important data of blade model 
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Figure 2. Eigen oscillation modes of the conservative 
system without Coriolis forces for a blade model with no 
flapping, no deflection and a blade angle feedback fac- 
tor of n = 5 = 1.667 rad/s 
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Figure 3. Imaginary part of the complex eigen frequencies 
5 a .junction of blade angle feedback factor and 
1,667 rad/s 
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Figure 5. Flapping and def 
n = 5 with a gust beginning 
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